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THERMOSPHERE DYNAMICS WORKSHOP II 


Preface 


The second Goddard workshop on thermosphere dynamics was held at the Ramada Inn in near- 
by Calverton, Maryland during October 3-5, 1984. Much of the discussion at the workshop involv 
ed the measurements from the Dynamics Explorer spacecraft (DE-1 and DE-2) and ground based 
measurements at radar facilities and optical observatories. Theoretical models were presented 
which amply attest to our maturing understanding of the processes that govern the structure and 
dynamics of the Earth’s thermosphere. The phenomenology ranged from auroral meteorology to 
equatorial tides and gravity waves. Some of the discussion was devoted to the thermosphere of 
Venus. The papers, abstracts, and extended abstracts included in this volume are in alphabetical 
order by author and cover most of the material that was presented. 

We express our thanks to the individuals who chaired the sessions and to the speakers of the in- 
vited and contributed papers for the high quality of their presentations. Our special gratitude is 
due to Dianna Gemma from Birch and Davis and Aleta Johnson (GSFC, Code 610) for coor- 
dinating the meeting. 
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APPLICATIONS OF A VENUS THERMOSPHERIC CIRCULATION MODEL 


S. W. Bougher, R. E. Dickinson, E. C. Ridley, and R. G. Roble 
National Center for Atmospheric Research 
Boulder, CO 80307 


A variety of Pioneer Venus observations suggest a global scale, day-to-night Venus ther- 
mospheric circulation. The two-dimensional hydrodynamic model of Dickinson and Ridley (1977) 
correctly predicted the gross characteristics of this largely symmetric circulation. However, it failed 
to calculate the observed cold nightside temperatures, and the exact phases and densities of the 
neutral constituents. This paper presents model studies of the dynamics and energetics of the 
Venus thermosphere, in order to address new driving, mixing and cooling mechanisms for an im- 
proved model simulation. 

The adopted approach has been to re-examine the circulation by first using the previous two- 
dimensional code to quantify those physical processes which can be inferred from the Pioneer 
Venus observations. Specifically, the model was used to perform sensitivity studies to determine 
the degree to which eddy cooling, eddy or wave drag, eddy diffusion and 15 /un radiational cool- 
ing are necessary to bring the model temperature and composition fields into agreement with 
observations. Three EUV heating cases were isolated for study. 

Global temperature and composition fields in good agreement with Pioneer data were obtained. 
Large scale horizontal winds < 220 m/s were found to be consistent with the observed cold 
nightside temperatures and dayside bulges of O, CO and C0 2 . Fine tuning required that an eddy 
coefficient <20% of previous one-dimensional models be used for nightside diffusion (K = 7.5xl0 6 
cm 2 /s). Very little eddy diffusion was required for the dayside (K<4xl0 6 cm 2 /s). Observed dayside 
temperatures were obtained by using a 7-19% EUV heating efficiency profile. The enhanced 15 /mi 
cooling needed for thermal balance is obtained using the best rate coefficient (K C02 _ 0 = 5x10“’ 3 
cm 3 /s) available for atomic O collisional excitation of CO 2 (0, 1, 0). Eddy conduction was not 
found to be a viable cooling mechanism due to the weakened global circulation. The strong 15 fim 
damping and low EUV efficiency imply a very weak dependence of the general circulation to solar 
cycle variability. Finally, the NCAR terrestrial thermospheric general circulation model (TGCM) 
was adapted for Venus inputs using the above two-dimensional model parameters, to give a three- 
dimensional benchmark for future Venus modelling work. 
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Figure 3. Global mean Venus one-dimensional thermospheric temperature 
profiles: eddy cooling tests. Here K t is the turbopause value of the eddy 
coefficient and Z t corresponds to the turbopause level. 
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Figure 4. Venus two-dimensional model temperature response to eddy cooling. HSYM corresponds to the 
Keating et al. (1984) model atmosphere, and VTS3 to the Hedin et al. (1983) empirical model. 
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Figure 5. Venus two-dimensional model response to viscous drag tests. The effects 
are shown for various adjustments of a Rayleigh friction coefficient (K RA ) (a) 

Too, (b) terminator u(z) profiles. 
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Figure 6. Venus two-dimensional model response to viscous drag tests (c) composition. 




Figure 7. Family of terminator horizontal wind profiles. Various drag profile shapes and K RA 
magnitudes were examined to obtain the best day-night distribution of model composition and 
temperature fields. 
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Figure 8. One-dimensional radiative transfer model reference 
temperatures: dayside and nightside mean profiles. These were 
calculated using the enhanced 15 /tm cooling (K^ - 0 = 5xlO~ 13 cm 3 /s) 
possible via atomic Ocollisional excitation of C0 2 (0, 1,0). Results are 
given for the BEST EUY case (7-19% profile) heating. 
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Best Parameter Adjustments for Each EUV Efficiency Case 


Parameter 
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Figure 9. Best parameter adjustments for each EUV efficiency case. The 7-19% case was isolated for further study since the 15/un 
cooling K cc , 2 ^ 0 rate adopted is comparable to that inferred by Sharma and Nadile (1981) from a terrestrial rocket measurement 
of C0 2 emission limb radiance above 100 km. 




Figure 10. Venus two-dimensional model temperature: BEST EUV case. Comparisons are made with 
Pioneer Venus model atmospheres (HSYM and VTS3) as well as with the adopted nightside 
temperature profile of Massie et al. (1983). Solar minimum fluxes of Heroux and Hinteregger (1978) 
(F10.7 = 80xl0~ 22 Wnr 2 Hz- 1 ) are used to derive a new noon temperature profile. Exospheric 
temperatures are seen to vary by nearly 70 °K from solar minimum to solar maximum conditions. 




Figure 11. Heating and cooling rates at 60° and 120° SZA: BEST EUV case. Weakened adiabatic 
heating and horizontal advection on the nightside (plus improved 15 /xm cooling) combine to give 
cooler nightside temperatures than Dickinson and Ridley (1977). 


Model and Empirical Data Set Comparisons at 150 km 



Averages Over Dayside (LT = 12-16 hr) 

2-D Model 
6/84 

Hedin et al 
H 

. (1983) 
H/2D 

Keating et al. (1984) 
K K/ 2D 

H/K 

n T 

1.39(10) 

1.01 

1.36(10) 

0.99 

1.02 

1.38(10) 

n o 

4.70(9) 

1.74 

3.56(9) 

1.32 

1.32 

2.70(9) 

n co 

2.32(9) 

1.05 

1.80(9) 

0.82 

1.29 

2.20(9) 

n co 2 

5.82(9) 

0.77 

7.13(9) 

0.95 

0.82 

7.52(9) 

V 

1.01(9) 

1.11 

1.04(9) 

1.15 

0.97 

9.06(8) 

V 

6.39(6) 

1.73 

4.57(6) 

1.24 

1.40 

3. 70(6) 

T 

271.3 

0.97 

252.5 

0.90 

1.07 

280 


Averages Over Nightside (LT = 19-24 hr) 


n T 

1.43(9) 

1.15 

1.39(9) 

1.12 

1.03 

1.24(9) 

n o 

1.12(9) 

1.14 

1.03(9) 

1.05 

1.09 

9.80(8) 

n co 

1.10(8) 

1.26 

1.09(8) 

1.25 

1.01 

8.73(7) 

n co 2 

9.70(7) 

0.89 

9.90(7) 

0.91 

0.98 

1.09(8) 

\ 2 

7.63(7) 

1.17 

9.74(7) 

1.50 

0.78 

6.51(7) 

V 

1.49(7) 

0.14 

3.11(7) 

0.29 

0.48 

1.08(8) 

T 

119.3 

0.98 

129.8 

1.06 

0.92 

122.1 


Figure 12. Model and empirical data set comparisons at 150 km. 
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Figure 13. Best EUV case contour plots (SZA vs Z): 

(a) Total temperature in deg K. 

(b) Horizontal velocity in cm/sec. 
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Figure 14. Solar minimum case contour plots (SZA vs Z): 

(a) Total temperature in deg K. 

(b) Horizontal velocity in cm/sec. 
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Figure 15. 60° SZA model helium (He) calculations. It is apparent from comparison of the one-dimensional 
and hydrodynamic model predictions for He at 100 km (Z ~ -8) that the self-consistent large-scale dynamics is 
crucial to a proper calculation of composition. One-dimensional models are insufficient. 




Figure 16. Venus two-dimensional model composition, helium: BEST EUV case. Noon and 
“nightside peak” profiles of helium are shown. The two-dimensional symmetric model ap- 
proximately simulates dynamical conditions for the observed 3-4 AM helium peak near LT 
= 23. 
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LOCAL TIME 


Figure 18. Equatorial slices (SZA vs Z) for the Venus three 
dimensional adaptation of the terrestrial TGCM: 

(a) T (total temperatures) 

(b) h (height field in km) 















DYNAMICS OF THE NIGHTTIME THERMOSPHERE AT ARECIBO 


R. G. Burnside and J. C. G. Walker 
Space Physics Research Laboratory 
Department of Atmospheric and Oceanic Science 
University of Michigan 
Ann Arbor, Michigan 48109 

C. A. Tepley 
Arecibo Observatory 
P.O. Box 995 

Arecibo, Puerto Rico 00613 

Incoherent scatter radar observations of the nighttime F layer at Arecibo, Puerto Rico, are used 
to determine the 0+ diffusion velocity at different altitudes. Further analysis allows evaluation of 
the neutral wind and the ion-drag force in the direction of the magnetic meridian. The local ac- 
celeration of the meridional neutral wind is also determined. The possibility of evaluating the 
viscous force using incoherent scatter radar data has also been investigated. Preliminary results in- 
dicate that, at certain times and at altitudes above about 350 km, viscous drag may be a signifi- 
cant term in the neutral equation of motion. Combining these results allows height profiles of the 
meridional pressure gradient to be deduced. The pressure gradient thus derived is compared with 
that determined from measurements of the horizontal temperature gradient and that given by the 
MSIS model atmosphere. 
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Geophysical parameters that may be measured or 
derived from incoherent scatter (radar) or 
optical (FPI ) observations at Arecibo 


parameter 

symbol 

vertical 

variation 

horizontal 

variation 

electron density 

N 

radar 

radar 

ion temperature 

T i 

radar 

radar 

electron temperature 

T e 

radar 

radar 

neutral temperature 

T n 

MSI S 

FPI 

meridional ion velocity 

v x 

radar 

radarf 

zonal ion velocity 

v y 

radar 

- 

meridional neutral velocity 

u x 

radarf 

FPI 

zonal neutral velocity 

— u y 

- 

FPI 


t The diffusion velocity, which may be derived from the 
observed slope of the electron density profile and the 
MS1S model atmosphere, is also required to evaluate 
these parameters. 


20.62 RSI ZM = 393.5 KM 


velocity (m s' 1 ) 



acceleration (cm s' J ) 


Figure 1. (a) Parallel ion velocity (v n ), diffusion velocity 
(v d ), and inferred southward neutral wind (u x ) on 17-18 
August 1982. (b) Derived ion and viscous drag 
accelerations, and the net southward pressure gradient. 
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Meridional ion and viscous drag 


16-17 August 1982. 17-18 August 1982 



H.SJ. 


Figures 2 and 3. Examples of meridional ion and viscous drag at 350km above 

Arecibo. 


25 







acceleration (ca s~2) velocity (n/s) height (ka) 


Arecibo 


fleridionai forces at 350k* 


flrcclbo 


Meridional forces at 350 ka 


16~17 August 1982 17-18 August 1982 7-8 October 1983 S~1Q October 1983 



H.S.T. as.T. 

Figures 4 and 5. h m = height of F layer peak. u x = southward neutral wind velocity. F id = inferred ion drag. 
F v isc = viscous drag. F x = inferred net southward pressure gradient. F msis = pressure gradient predicted by 
MSIS-83 model atmosphere. 
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Figure 6. Observed meridional gradients in ion temperature from incoherent 
scatter measurements at Arecibo. 


300-350 ka 



. <b) 

\ i r~ 

Feb. 1982 

~~i — — i — 

— T 1 — 

1 

< 

u. -I.., 

i m-” 




- 

\\ 1 1 

\ 1 

\ 

- 

1 

* « < 

■ t 1 — 

i t 

. . 1 



Figure 7. Monthly mean values of meridional ion temperature gradient for four 

different seasons. 
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Figures 8 and 9. Examples of meridional temperature and density gradients at Arecibo, and the net southward 

pressure gradient predicted by the MSIS-83 model atmosphere. 





MEDIUM AND LARGE-SCALE VARIATIONS OF DYNAMO-INDUCED ELECTRIC 
FIELDS FROM AE ION DRIFT MEASUREMENTS 

W. R. Coley and J. P. McClure 
Center for Space Sciences, Physics Program 
The University of Texas at Dallas 
Richardson, Texas 75083 

Current models of the low latitude electric field are largely based on data from incoherent scat- 
ter radars. We are extending these observations through the addition of the rather extensive high 
quality electric field measurements from the Ion Drift Meter (IDM) aboard the Atmosphere Ex- 
plorer (AE) spacecraft. We present here some preliminary results obtained from the Unified 
Abstract files of satellite AE-E. This satellite was active from the end of 1975 through June 1981 
in various elliptical and circular orbits having an inclination near 20°. The resulting data can be ex- 
amined for the variation of ion drift with latitude, longitude, season, solar cycle, altitude, and 
magnetic activity. The results presented here will deal primarily with latitudinal variations of the 
drift features. 

Figure 1 shows data from a single AE-E orbit. Of particular interest is the upward ion drift 
enhancement at 2:38 UT. This corresponds to 18:30 solar local time (SLT). This enhancement of 
the upward drift after sunset prior to reversal is a classic pattern that is often seen both by 
ground-based radar and satellite. However, on many occasions this prereversal enhancement 
(PRE) is absent. In order to investigate the reasons for these and other variations of the observed 
drift features it is of course necessary to deal with large quantities of data. 

Figures 2a and 2b show the type of data presentation used in this study. Figure 2a is a “mass 
plot” of vertical drift data from AE-E obtained from the 15-second average values of the Unified 
Abstract files. All vertical drifts in the file between -20° and 20° dip latitude (DLAT) for the 
August-September 1978 time period are plotted versus SLT. In all there are some 6987 points. 
Figure 2b shows the same data averaged into 30-minute bins. The dotted line represents the 24-hr 
average vertical drift. For comparison, Figure 3 shows some incoherent scatter radar results from 
Jicamarca for the July-August 1968-1971 time period (Fejer et al., 1979). In both cases we see the 
same general features: downward flow at night reversing at 0600 SLT, a daytime peak followed by 
a decline to a prereversal enhancement between 19 and 20 SLT, and a reversal to downward flow 
at night. 

Figures 4a and 4b cover May-July 1979 from -20° to 20° DLAT. The satellite was in a near- 
circular orbit during this period within altitudes of 450-475 km. There are some 17,000 points plot- 
ted on Figure 4a, most falling in a band some 75 m/s wide. The exception is the large spread in 
velocity in the post-sunset time period. 

The next four figures take the data set of Figure 4 and break it up into 10 degree wide bins of 
dip latitude. Figures 5 through 8 cover the ranges 10° to 20° DLAT, 0° to 10° DLAT, - 10° to 
0° DLAT, and -20° to 10° DLAT, respectively. The raw data plots indicate generally a much 
smaller spread in the vertical ion velocity at any given local time. The averaged plots indicate a 
systematic variation of the average velocity with latitude. The largest average upward velocity is 
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found in the most northerly latitude bin (summer hemisphere). The most southerly bin (winter 
hemisphere) has a net downward velocity. Other features also show latitude dependence; the even- 
ing prereversal enhancement is strong in the summer hemisphere and weak or nonexistent in the 
winter hemisphere. 

Figures 9 and 10 show a similar pattern for the Nov 77 - Jan 78 time period. Figure 9 shows the 
10° to 20° DLAT bin (winter hemisphere). Figure 10 shows the -20° to - 10° DLAT bin (sum- 
mer hemisphere). As before there is a generally upward bias of the drift in the summer 
hemisphere and a downward bias in the winter hemisphere. 

One possible interpretation of this pattern is that the dynamo induced ion velocities are being 
modified by neutral wind flow, the meridional component of which flows generally from the sum- 
mer to winter hemisphere under solstice conditions. Dachev and Walker, 1982, calculated the ver- 
tical ion drift velocity imposed by dynamo electric fields and the zonal and meridional winds in 
the 19-22 SLT range (Figure 11). These calculated patterns give good qualitative agreement with 
the data presented. This model also predicts longitudinal variations in the ion drift. This has not 
yet been tested for in the AE-E data set. Such tests should provide a useful check on this propos- 
ed mechanism. The preliminary results presented here indicate that IDM data from the AE and 
the more recent Dynamics Explorer B spacecraft should continue to disclose some interesting and 
previously unobserved dynamical features of the low-latitude F-region. 
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Figure 1. Data from AE-E orbit 13406 showing the sunset behavior of N, and V v at low dip latitudes at 315 km altitude in the 
Atlantic longitude sector on 22 April 1978. The post-sunset enhancement and reversal of V y is clearly visible in this example of 
our routinely processed data. 
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Figure 4b. May - July 1979 - 20 ° to 20 ° DLAT 
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Figure 6a. May - July 1979 0° to 10° DLAT 
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Figure 7a. May - July 1979 
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Figure 7b. May - July 1979 - 10° to 0° DLAT 
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Figure 8a. May - July 1979 -20 
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Figure 10b. Nov 1977 - Jan 1978 -20° to - 10° DLAT 
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Figure 11. Contours of vertical ion drift velocity calculated for D 
months (northern hemisphere winter). A realistic model of the 
magnetic field at 300-km altitude is used for magnitude, inclina- 
tion, and declination. The electric field is assumed to be constant 
and eastward at a value of 0.6 mV nr 1 . The wind field varies 
linearly with geodetic latitude as indicated on the right-hand side 
of the figure. The dashed line denotes the geomagnetic equator, 
and the letter J refers to the Jicamarca Observatory. Regions of 
maximum upward drift velocity are shaded. Their locations cor- 
respond to the distribution of depleted plasma regions observed 
at winter solstice and shown in Figure 6. 
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Figure 12. Vertical ion drift velocity calculated for northern 
hemisphere summer (J months). The wind field is the same as 
that used in Figure 11, with appropriate seasonal changes. 
Regions of maximum upward drift are shaded; their distribution 
corresponds closely to the observations of depleted plasma 
regions shown in Figure 7. 
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ATOMIC NITROGEN DENSITIES NEAR THE POLAR CUSP 


M. J. Engebretson and J. T. Nelson 
Department of Physics 
Augsburg College 
Minneapolis, MN 55454 

The Neutral Atmosphere Composition Spectrometer (NACS) on board the Dynamics Explorer-2 
spacecraft sampled several major and minor thermospheric gases including atomic nitrogen. We 
present here a selection of passes over the polar cusp that provide a quantitative measure of the in- 
crease of N densities due to soft particle precipitation occurring in this region. Increases in N den- 
sities are frequently observed, but are smaller than accompanying increases in N 2 densities. Our 
observations support earlier studies suggesting that 1) N densities increase more rapidly than O den- 
sities during periods of high solar EUV flux and 2) N densities are larger in the summer 
hemisphere than in the winter hemisphere. A series of passes in February 1983, late in the lifetime 
of DE-2, indicated N densities at 2000 km altitude which were a factor of = 2 larger near the 
southern cusp than near the northern cusp. 
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TABLE I 


Atomic nitrogen densities in or near the polar cusp late in the 
lifetime of DE-2. 

ORBIT DATE F Kp CUSP ALTITUDE (km) DENSITY (cm" 3 ) 


Northern Hemisphere 


8543 

83047 

91.9 

4 

Y 

236 

6.4 x 10 

8553 

83048 

94.2 

3 

Y 

238 

6 x 10 6 

8562 

83048 

94.2 

3 

N 

186 

1.1 x 10 7 


Southern Hemisphere 


8436 

83040 

125.3 

3+ 

Y 

236 

1.4 

X 

10 

8514 

83045 

94.1 

3 

N 

230 

1 

X 

10 

8551 

83048 

94.2 

4- 

Y 

205 

2.4 

X 

10 

8574 

83049 

99 

3 

N 

199 

2.8 

X 

10 


52 




19 / 2/84 


DE-B MACS 


ORBIT NO. : 1630 MASS : 30 40 14 4 

PROCESSED : 7/17/83 SYMBOL *. $ A N 4 

MAF KEY : D 81325 56817325 81325 58785293 B NACS COMP2001 



10 ** 4 -| 1 1 1 1 1 — I 1 1 1 1 1 1 1 1 1 1 1 1 1 


UT (SEC) 

56817 

57317 

57817 

58317 

58817 

ALT < KM) 

437. 

324. 

304. 

388. 

541. 

LAT 

67. 

33. 

-1. 

-34. 

-67. 

LONG 

43. 

41. 

39. 

37. 

35. 

LSTCHRS) 

18.9 

18.9 

18.9 

18.9 

18.9 


Figure 1. Ion source number densities of NO, Ar, N 2 and He measured by the Neutral Atmospheric Composition 
Spectrometer (NACS) of the Dynamics Explorer-2 satellite on orbit 1630, November 21, 1981. Number densities 
are plotted on a logarithmic scale as a function of universal time (UT). Also shown on the abscissa are satellite 
altitude, geographic latitude and longitude, and local solar time. 
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Figure 3. Ambient number densities of N, O, N 2 , and He determined from DE-2 NACS data for orbit 1630, 
November 21, 1981. Number densities are plotted on a normalized logarithmic scale as a function of universal 
time. The densities shown in the graph are to be multiplied by the appropriate scale factor shown at the top of 
the figure to give the correct ambient densities. 
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Figure 6. Ambient number densities of N, O, N 2 , and He determined from DE-2 NACS data for orbit 574, 

September 11, 1981, as in Figure 3. 
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Figure 8. Ambient number densities of N, O, N 2 , and He determined from DE-2 NACS data for orbit 8551, 

February 17, 1983, as in Figure 3. 
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SIMULATION OF THE THERMOSPHERIC TIDES BY USE OF THE 
NCAR THERMOSPHERIC GENERAL CIRCULATION MODEL 

C. G. Fesen, R. E. Dickinson, and R. G. Roble 
National Center for Atmospheric Research 
P.O. Box 3000 
Boulder, CO 80307 

Numerical calculations of the thermospheric tidal winds and temperatures were produced by use 
of the NCAR Thermospheric General Circulation Model (TGCM). The present calculations were 
restricted to solar minimum equinox conditions. The effects of viscosity, conductivity, diffusion, 
ion drag, winds, and temperature gradients were included. The semidiurnal propagating waves ex- 
cited by heating in the lower atmosphere were modelled by use of the classical tidal perturbations 
as lower boundary conditions. The TGCM was tuned by adjusting the tidal forcing term until 
calculated semidiurnal winds and temperatures matched incoherent scatter observations. The tidal 
TGCM results are consistent with previous theoretical work and successfully reproduce high 
altitude temperature and velocity data, but give somewhat lower magnitudes for velocities and 
temperatures near 160 km than are seen observationally. Future plans call for similar studies being 
done for solar maximum equinox conditions, and then the solstice cases. These studies would 
benefit enormously from additional data, particularly if higher order modes are to be accurately 
included. 
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AMPLITUDE 


SEMIDIURNAL HOUGH FUNCTIONS 



Figure 1. Classical Hough functions for the semidiurnal tides. Note increasing latitude 

structure for higher order modes. 


64 


G (Z, 0, t) = E L n (Z) Q n (0) e 1 <<* + 
n 

FOR SEMIDIURNAL TIDES 
2 7T 

y 2 

S = 2 

$97 = { $ 2 , 2g 7 0 2, 2 e ' k2 ' 2 297 + 

$2 4 g7 0 2 4 S' ^2, 4 297 } * S' * 47rt + 

Figure 2. Classical solutions for any 
tidal variable as a function of altitude, 
latitude, and time can be written as 
shown in the top line of the figure. 
The geopotential <i> at 97 km can be 
written as shown. 
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1. SPECIFY THE GEOPHYSICAL CONDITIONSj 

2. START WITH EQUILIBRIUM SOLUTIONS OF THE TGCM FOR 
SOLAR HEATING ALONE FOR THE APPROPRIATE CONDITIONSj 

3. GUESS THE AMPLITUDE AND PHASE FOR THE 2,2 AND 2,4 
GEOPOTENTIAL AT 97 KM, AND CALCULATE THE TOTAL TIDAL 
GEOPOTENTIAL AT 97 KMj 

4. CALCULATE U, V, T AT 97 KM FROM THE GEOPOTENTIAL AT 97 KM 
USING FORMULAE FROM CHAPMAN AND LINDZEN AND USE THESE AS 
LOWER BOUNDARY CONDITIONS; 

5. INTEGRATE THE MODEL EQUATIONS UNTIL EQUILIBRIUM IS 

achieved. This typically takes 10k min Cray cpu; 

6. FOURIER DECOMPOSE THE U, V, AND T FIELDS INTO ZONAL 
WAVENUMBERS 1 AND 2j 

7. COMPARE THE TGCM U, V, AND T AMPLITUDES AND PHASES AT 
18° AND 45° WITH OBSERVATIONS FROM ARECIBO, MILLSTONE 
HILL, AND SAINT SANTINj 

8. MODIFY THE INITIAL GUESSES FOR THE 2,2 AND 2,4 AMPLITUDES 
AND PHASES OF THE GEOPOTENTIAL AND REPEAT FROM STEP 4 
UNTIL REASONABLE AGREEMENT IS OBTAINED BETWEEN THE MODEL 
AND OBSERVATIONS. 

Figure 3. Steps in the calculation of the thermospheric winds and temperatures including the 

effect of tidal excitation. 
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Figure 5. The incoherent scatter radar observations to which the TGCM winds and 

temperatures were compared. 
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BEST FIT TO OBSERVATIONS FOR 

2, 2 AND 2, 4 GEOPOTENTIAL AT 97 KM: 

A 2 2 = 4 x 10 4 m 

A 2 4 = 3 x 10 4 m 

</> 2 2 ~ “ 0-5 hr 

0 2 , 4 = “2.5 hr 

Figure 6. The amplitudes and phases for the 
2, 2 and 2, 4 geopotential at 97 km which 
best fit the observations shown in Figure 5. 
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Figure 8. The TGCM north-south velocity at 18° and 45°. 
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Figure 9. The temperature observations at 18° and 45°. 
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Figure 10. The TGCM temperatures at 18° and 45°. 
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OBSERVATIONS 



Figure 1 1 . The east-west velocity observations at 45 °. 
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Figure 12. The TGCM east-west velocities at 18° and 45°. 
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Figure 14. The TGCM diurnal north-south velocity results at 45°. 
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TIDAL COUPLING WITH THE LOWER ATMOSPHERE 
(Invited Review) 


Jeffrey M. Forbes 

Department of Electrical, Computer, and Systems Engineering 
Boston University, Boston, MA 02215 


1. INTRODUCTION 

The purpose of this paper is to review the various ways in which propagating tidal components 
excited in the mesophere and below affect the structure of the thermosphere and ionosphere above 
100 km. Dynamo effects are not treated here, as they will be addressed separately in the paper by 
A. D. Richmond. 

We begin by examining the physical processes affecting the propagation of upward propagating 
tides, and how they are interrelated in the context of a numerical model. Propagating diurnal and 
semidiurnal tides which reach thermospheric heights are excited primarily by insolation absorption 
by tropospheric water vapor (0-15 km) and stratospheric/mesospheric ozone (40-60 km), respec- 
tively. Simulation of these oscillations requires consideration of mean zonal winds and meridional 
temperature gradients, and the damping effects of turbulent and molecular dissipation, radiative 
cooling, and ion drag. These effects must be considered on a spherical rotating atmosphere ex- 
tending from the ground to above 300 km, as they are in the model developed by Forbes (1982 a, 
b) depicted schematically in Figure 1 . 

2. WINDS AND TEMPERATURES 

Figure 2 illustrates amplitudes and phase vertical structures at 0°, 18°, 42°, and 60° latitude for 
the solar diurnal westerly wind at equinox from the Forbes model. In-situ EUV excitation above 
90 km is included in these simulations. The following features are worth noting: 

(1) Below 100 km at low latitudes the exponential amplitude growth and phase progression (X z 
® 30 km) with height are characteristic of the (1, 1) diurnal propagating tide. The (1, 1) mode at- 
tains its peak amplitudes near 1 10 km and decays rapidly above this height due to molecular 
dissipation. 

(2) Below 100 km at high latitudes the relative absence of amplitude growth and phase progres- 
sion with height is indicative of the (1, -2) trapped mode. Superposition of the (1, 1) and (1, -2) 
modes accounts for the illustrated changes in vertical structure of the diurnal tidal winds and 
temperatures. 

(3) Amplitudes and phases of u, v, and <5T are asymptotic to constant values above 200 km. 

This behavior is consistent with the dominance of diffusion in the upper thermosphere, and with 
the condition that there be no sources of heat or momentum in the upper thermosphere. 

(4) Diurnal tidal oscillations in the 90-150 km region receive about equal contributions from up- 
ward propagating and in situ excited components. 
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Item (4) is examined in detail in Figure 3 where the northerly velocity at 18° latitude is separated 
into relative contributions due to the (1, 1) propagating tide (predominant below 150 km) and that 
excited in-situ by EUV and UV solar radiation absorption (predominant above 150 km). Note the 
transition from a 30-km vertical wavelength phase progression with height below 150 km indicative 
of the (1, 1) mode to phase and amplitude constancy with height indicative of fast molecular dif- 
fusion and in situ excitation. 

Figures 4 and 5 illustrate amplitude and phase vertical structures at 0°, 18°, 42°, and 60° 
latitude for the solar and lunar semidiurnal westerly winds at equinox from the Forbes model. For 
the solar component the transition to shorter vertical wavelengths between 80 and 100 km, the ef- 
fects of dissipation on the upward propagating components between 110 and 150 km, and the 
asymptotic behavior characteristic of the upper thermosphere are clearly illustrated. The upper 
thermosphere velocity and temperature fields, which typically range between 10-50 m sec -1 and 
20-40 K with the larger values at low latitudes, originate with about equal weighting from three ex- 
citation sources: (1) in situ EUV excitation, (2) ion-drag momentum coupling with the diurnal tidal 
winds, and (3) upward propagating modes excited below 100 km. Although the lunar gravitational 
excitation consists of some (2, 4) in addition to the predominant (2, 2) forcing, the excitation of 
the higher-order (2, 4) and (2, 5) modes are due almost exclusively to mode coupling due to 
mesospheric mean winds and meridional temperature gradients. Winds (temperatures), in fact, 
reach amplitudes of order 10-15 m sec 1 (10-15 K) in the lower thermosphere, and 5-10 m sec 1 (5-10 
K) in the upper thermosphere, and may thus account for a significant portion of day-to-day 
variability reported in measurements of the solar semidiurnal tide. 

The joint presence of molecular viscosity, thermal conductivity, anisotropic ion drag, and rota- 
tion on a sphere renders the viscid tidal equations inseparable with respect to height and latitude, 
whereas in an inviscid atmosphere where the background temperature is independent of latitude 
the equations are separable, and classical tidal theory applies. In classical tidal theory the eigen- 
solutions (Hough functions) of Laplace’s tidal equation define the horizontal structures of each 
mode, and the eigenvalues (equivalent depths) fix each mode’s vertical structure. Thus, besides 
alteration of the vertical tidal structures from exponential growth (for propagating tides) or decay 
(for trapped tides) to asymptotically constant solutions in the upper thermosphere, the region 
where x, the ratio of the wave period to the dissipative time scale, approaches unity is also 
characterized by a transition from tidal solutions that are separable with respect to height and 
latitude to one in which vertical structures for a particular ‘modal extension’ into the thermosphere 
vary with height. This behavior is illustrated in Figures 7 and 8, which depict the horizontal shapes 
of the (2, 2) and (2, 4) Hough mode extensions (HME) of the semidiurnal temperature oscillation 
at various heights. Note that the node at 15° latitude for (2, 4) disappears and the (2, 2) horizontal 
shape broadens considerably at progressively greater heights in the thermosphere. 

An illustration of how upward-propagating tides affect the local-time structure of the thermo- 
sphere is illustrated in Figure 9, taken from Garrett and Forbes (1978). These authors superimposed 
in-situ diurnal simulations with theoretical semidiurnal ‘Hough mode extensions’ calibrated to 
observational data. Note the high degree of structure exhibited below 200 km. The hour of max- 
imum amplitude shifts to earlier times with increasing altitude, indicating the presence of upward 
propagation waves. A 12-hour period is particularly evident below 150 km. 

At solar minimum, the upward propagating components can be expected to exert a greater in- 
fluence on thermospheric structure, due to the relatively smaller contribution from in-situ EUV 
sources. An example of how this influence can be reflected in midlatitude exospheric temperatures 
is illustrated in Figures 10 and 11, where temperatures measured at Millstone Hill during 1974 
(SSMIN) and 1980 (SSMAX) are depicted. Note the predominance of a diurnal component at 
SSMAX whereas the SSMIN curve exhibits structure containing semidiurnal and terdiurnal 
components. 
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Another manifestation of upward propagating tides occurs via the nonlinear coupling between 
the semidiurnal thermospheric oscillation (which originates via in-situ as well as lower atmosphere 
contributions) and the diurnal variation of ion drag to produce a terdiurnal component to the 
thermospheric temperature and density, particularly at low latitudes where the coupling is strong. 
The nonlinear coupling between tidal winds and ion drag has been invoked by Mayr et al (1974) to 
explain spurious anamalous increases in the equatorial neutral temperature around midnight ob- 
served by the NATE experiment on AE-E, as illustrated in Figure 12. The semidiurnal component 
receives strong contributions from (a) diurnal winds interacting with the diurnal ion drag and (b) 
upward propagating modes excited by 0 3 insolation absorption. The semidiurnal fields interacting 
with the diurnal component of ion drag, in turn, generate a substantial terdiurnal component. The 
synthesis of all contributions can yield the signature of a midnight temperature maximum in the 
upper thermosphere as illustrated in Figure 12. 

3. ACCELERATION AND HEATING OF THE LOWER THERMOSPHERE DUE TO 

DISSIPATING TIDAL WAVES. 

In the zonal mean, the atmosphere can be accelerated and heated by the deposition of momen- 
tum and thermal energy by the so-called ‘eddy’ or ‘perturbation’ motions of the atmosphere. Since 
gravity wave and tidal amplitudes grow exponentially with height, it has often been suggested that 
these motions might contribute significantly to the mean momentum and energy budget of the 
lower thermosphere. These effects enter in the zonal mean momentum equation as a divergence of 
the eastward eddy momentum flux: 


1 d 1 d 

F u = (u'v ' cos 2 4>) p u 'w ' 

acos 2 0 d0 p dz 

and in the thermal energy equation as a divergence of the thermal eddy momentum flux: 

1 d 1 d 

F T (v ' cos 0) p w , 

acos 0 30 z p 3z z 


where 

a = radius of the earth 

0 = latitude 

z = altitude 

p = pressure 

u ' = perturbation westerly velocity 
v ' = perturbation northerly velocity 
w ' = perturbation vertical velocity 
3/ = perturbation geopotential 

Miyahara (1978) has investigated the deposition of mean momentum and heat in the lower 
thermosphere connected with the dissipating (1,1) and (2, 4) tidal modes using the above equa- 
tions. In Figure 13 Miyahara’ s calculation of F u and the resulting mean zonal wind are illustrated. 
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fhe (1, 1) mode is apparently capable of producing an easterly jet ( — 60 m sec -1 ) in the equatorial 
lower thermosphere, and a westerly flow of order 30 m sec -1 at midlatitudes. Results for the (2, 4) 
mode in Figure 14 indicate mean zonal winds of order 10-15 m sec 1 , which are smaller but not 
negligible compared with the (1, 1) mode. The zonal flow generated by the dissipating tidal modes 
is comparable to the flow generated by direct solar heating in this altitude regime as computed by 
Dickinson et al (1975). 

Recent calculations by Groves and Forbes (1984) indicate that the effects of tidal heating on 
thermospheric temperature structure appear secondary to the influence of tidal accelerations on the 
zonal mean wind field of the 90-150 km region. 

4. TIDAL EFFECTS ON COMPOSITION 

Mayr and Harris (1977) and Forbes (1978) have investigated tidal variations in thermospheric O, 
0 2 , N 2 , Ar, He and H using models that take into account the effects of tidal temperatures, 
horizontal and vertical tidal winds, photo- and ion-chemistry, exospheric transport, and thermal 
diffussion. Hydrogen tidal variations are dominated by vertical flow due to lateral transport in the 
exosphere, but wind induced diffusion is the single most important process for causing deviations 
from diffusive equilibrium (temperature-dominated) solutions of tidal variations in O, O z , N 2 , Ar, 
and He in the thermosphere. The effects of ion- and photo-chemistry on the variations of O and 
0 2 , and exospheric transport on He, are found to be of secondary or negligible importance above 
120 km. Measurements of neutral composition and temperature aboard the AE-E satellite have 
been analyzed to determine the semidiurnal and terdiurnal variations of O, N 2 , He and Ar from 
145-295 km (Hedin et al, 1980). The semidiurnal variations of O and N 2 are illustrated in Figures 
15 and 16, respectively, along with predictions from the MSIS model, the Forbes model, and the 
Mayr et al (1979) model. For O variations, the Forbes and Mayr et al models predict the overall 
phase and ampltitude structures quite well, with some overestimate of amplitude. Surprisingly, the 
empirical MSIS does not fit the measurements as well as the theoretical models. On the other 
hand, the MSIS model provides a much better fit to the semidiurnal amplitude of N 2 whereas the 
other models overestimate its amplitude. The semidiurnal phases for N 2 are adequately reproduced by 
Forbes and MSIS models, but the Mayr et al model yield phases about 2 h too late above 200 km. 

It is well known that winds exert an important influence on the height and density of the F-layer 
peak. A pronounced example is the dynamic behavior of the nighttime Arecibo ionosphere, par- 
ticularly the so-called ‘midnight collapse’ phenomenon, which receives its primary drive from the 
semidiurnal tide propagating upward from the lower atmosphere. A simulation of the Arecibo 
ionosphere from the model of Crary and Forbes (1984) is shown in Figure 17. The post-midnight 
descent (or ‘collapse’) of the F-layer peak is shown by these authors to be precipitated by a sudden 
abatement of southward winds determined by a reinforcement of the in-situ excited diurnal com- 
ponent, an upward propagating semidiurnal component, and a terdiurnal component which these 
authors conjecture to arise from ion-drag coupling between the semidiurnal wind and the diurnal 
component of ion drag. In addition, the bunching of contours in the bottomside F-region as il- 
lustrated in Figure 17 is due to the vertical shear of the semidiurnal wind field, and may produce 
electron density gradients sufficient to excite the ExB gradient drift instability. 

Acknowledgements. The preparation of this review was supported by NSF Grant ATM - 8319487 
to Boston University. 
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Figure 1 . Schematic describing the numerical model of atmospheric tides 
developed by Forbes (1982 a, b). 
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Figure 3. Northerly velocity at 18° latitude due to in 
situ excitation, the propagating (1, 1) Mode, and the 
sum of these components. 
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Figure 4. Solar semidiurnal westerly winds from the Forbes (1982 a, b) model. 
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Figure 6. Hough mode decompositions of solar and lunar semidiurnal temperatures excited in the 

lower thermosphere. 
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Figure 7. Normalized temperature oscillation amplitude and phase 
vertical structures of the (2, 2) Hough Mode Extension at 100, 

150, and 300 km for T 0 = 1000 K with normalizing factors equal to 
1.00, 7.00, and 4.87, respectively. 
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Figure 8. Normalized temperature oscillation amplitude and phase 
vertical structures of the (2, 4) Hough Mode Extension at 100, 

150, and 300 km, for T 0 = 1000 K with normalizing factors equal 
to 1.00, 2.91 and 1.88, respectively. 
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Figure 9. Wind and temperature contours from the model of 
Garrett and Forbes (1978). 
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Figure 10. Example of Millstone Hill neutral exospheric temperature data during a 
magnetically quiet day at solar minimum. 
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Figure 11. Example of Millstone Hill Neutral Exospheric Temperature Data During 
a Magnetically Quiet Day at Solar Maximum. 
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Figure 12. Comparisons between OGO-6 and 
MSIS models at 250 km and fourier analyzed in- 
Situ T g measurements from the NATE 
experiment on AE-E (Mayr, et al., 1979). 
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Figure 15. Atomic oxygen semidiurnal variations (Hedin, et al., 

1980) 



Figure 16. Molecular nitrogen semidiurnal variations (Hedin, et al., 

1980). 
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Figure 17. Log electron density contours and height of F-layer maximum density (dashed line) simulated for 
equinox conditions over Arecibo using the model of Crary and Forbes (1984). 




THE EFFECT OF BREAKING GRAVITY WAVES ON THE DYNAMICS AND 
CHEMISTRY OF THE MESOSPHERE AND LOWER THERMOSPHERE 

(Invited Review) 

Rolando R. Garcia 

National Center for Atmospheric Research* 

P.O. Box 3000 
Boulder, Colorado 80307 

The influence of breaking gravity waves on the dynamics and chemical composition of the 
60-110 km region is investigated with a two-dimensional model that includes a parameterization of 
gravity wave momentum deposition and diffusion. 

The dynamical model is that described by Garcia and Solomon (1983) and Solomon and Garcia 
(1983) and includes a complete chemical scheme for the mesosphere and lower thermosphere. The 
paramaterization of Lindzen (1981) is used to calculate the momentum deposited and the turbulent 
diffusion produced by the gravity waves. 

It is found that wave momentum deposition drives a very vigorous mean meridional circulation, 
produces a very cold summer mesopause and reverses the zonal wind jets above about 85 km (Fig. 
1). The momentum deposition and turbulent diffusion are much weaker at equinox than at solstice 
(Fig. 2). In Figure 3 the seasonal variation of the turbulent diffusion coefficient is consistent with 
the behavior of mesospheric turbulence inferred from MST radar echoes (Balsley et al., 1983; see 
Fig. 3 a, b). 

The large decreases in turbulent diffusion at the equinoxes have profound effects on the 
distribution of certain chemical trace species. In particular, atomic oxygen in the lower thermo- 
sphere increases at the equinoxes because loss through diffusion into the sink region below = 90 
km is reduced. This increase in atomic oxygen is reflected in the greater intensity of the 5577 A green 
line airglow emission. The model reproduces very well the observed seasonal variation of the green 
line airglow (Fig. 4). 

Ozone in the upper mesosphere is destroyed efficiently by HO x (H + OH + H0 2 ) catalysis. 

HO x abundance at the 80 km level depends sensitively on the concentration of water vapor, which 
in turn is determined by upward diffusion from the lower mesosphere. At equinox, the diffusive 
flux of H 2 0 is reduced, HO x abundances decrease and 0 3 increases in the upper mesosphere. The 
1.27 n m 0 2 airglow is a proxy for 0 3 abundances at 80 km. Fig. 5 shows that the model 
reproduces the major features of the seasonal variation of 1 .27 g. m airglow observed by the SME 
satellite (Thomas et al., 1984). 

In conclusion, the large degree of consistency between model results and various types of 
dynamical and chemical data supports very strongly the hypothesis that breaking gravity waves 
play a major role in determining the zonally-averaged dynamical and chemical structure of the 
60-110 km region of the atmosphere. A detailed account of this work is given in Garcia and 
Solomon (1984). 


*The National Center for Atmospheric Research is sponsored by the National Science Foundation 
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Figure 1. Zonally averaged circulation and temperature for model northern hemisphere winter solstice. 
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Figure 3. (a) Height distribution of MST radar echoes at Poker Flat (65 °N) as a function of 
season. Cross-hatched areas denote periods and altitudes for which no data are available. (From 
Balsley et al., 1983). (b) Vertical profiles of radar echo signal-to-noise ratio for summer (June 
18- August 21, 1980 average) and fall (September 16-October 23, 1980). (Balsley et al., 1983). 

(c) Seasonal variation of the turbulent diffusion coefficient at 61 ° north latitude. 
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Figure 4. Seasonal variation of the atomic oxygen green line (5577 A) emission intensities observed by 
the ISIS satellite at 35 °N adapted from (Cogger et al., 1981) compared with model calculations. 
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Figure 5. Top panel: Ratio of 0 3 observations by 
the Solar Mesosphere Explorer satellite (Thomas 
et al., 1984) for one week averages about day 270 
to day 1 (southern hemisphere spring/ summer; 
northern hemisphere fall/winter). Bottom panel: 
Same, as computed by the model. 
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THERMOSPHERIC NITRIC OXIDE AND ITS ROLE IN 
THERMOSPHERIC DYNAMICS AND COMPOSITION 

J-Cl. Gerard 

Institut d’Astrophysique, Universite de Liege - B4200 Cointe-Liege Belgium 

and 

R. G. Roble 

National Center for Atmospheric Research - P.O. Box 3000 - Boulder Colorado 80307 

The global distribution of thermospheric nitric oxide has been measured under different 
geophysical conditions by the Atmosphere Explorer (AE) and Solar Mesosphere Explorer (SME) 
satellites. These observations show a great variability in the NO concentration but the existence of 
a latitudinal gradient is clearly evidenced by statistical maps of the NO meridional distribution. A 
two-dimensional zonally averaged chemical-dynamical model has been used to investigate the im- 
portance of nitric oxide 5.3 jurn cooling and its role on thermospheric temperature, dynamics and 
major gas composition. 

For this purpose, a 2-D background atmosphere code and an odd nitrogen code have been 
coupled and run to steady state. The NO distribution obtained is in good agreement with the AE- 
D global picture calculated at solstice for solar minimum activity. The importance of each term in 
the thermodynamic equation has been studied spatially. It is found that the NO I.R. cooling term 
competes with conduction in the upper thermosphere and reaches its maximum value near 200km 
at high summer latitudes. 

The primary effect of including the NO cooling term is to increase the temperature in the upper 
winter thermosphere and decrease it in the other hemisphere, if the global average temperature is 
fixed. In general, the strength of the circulation is decreased and the thermal gradient is smoother. 
As a consequence of these changes in the wind and temperature fields, the 0 2 , N 2 and O densities 
also react to the effect of the NO cooling. The importance of this process depends on the level of 
solar activity which controls the NO distribution. 
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Figure 1. Log 10 of horizontal column densities (cm- 2 ) of NO observed with the UVNO experiment on the 

Atmosphere Explorer-D (AE-D) satellite. 
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Figure 2. Latitudinal distribution of log 10 of the NO peak concentration measured from AE-D (“spinning 

orbits”). 
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Figure 3. Comparison between NO peak densities obtained with the 2-D odd nitrogen model 
and AE-C, AE-D and high latitude rocket measurements. 
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THERMOSPHERIC DYNAMICS-NO COUPLING 



Figure 4. Block diagram of the coupling between the thermospheric structure and the nitric oxide concentration. 



MODEL EQUATIONS 
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(- p div v = 0 

dt 

THERMODYNAMIC: 

dT 

pc„ — = - div <£>u - P div u + Q 
dt 


MOMENTUM: 

K -j 

— = -2fi xu — grad P + g + v V 2 v - \v 
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HYDROSTATIC: 

9P _ P 

dz H 

PERFECT GAS LAW: 

P = p k T/m 


CHEMICAL: 

d (ri: rri: p) 
= Pi 


at 


- I, ( rij m, p ) - div <£, 


</>i = n i m i P v + 4>diff 


Figure 5. Equations of the 2-D model of the 
thermospheric structure. 


110 



NITRIC OXIDE I. R. COOLING 
L N o = h v [NO] v = j A-, o 
= h« A 1f 0 w [NO] V=0 

WITH: 

ki , 0 [ 0 ] 

k-i, o [O] + g 

k 1i0 : QUENCHING COEFFICIENT OF N0(v=1)BY 0 - 6.5 x 10 11 cm 3 ^ 1 
A 1 o = 13.3 s- 1 

Figure 6. Expression, for the heat loss by NO (v= 1 — 0) infrared transition. 
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Figure 7. Terms of the thermodynamic equation. 
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SOLAR ACTIVITY CONTROL 



Figure 8. Block diagram showing the response of various processes in the model to solar activity. 
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Figure 9. Contour plots of (a) the meridional wind (m.s- 1 ), (b) zonal wind (m.s -1 ), (c) perturbation 
temperature (K), (d) vertical wind (cm s~*) from the 2-D chemical-dynamical model. 


114 




ALTITUDE (km) 


356 


LATITUDE (degrees) 



LATITUDE (degrees) 


Figure 10. Contour plots of log of (a) 0 2 , (b) O and (c) N 2 from the 2-D dynamical-chemical 

model, all in cm- 3 . 










Figure 12. Contour plots displaying percentage variations of (a) 0 2 , (b) O and (c) N 2 when model 
calculations including NO IR cooling are compared to those without NO IR cooling. 






GLOBAL LARGE SCALE STRUCTURES IN THE F REGION 


S. H. Gross 

Polytechnic Institute of New York 
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Measurements of neutral densities, temperatures and ion density and temperature by the At- 
mosphere Explorer C and E satellites exhibit fluctuations made up of coherent waves of scale size 
from hundreds to thousands of kilometers. The fluctuations in the AE-C data are found at times 
of low magnetic activity to extend from the auroral region in one hemisphere across the equator to 
the corresponding region in the other hemisphere with little change in apparent wavelength. The 
fluctuations in the AE-E data similarly exhibit at times of low magnetic activity nearly uniform 
wave structure encircling the equatorial belt. Typical examples of waves of more than 1000 km ap- 
parent wavelength in both the AE-C and AE-E data will be presented. These observations during 
quiet times may be explained by a system of quasi-stationary waves excited by the auroral belts 
that are distributed about the earth possibly with tidal or sub-tidal periods and scale lengths deter- 
mined by the structure of the belts. According to this interpretation, TID’s are manifestations of 
transient wave systems moving past a ground observatory that result from changes in magnetic ac- 
tivity with periods and scale sizes associated with the regions of enhanced auroral activity. Some 
evidence suggesting tidal-like stationary waves from the Millstone Hill incoherent scatter radar is 
presented. 

The first eight figures that follow are for AE-C. The next four are for AE-E and the last is for 
Millstone Hill incoherent scatter radar measurements. It will be seen for the fluctuations in the 
data measured by AE-C that there is little variation in the scale size of the structures while the 
satellite is moving from one hemisphere to the other. Since the orbit is inclined 68°, there are por- 
tions of the orbit when the satellite is at high latitudes moving almost west to east. There is no 
reason why the behavior of the wave structure for these parts of the orbit should be the same as 
the parts in which the satellite is moving mostly from north to south or from south to north. Fur- 
thermore, as one might expect with fixed wavelength waves, there is no evidence of exponential at- 
tenuation with distance as might occur due to loss processes. The mode of propagation, if these 
are propagating waves, must be such as to incur low losses. 

Figure 1 is a plot of AE-C measured data on January 20, 1975. The orbit was nearly circular, 
inclined 68° and the satellite was at an altitude of about 250 km. Densities of atomic oxygen, 
nitrogen, helium and the electron density (Ne in the figure) are plotted in the upper graph. The 
lower graph contains the electron temperature (dotted line) and the ion temperature (heavy line). 
These are 15 seconds, United Abstract file data. (All the AE examples shown here are from these 
files, including AE-E). Magnetic activity was quite moderate on this day. 

Figure 2 is a plot of fluctuations obtained for the data of Figure 1. The data are normalized, 
detrended and passed through a high pass filter that passes all waves with wavelengths less than 
2600 km. 2000 km waves are evident in the figure. Three pairs of fluctuations are shown. The solid 
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line in each is for oxygen for comparison, the other in each, shown dotted, is as indicated. The 
three nearly even peaks in the oxygen fluctuations near the start from just before 29760.0 seconds 
UT to about 29760.0 seconds are for the part of the orbit where the satellite is at high latitudes 
moving west to east. From 29760.0 seconds to about mid-way between 30585.0 and 31410.0 
seconds there are five nearly even spaces between peaks in the oxygen fluctuations. Similar 
behavior is evident in the nitrogen fluctuations and the phase shifted ion temperature fluctuations. 
During this time the satellite was moving southward from just below 60 °N to -30 °S. Beyond, there 
is evidence of nearly equal spacing, but the magnitude is diminished. At 31410.0 seconds the 
satellite is again traversing mostly from west to east in the southern hemisphere. Exponential 
damping is not evident. 

Figure 3 contains nitrogen fluctuations for AE-C measurements on January 19, 1975, a day of 
moderate magnetic activity in which k p , did not exceed 3 0 . The same filter as for Figure 2 is 
employed. The satellite is at high latitudes moving mostly from west to east until just before 

62481.0 seconds UT when it starts moving southward. There are roughly evenly spaced waves of 
2000 km scale size until a time midway between 64173.0 and 65019.0 seconds UT. At this point 
the satellite has reached southern high latitudes and is moving mostly from west to east before 
reversing its path and returning to the nothern hemisphere. 

Figure 4 is for the same orbit as Figure 3, but the data were passed through a high pass 
filter that passes only waves with wavelengths less than 1300 km. Waves of the order of 900 km 
are evident. The north to south region is again from about 62481.0 seconds to midway between 

64173.0 and 65019 seconds. More uniformity appears in the region from this midpoint (roughly 
-60 °S) to earlier times, a point somewhat to the left of 63327.0 seconds (at about 30°N). 

Figure 5 is for another day of AE-C data, February 17, 1975. The densities of nitrogen (top 
curve), helium (second highest), argon (heavy line) and electron density (thin curve varying about 
the argon curve after 36232.5 seconds) are plotted. Two auroral zone heating signatures are evi- 
dent, one very distinct on the left at high northern latitudes, centered somewhat to the left of 
34856.3 seconds, and a second, less distinct, at high southern latitudes centered to the left of 
37608.8 seconds. The heating in both regions may act as wave sources. The k p for this day did not 
exceed 3 + . 

Figure 6 is a plot for the fluctuations of nitrogen (solid line) and argon (dotted line) for the data 
in Figure 5. It was obtained by passing the data through a bandpass filter allowing only waves in 
the range 1100 to 2600 km to pass. Strong waves to the left of 34971.0 seconds are for the nor- 
thern auroral region. Similar strong waves are observed in the southern auroral region to the left 
of 37509.0 seconds. The equator is crossed at about 35817.0 seconds. The pattern looks like that 
of two interfering wave systems. High latitude west to east traversal in the northern hemisphere is 
for the time interval 34125.0 to 34971.0 seconds, and somewhat to the right of 36663.0 to 
somewhat past 37509.8 seconds for the southern hemisphere. 

Figure 7 contains the cross correlation function for the fluctuations of Figure 6. The large 
negative peak of -.8 shows that argon and nitrogen are well correlated and that a coherent wave 
system exists. 

Figure 8 contains fluctuations for the data of Figure 5 but for a bandpass filter passing waves in 
the range 800 to 1300 km. These are roughly 1000 km waves of nitrogen fluctuations. The 
southern region at about 36663.0 seconds is stronger by a factor of about 2. The fluctuations 
resemble that of two interfering waves, perhaps from two sources. Fairly evenly spaced waves are 
evident between 34971.0 and 36663.0 seconds. 

Figure 9 is a plot similar to that of Figure 1, but is for AE-E measurements on August 1, 1977. 
AE-E was in a nearly circular orbit inclined 20° at an altitude of about 260-270 km. The upper 
graph contains densities of oxygen (upper curve), nitrogen (second highest), helium (third curve) 
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and the electron density (lowest). The lower graph contains ion temperature (dotted line) and 
neutral gas temperature (heavy line). Magnetic activity was moderate; k p did not exceed 3 . The 
satellite was mainly moving from west to east in an equatorial belt between ±20°. 

Figure 10 contains plots of four pairs of fluctuations for the data of Figure 9 as obtained from 
a high-pass filter passing only waves with wavelengths less than 2600 km. The scale size of the 
waves is about 1800 km. The solid line in all four graphs is for oxygen, the others (dotted) as in- 
dicated. The time shown is for about 3/4 of an entire orbit, limited due to the absence of data 
prior to and after the orbit. Nearly uniform waves are evident girdling this much of the equatorial 
belt. 

Figure 11 is for AE-E for November 23, 1977. The satellite was in the same orbit as in Figure 9. 
Neutral and ion winds are shown as well. Magnetic activity was low for this day, k p did not ex- 
ceed 1 0 . 

Figure 12 is a plot of fluctuations for the data of Figure 11. Here the data were passed through a 
bandpass filter allowing only waves in the range 2600 to 5200 km to pass. The waves in the figure 
are roughly of 2800 km in scale size. The upper graph is for neutral temperature and nitrogen den- 
sity fluctuations. Its cross correlation function is to its right and shows a negative peak of about 
-.85. Fluctuations of horizontal (roughly north-south) neutral winds and nitrogen density are 
shown in the lower graph. Its cross correlation function is to its right and contains a negative peak 
of about -.68. The time period covers about one half of an orbit, the limitation being the extent of 
the measurements. 

Figure 13 is for the data of Figure 11, but the fluctuations are obtained from a high-pass filter 
that passes all wavelengths less than 1300 km. These waves are about 1000 km in scale size and the 
plot covers about half the equatorial belt. The uniformity of these waves is quite remarkable. The 
solid line is for nitrogen and the dotted for the neutral temperature. 

Figure 14 is a plot of data obtained from a series of Millstone Hill incoherent scatter radar 
measurements while in an elevation scan mode from due south to due north. The scan is com- 
pleted in about 20 minutes, but the scan to scan period is 45 minutes. Shown here are the so-called 
raw fluctuations for individual scans plotted with separations so that they can be seen relative to 
each other. Not all scans are shown. The plots are to a scale which is not shown. The abscissa is 
the geomagnetic latitude. The Millstone Hill station is at 54°. The curves shown are constant 
altitude electron densities at 350 km. The day and beginning time of each scan is shown on the left 
of each curve. The nature of the variation appears to change from flat to concave to flat to con- 
vex, etc., with a period that may relate to the diurnal or semi-diurnal period and for a half 
wavelength of about 1000 km, much like the behavior of a stationary wave. There is an unfor- 
tunate four hours 45 minutes gap between th 03:11 and 07:55 curves due to a radar malfunction. 
The change in shape for these two times is quite interesting. Whether these variations can be inter- 
preted as evidence of tidal like, stationary variations is highly speculative at this time. 

Measurements on other days are also somewhat like these curves. 
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Figure 1. 
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THE MEASURED MOTIONS INSIDE EQUATORIAL PLASMA BUBBLES 

W. B. Hanson 
Center for Space Sciences 
P.O. Box 688 

The University of Texas at Dallas 
Richardson, TX 75080 

D. K. Bamgboye 
Physics Department 
University of Ilorin 
Ilorin, Nigeria 

A preliminary study of the vertical and north-south horizontal ion motions in plasma bubbles in 
the near-equatorial ionosphere utilizing drift meter data from Atmosphere Explorer E is presented. 
High-resolution data show that the vertical ion velocity in some bubbles increases approximately 
linearly with (N 0 -N)/N, where N 0 is the background ion concentration and N is the bubble ion 
concentration. At sufficiently large N 0 /N the vertical ion velocity saturates, but often at a value 
substantially larger than the ratio of the gravitational acceleration to the ion neutral collision fre- 
quency (g/v in ), which is the nominal collision-dominated velocity limit for cylindrical bubbles. 
These larger than nominal velocities may result from background eastward electric fields and/or 
from a vertically elongated bubble cross section. The unanticipated observation that large 
poleward horizontal drifts accompany these vertical drifts seems to follow naturally from a 
redistribution of plasma along flux tubes as the plasma convects from the bottomside of the F 
region to high altitudes. 
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Figure 1. Ion number density and drift data (absolute points only) from the ion drift meter showing (in the lower 
trace) the density depletions (bubbles) and their horizontal (middle trace) and vertical (solid dots) velocities. The 
horizontal velocities are in a direction normal to the satellite velocity vector. The data are from AE-E orbit 20496 
on July 21, 1979, when Ap = 7, Kp = 2, and F 10 7 = 144. 
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Figure 2. High-resolution version of the portion of Figure 1 marked HR. Note how depletions and 
velocities that are barely visible in Figure 1 can be quite clearly identified and measured in this 
format. Alternate 4-s segments of pitch (vertical) and yaw (horizontally) velocity data are shown in 
the upper trace. 
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(b) 



Figure 3. Vertical ion velocity data from the bubble labeled A in 
Figure 2 (a) plotted versus (N 0 -N)/N 0 and (b) plotted versus 
(N 0 -N)/N. 
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Figure 4. A plot of V v versus (N 0 -N)/N for one of the very low concentration bubbles of Figure 1. Note the 

broken scale for AN/N. 
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Figure 5. Plots versus altitude of the ion-neutral collision frequency (v in ), the ratio of gravitational 
acceleration to v in , and the scale, L, that divides the collisional and collisionless regimes of the 
Rayleigh-Taylor gravitational instability. v in and g/v in are plotted for three different levels of solar 
activity, the MSIS model (Hedin, 1983) and collision cross sections from Banks and Kockarts 
(1973) were used to evaluate v in . The heavy dashed line (Lj) defines v 2 in = 4g/L for F(10.7) = 
175, and the light dashed line (L 2 ) does the same for F(10.7) = 75. 








Figure 6. Vertical velocity variations with dip latitude for 
maximum depletions at fixed heights of 270, 405, 425, and 450 
km. The number by each data point refers to the number of 
bubbles having that point as their mean maximum depletion 
velocity (low-resolution data). 
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Figure 7 . (a) The scatter plot and (b) the mean variation of the vertical velocity at maximum 
depletion as functions of dip latitude. In Figure 7b the standard error of each point is shown by a 
vertical bar. The data refer to bubbles with depletions of two orders of magnitude or less. 
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Figure 8. Variation of vertical velocity (at maximum depletion) with altitude 
at dip latitudes 7 °, 8°, 10.5°, 12 °, 13°, and 14°. 
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Figure 9. Plots of the horizontal ion velocity, V H , versus 
depletion amplitude for the bubble labeled B in Figure 2. Figure 
9a uses AN/N 0 and Figure 9b uses AN/N for the abscissa. 
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Figure 10. Plot of the maximum depletion horizontal velocity observed from the high-resolution 
data of AE-E orbit 20496 of July 21, 1979, as functions of dip latitude. 
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Figure 11. (a) Values of the north-south horizontal velocity at 
maximum depletion and (b) their mean variation, plotted as 
functions of dip latitude for bubbles with depletions of two 
orders of magnitude or less. In Figure lib the standard error 
of each mean point is shown as a vertical bar. 
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Figure 13. (a) The mean resultant ion velocity at maximum depletion as a 
function of dip latitude for bubbles from all the data with two-decade 
depletion or less, (b) The variation of the angle between the resultant 
velocity and the horizontal. The dashed lines give 4 > if there were no field- 
aligned motion. 


150 


(uj*) M V 



Figure 14. A plot of the altitude difference, Ah, between a satellite at radial distance r G (6700 km) 
and a given dip latitude to the equatorial crossing point of its dipole field line, Ah = r c tan 2 X = 
r eq sin 2 X. 
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UT/LONGITUDINAL VARIATIONS OF COMPOSITION IN DE DATA 


A. E. Hedin 
Code 614 

NASA/Goddard Spaceflight Center 
Greenbelt, MD 20771 

Composition data obtained by the quadrupole mass spectrometer (NACS) carried aboard the DE- 
B satellite have been examined for variations in the polar regions during magnetically quiet condi- 
tions. The corresponding predictions of the MSIS-83 model are often used for comparison. The 
MSIS-83 model is used to emphasize the variations of interest by suppressing all other variations. 

A persistent enhancement of N 2 density and a depletion of He density are present in the vicinity of 
the magnetic poles with maximum density response in the morning hours (magnetic) on average. A 
universal time (UT) variation in average density levels is evident near both the geographic and 
magnetic poles. There are systematic morphology changes with UT and between summer and 
winter which are qualitatively consistent with the simple concept that thermospheric heating effects 
are shifted or spread in the downwind direction of the global circulation systems driven by EUV 
and magnetospheric sources. The magnitude of the UT variations is larger in the southern 
hemisphere and is larger in local winter than summer. 
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Figure 1 . Contour plot in geomagnetic latitude (-90 to -30 degrees) and magnetic local time 
coordinates of southern hemisphere N 2 data and corresponding MSIS model results for the same 
data distribution: (a) logarithm of N 2 data divided by the MSIS model omitting the UT/longitude 
terms of the model; (b) same as panel (a) using full MSIS model N 2 densities in place of data; (c) 
and (d) same as panels (a) and (b) but dividing by MSIS model omitting the UT/longitude, local 
time, and time independent latitude terms. 
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Figure 2. Contour plot in geomagnetic latitude (-90 to -30 degrees) and magnetic local time 
coordinates of southern hemisphere He data and corresponding MSIS model results for the same 
data distribution: (a) logarithm of He data divided by the MSIS model omitting the UT/longitude 
terms of the model; (b) same as panel (a) using full MSIS model He densities in place of data; (c) 
and (d) same as panels (a) and (b) but dividing by MSIS model omitting the UT/longitude, local 
time, and time independent latitude terms. 
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Figure 3. Contour plot in geographic latitude (-90 to -30 degrees) and longitude coordinates of 
southern hemisphere He data and corresponding MSIS model results for the same data 
distribution: (a) logarithm of He data divided by the MSIS model omitting the longitude terms of 
the model; (b) same as panel (a) using full MSIS model He densities in place of data; (c) and (d) 
same as panels (a) and (b) but dividing by MSIS model omitting the longitude and time 
independent latitude terms. The x indicates the geomagnetic pole and square the dip pole. 
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Figure 4. Contour plot in geomagnetic latitude (-90 to -30 degrees) and 
magnetic local time coordinates of southern hemisphere He data for six 
hour UT intervals and corresponding MSIS model results for the same 
data distribution: (a) logarithm of He data divided by the MSIS model 
omitting the UT/longitude terms of the model for 0 to 6 hrs UT; (b) 
same as panel (a) using full MSIS model He densities in place of data; 
(c) and (d) same as panels (a) and (b) for 6 to 12 hours UT; (e) and (f) 
for 12 to 18 hrs UT; (g) and (h) for 18 to 24 hrs UT. 
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Figure 5. Contour plot in geographic latitude (-90 to -30 degrees) and 
longitude coordinates of southern hemisphere He data for six hour UT 
intervals and corresponding MSIS model results for the same data 
distribution: (a) logarithm of He data divided by the MSIS model 
omitting the UT/longitude terms of the model for 0 to 6 hours UT; (b) 
same as panel (a) using full MSIS model He densities in place of data; 
(c) and (d) same as panels (a) and (b) for 6 to 12 hours UT; (e) and (f) 
for 12 to 18 hrs UT; (g) and (h) for 18 to 24 hrs UT. 
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Figure 6. Contour plot in geographic latitude (-90 to -30 degrees) and longitude coordinates of 
southern hemisphere He data for six hour UT intervals. Contours indicate the logarithm of He 
data divided by the MSIS model omitting the UT/longitude and time independent latitude terms 
of the model for (a) 0 to 6 hrs UT; (b) 6 to 12 hrs UT; (c) 12 to 18 hrs UT; and (d) for 18 to 24 
hrs UT. 
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Figure 7. Contour plot in geographic latitude (-90 to -30 degrees) and longitude coordinates of 
southern hemisphere He data: (a) logarithm of He data divided by the MSIS model omitting the 
longitude terms of the model for summer data; (b) same as panel (a) but dividing by MSIS model 
omitting the longitude and time independent latitude terms; (c) and (d) same as (a) and (b) but for 
winter data. 
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Figure 8. Contour plot in geographic latitude (-90 to -30 degrees) and longitude coordinates of 
northern hemisphere He data: (a) logarithm of He data divided by the MSIS model omitting the 
longitude terms of the model for summer data; (b) same as panel (a) but dividing by MSIS model 
omitting the longitude and time independent latitude terms; (c) and (d) same as (a) and (b) but for 
winter data. 


161 



PEAK AMPLITUDE PEAK AMPLITUDE 




Wi Eq Su Wi Eq Su Wi Eq Su 


Figure 9. Magnitude of the polar peak, angular distance from the geographic pole, and amplitude 
of the UT variation at the geographic pole as a function of season for He and N 2 data. N 
indicates northern hemisphere and S indicates southern hemisphere. 
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CONVECTION PATTERN MORPHOLOGY AND VARIATIONS 

(Invited Review) 

R. A. Heelis 
Center for Space Science 
P.O. Box 688 

The University of Texas at Dallas 
Richardson, TX 75080 

A fairly straightforward consideration of the interaction between a southward interplanetary 
magnetic field and the earth’s magnetic field will result in the prediction of a two-cell convection pat- 
tern in the high latitude ionosphere. This is shown schematically in Figure 1 where the antisunward 
flow at high latitudes results from the application of the solar wind electric field to the ionoshpere 
and the return sunward flow results from an electric field generated in the plasma sheet to ensure 
continuity. 

This convection pattern can be quite easily characterized in terms of the radius of the approx- 
imately circular region containing the anti-sunward flow, called the polar cap and the maximum 
potential difference applied across this region. The polar cap potential difference can be shown to 
be a sensitive function of the orientation of the interplanetary magnetic field and consideration of 
magnetic flux conservation can predict a dependence of the polar cap radius on the polar cap 
potential (Siscoe, 1982). The size of the polar cap is predicted to vary as the 3/16 power of the 
potential difference. This dependence is not strong enough to be unequivocably demonstrated by 
the data. However, Figure 2 shows the dependence of the polar cap potential drop on the ratio of 
electric fields inside and outside the merging region (Reiff et al., 1981). Merging theory simply ex- 
presses this ratio as the cosine of the angle between the IMF and a plane perpendicular to the sun- 
earth line (i.e. Bz/Bx) (Hill, 1974). 

It can be seen from Figure 2 that the polar cap potential increases as the IMF becomes more 
southward and that a relatively large potential drop remains even when the IMF turns northward. 
This residual potential drop has been attributed to a ‘viscous interaction’ process at the flanks of 
the magnetosphere that exists in addition to the dayside merging process that operates when the 
IMF has a southward component. A more detailed analysis of the behavior of the potential dif- 
ference when the IMF turns northward is shown in Figure 3 (Wygant et al., 1983). It displays an 
almost exponential decay, requiring some 2-4 hours to reach a residual value near 15 Kv. 

In addition to a variation in the polar cap potential drop as a function of the IMF orientation 
examination of the signatures of ionospheric convection from satellites show that the cell 
geometry or cell shape is also dependent on the IMF. First observations of the flow geometry near 
local noon showed that substantial east- west flows could occur in this region. An arbitrary inspec- 
tion of the data also showed that this flow could be predominantly to the east or predominantly to 
the west. Examples of these differing flow signatures are shown for almost identical orbital condi- 
tions in Figures 4 and 5. A compilation of such data without strict consideration of the IMF 
orientation leads to the concept of a restricted region of sunward flow near local noon sometimes 
called the ‘throat’. Figure 6 shows a schematic of such a flow geometry generated from a 
mathematical model. Studies of the electric field near the dawn and dusk convection reversals had 
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already revealed a dependence of electric field magnitude on the IMF (Heppner, 1972). It was 
known that the antisunward flow on the poleward side of the convection reversal was larger on 
the dusk (dawn) side than on dawn (dusk) side if the IMF By component was negative (positive). 

A rotation of the throat region was postulated to account for the previously known dependence 
of the dawn/dusk flow speeds on the IMF but it was recognized that the data base of convection 
signatures should be examined for different IMF orientations and that even then the single cuts 
available from satellites like AE and DE could not define the entire high latitude convection 
geometry. Examination of the AE data base shows that reproducible signatures of the convection 
pattern can be obtained if we restrict our attention to times when the IMF has a southward com- 
ponent (Bz < 0) and a normal garden hose orientation (BxBy <0). Under these conditions the 
pattern can be generally characterized by dividing it into signatures seen when By is positive and 
when By is negative. 

Figure 7 shows examples of the signature seen in the northern high latitude ionosphere when By 
is positive (and Bx is negative). The most reproducible feature is that the flow near local noon is 
directed toward dawn at all invariant latitudes between 70 and 80 degrees. If we trace the position 
of the zero potential line that separates cells in which the flow is either clockwise around the dawn 
side or anti-clockwise around the dusk side, we note that it shows a tendency to be displaced 
slightly to the dawn side of local noon. This results in the dusk convection cell being relatively 
large and almost circular in shape while the dawn cell is apparently more cresent shaped and abuts 
the dusk cell. A more detailed description of the convection geometry for this orientation of the 
IMF has been obtained by using simultaneous data from the DE satellite and ground based radars 
at Chatanika and Millstone Hill. This data set provides multiple data points around the high 
latitude region at which the electrostatic potential and flow direction are specified. The results of 
an examination of this data set are shown in Figure 8 and confirm our expectations. 

Figure 9 shows examples of the northern hemisphere high latitude convection signature seen 
when the IMF By is negative (and Bx is positive). Contrasting these signatures with those shown in 
the previous two figures it can be seen that the flow with an anti-sunward component near local 
noon is directed toward dusk rather than toward dawn. In extreme cases when the magnitude of 
By greatly exceeds Bx the flow may be directed eastward as shown in the last two examples. In all 
cases however the flow signature is consistent with a more circular perimeter to the dawn cell with 
a cresent shaped dusk cell that surrounds it. To my knowledge a study of multipoint 
measurements to confirm such a flow geometry has not yet been completed. Figure 10 shows 
schematically the dayside high latitude flow geometry that exists in the northern hemisphere for 
different configurations of a southward IMF. 

When the IMF has a northward component (Bz > 0) the observed convection velocities are 
usually much more structured both in the summer and in the winter hemispheres. Nevertheless 
large scale convection signatures can often be recognized and reproducible features emerge with a 
dependence on the IMF orientation. Figure 1 1 shows the most easily identified feature from S3-2 
electric field data (Burke et al., 1982). Here sunward convection exists at the very highest magnetic 
latitudes where anti-sunward convection would normally be expected if the IMF had a southward 
component. Figure 12 shows the ion convection signature observed by DE-B during times of 
northward IMF. Similar regions of sunward convection are seen in this data base but additional 
features dependent on By can also be identified. Regions of anti-sunward convection on one or 
both sides of this high latitude sunward convection allow the identification of one and sometimes 
two convection cells within a region that can be associated with open field lines. These convection 
cells are labelled I and II in the figure. In addition to the convection cells within the polar cap 
there usually exists additional convection cells at lower latitudes that circulate in the same sense as 
the cells existing during times of southward IMF. These cells are labelled III and IV in the figure. 
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Examination of a large data base shows a dependence of the number of convection cells and 
their sense of circulation on the By component of the IMF. Figure 13 shows schematically these 
dependences for the high latitude northern hemisphere. When By is close to zero two convection 
cells apparently exist within the polar cap. They occupy the dawn and dusk sides and circulate 
anti-clockwise and clockwise respectively to produce sunward flow at the highest latitudes. At 
lower latitudes cresent shaped cells form about each of the polar cap cells and circulate in the manner 
expected for a southward IMF. When By is negative the dawn-side polar cap convection cell tends 
to disappear giving way to a slightly larger dusk-side cell and more turbulent flow on its dawn . 
side. The two lower latitude convection cells remain. When By is positive the polar cap dawn-side 
cell remains and perhaps expands and more turbulent flow exists on its dusk side. Again the lower 
latitude convection cells remain. The completeness of the convection cells is extremely difficult to 
observe since the nightside high latitude convection is extremely turbulent during times of north- 
ward IMF. The convection pattern I have described is at least understandable in terms of solar- 
wind/magnetosphere interaction in which open field lines are recirculated in the polar cap and a 
viscous interaction process exists near the flanks of the magnetosphere’s equatorial plane giving rise 
to the lower latitude convection cells. 
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Figure 3. 
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INTERACTION OF THE EQUATORIAL MIDNIGHT PRESSURE BULGE AND 
THERMOSPHERIC ZONAL WINDS 


F. A. Herrero 
Code 614 

NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 


The zonal component of the thermospheric neutral wind in the equatorial region is found to 
flow eastward all night, having a minimum near midnight. The data, obtained by Spencer et al. 
(1981) on the Dynamics Explorer-2 (DE-2) satellite, and reported more recently by Wharton et al. 
(1984), indicate that the minimum feature is significant. The minimum is illustrated in Figure 1 
which shows a 4th order Fourier series fit to the data points reported by Wharton et al. (1984). 

The table below the plot gives the amplitudes and phases of the four harmonic components used 
with the estimated errors for each amplitude. One asks whether this behavior is due to the local 
passage of the equatorial midnight pressure bulge associated with the midnight temperature 
anomaly. If so, it is to be exptected that other data characterizing the midnight temperature 
anomaly should show consistency with this observation, especially in view of the large changes 
observed in the zonal velocity in the midnight sector. Consistency in momentum conservation may 
be checked using available data from previous independent experiments. Such data, taken under 
similar conditions of solar activity, has been substituted into the momentum equation and found 
to be consistent using a priori estimates of the effect of viscosity. 

The DE-2 data provide the zonal wind Ux, and its local time derivative 3U x /3t. The pressure 
gradient is obtained from the average equatorial nighttime neutral temperatures measured on AE- 
E (Herrero and Spencer, 1982) and the MSIS neutral density (Hedin, 1983), and this is shown in 
Figure 2. Specifying the ion-drag requires in addition the eastward ion-drift V ix and the ion density 
n ; at the altitude of interest, 350 km in this case. The eastward ion-drift is well known near the 
equator from the Jicamarca measurements of 1970-71 (Woodman, 1972; Fejer et al., 1981; see 
Figure 3). Solar activity for that period was comparable to solar activity during the DE-2 
measurements of 1981-82, essentially one solar cycle apart. The average ion density variation may 
be represented by Chiu’s empirical model (Chiu, 1975), and this is plotted in Figure 4. Basically, 
zonal momentum is balanced as the pressure gradient dp/dx of the neutral gas is opposed by ion- 
drag and viscosity with the difference appearing as a local rate of change dUx/dt. That is, 


aux 

at 


1 9p + ^ a 2 u x 

P 3x p 3z 2 


Vin (U x - V ix ) 


where the ion-drag term is characterized by the ion frequency v in and the difference in velocities 
U x - V ix . ii is the coefficient of viscosity and p the neutral mass density (see, for example, 
Rishbeth and Garriott, 1971). The similarity between the nighttime zonal winds and the eastward 
ion-drift is worth noting as it is responsible for a drastic reduction in the ion-drag term. U x and 
V ix both reach their highest maxima togehter at 2100 hours LT, and pass through a minimum 
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shortly after midnight and a secondary maximum between 0300 and 0400 hrs LT. This similarity 
makes it possible to approximate V ix in terms of U x . From Figures 1 and 2, V ix ~ 0.8U X . Thus, 
the ion-drag term may be written as U x /r ( where r { = 1/0.2 v in , where v in is proportional to n;. 
The viscosity term is the only component of the equation that cannot be characterized by actual 
data. Therefore, it can be determined in terms of the others, and its magnitude checked for con- 
sistency against estimated values. The table below gives the magnitudes of the terms entering the 
momentum equation for the four time regions shown in Figure 2. The values of U x , dU x /dt, Ni 
and -1 / p' (dp/dx) used correspond to the times near the middle of each region. The viscosity is 
found from the net effect of the terms in the Table and compared to the a priori estimates. In 
regions I, II, and III the net effect gives a viscosity value varying between -.005 and + .005 m/s 2 . 
This amounts to 30% or less of the dominant term in each case, and perhaps is not indicative of 
the average behavior. However, in the early morning hours (section IV) the viscosity term should 
account for .01 m/s 2 , and the fact that it is positive here may be significant and consistent with 
the simultaneous reversal in the pressure gradient shown in Figure 2. 

A priori estimates of the viscosity term follow from approximating this term using a “viscosity” 
scale height H v . This gives a viscosity decay time t v ~ p H 2 /p. Previous H v estimates (Rishbeth, 
1972) indicate values of the order of 100 km which are consistent with the numbers obtained here. 


Table 1 



I 

II 

III 

IV 


(1930 LT) 

(2230 LT) 

(0130 LT) 

(0430 LT) 

Tj (hrs) 

0.81 

2.1 

6.2 

7.2 

U x /Ti (m/s 2 ) 

0.035 

0.014 

0.004 

0.002 

dU x /dt (m/s 2 ) 

0.020 

-0.015 

0.010 

-0.025 

-1/p (dp/dx) (m/s 2 ) 

0.050 

0.004 

0.015 

-0.012 
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Average altitude 350 km 
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AMPLITUDES AND PHASES 


ORDER 

AMP(m/s) 
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-1.5 

2 
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-10.8 


Figure 1 
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ION NEUTRAL COUPLING IN THE HIGH LATITUDE THERMOSPHERE: PART I 

T. L. Killeen 

Space Physics Research Laboratory 
The University of Michigan 
Ann Arbor, MI 48109 

Measurements of the neutral wind in the polar F-region from Dynamics Explorer-2 (DE-2) have 
been used to illustrate asymmetries in the neutral circulation that are dependent on the sign of the 
B y component of the interplanetary magnetic field (IMF). Individual DE-2 orbits and averaged 
data sets from different Universal times are presented. The data are categorized according to the 
sign of the hourly averaged IMF B y component measured by ISEE-3 for the hour preceding the 
DE-2 measurement. The major features observed are: 1) an asymmetry in the polar cap neutral 
flow velocity with the region of most rapid anti-sunward flow shifting from the dawn-side to the 
dusk-side of the polar cap as B y changes from positive to negative; 2) a shift in magnetic local 
time of the region of entry of neutral gas into the polar cap from a location on the dawn-side of 
the noon-midnight meridian for B y positive to one more biased towards the dusk-side for B y 
negative; 3) an enhancement in the velocities associated with the dawn, anti-clockwise neutral 
vortex for B y negative relative to those observed for B y positive. The B y neutral wind asymmetries 
can be explained by similar asymmetries, previously observed, in the polar ion convection pattern. 
They imply a direct causal relationship between solar wind/magnetosphere coupling and neutral 
thermospheric dynamics. 
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Figure 1. Average neutral wind vectors for December 1981 plotted in geomagnetic polar 
coordinates (magnetic latitude and local time) obtained for the specified UT intervals with a) 
B y negative and b) B y positive. 
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Figure 2. Vector neutral wind measurements for individual orbits of DE-2; a) orbit 1861, b) orbit 1755, c) orbit 
1853 and d) orbit 1749. The winds are plotted in geomagnetic polar coordinates. The solar terminator is indicated 
by the curved hatched line. Where no FPI data are available, the WATS measurements are indicated by the bars 
plotted at right angles to the track of the satellite. The hourly averaged IMF values from ISEE-3 taken for the 
hour preceding the pass are shown at left with the Kp and Ap indices. The wind scale is given at bottom right. 




ION-NEUTRAL COUPLING IN THE HIGH LATITUDE THERMOSPHERE: PART II 


T. L. Killeen 

Space Physics Research Laboratory 
The University of Michigan 
Ann Arbor, MI 48109 


On 24th November, 1982, The North-South (B z ) component of the Interplanetary Magnetic 
Field (IMF) became positive for a period of about 1 1 hours reaching a relatively large and steady 
value of —25 nT. During this rare occurrence, the Dynamics Explorer-2 (DE-2) spacecraft was in 
a configuration that enabled the dynamics of both ionic and neutral species of the high-latitude 
F-region to be mesaured simultaneously along the track of the polar- orbiting satellite. Results from 
two Northern (winter) polar passes of DE-2, extracted from a larger data set, are shown to il- 
lustrate the response of the neutral F-region to ion drag forcing arising from a configuration of 
ion convection characteristic of strongly northward IMF. The measured neutral winds differ ap- 
preciably from those more commonly observed for periods of southward IMF. The multi-cellular 
ion drift pattern associated with positive B z is observed to drive a similar but less structured and 
weaker neutral wind configuration in the winter polar cap. Major features of the ion drift pattern 
are mimicked by the neutral circulation but smaller-scale and more irregular structures of ion flow 
are not. This is ascribed to the relatively long time constant (few hours) for momentum exchange 
between the ion and neutral gases. The results demonstrate that sunward flow of neutral gas can 
be established and maintained by ion drag in the central polar cap for positive B z . 
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Figure 1 . Variation of the Interplanetary Magnetic Field as measured by ISEE-3 during 
24th November, 1982. The times corresponding to the orbits studied are indicated. 
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Figure 2. Neutral wind and ion drift vectors for orbit 7210 plotted in geographic polar coordinates (latitude, pole 
to 40 °N; local solar time). The neutral winds are coded by the yellow (white) arrows, and the ion drifts are 
coded by the red (gray) bars. The curved line represents the solar terminator (90% solar zenith angle). The 
symbol N refers to the location of the invariant pole at the UT of the pass. Note the scale difference of a factor 
of 2 indicated at lower right. 
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Figure 3. Neutral wind and ion drift vectors for orbit 7212 plotted in geographic polar coordinates (latitude, pole 
to 40 °N; local solar time). The neutral winds are coded by the yellow (white) arrows, and the ion drifts are 
coded by the red (gray) bars. The curved line represents the solar terminator (90% solar zenith angle). The 
symbol N refers to the location of the invariant pole at the UT of the pass. Note the scale difference of a factor 
of 2 indicated at lower right. 
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Figure 4. Zonal neutral winds and ion drifts measured on DE-2 for a) orbit 7210 and 
b) orbit 7212 plotted as a function of latitude along the track of the satellite. Neutral 
zonal winds are plotted according to the scale shown on the left, ion zonal drifts 
according to the scale on the right. 
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GEOSTROPHIC ADJUSTMENT IN A SHALLOW- WATER NUMERICAL MODEL 
AS IT RELATES TO THERMOSPHERIC DYNAMICS 

M. F. Larsen 1 and I. S. Mikkelsen 2 

'Dept, of Physics and Astronomy 
Clemson University 
Clemson, SC 29631 

2 Geophysical Division 
Meteorological Institute 
Copenhagen, Denmark 

The theory of geostrophic adjustment and its application to the dynamics of the high-latitude 
thermosphere have been discussed by ourselves (Larsen and Mikkelsen, JGR , 1983; Mikkelsen and 
Larsen, JGR , 1983) and Walterscheid and Boucher ( JAS , 1984) in previous papers based on a 
linearized treatment of the fluid dynamical equations. However, a linearized treatment is only 
valid for small Rossby numbers given by Ro = V/fL, where V is the wind speed, f is the local 
value of the Coriolis parameter, and L is a characteristic horizontal scale for the flow. For typical 
values in the auroral zone, the approximation is not reasonable for wind speeds greater than 25 
m/s or so. We have developed a shallow-water (one layer) model that includes the spherical 
geometry and full nonlinear dynamics in the momentum equations in order to isolate the effects of 
the nonlinearities on the adjustment process. A belt of accelerated winds between 60° and 70° 
latitude was used as the initial condition. We found that the adjustment process proceeds as ex- 
pected from the linear formulation, but that an asymmetry between the response for an eastward 
and westward flow results from the nonlinear curvature (centrifugal) terms. In general, the 
amplitude of an eastward flowing wind will be less after adjustment than a westward wind. For in- 
stance, if the initial wind velocity is 300 m/s, the linearized theory predicts a final wind speed of 
240 m/s, regardless of the flow direction. However, the nonlinear curvature terms modify the 
response and produce a final wind speed of only 200 m/s for an initial eastward wind and a final 
wind speed of almost 300 m/s for an initial westward flow direction. Also, less gravity wave 
energy is produced by the adjustment of the westward flow than by the adjustment of the 
eastward flow. The implications are that the response of the thermosphere should be significantly 
different on the dawn and dusk sides of the auroral oval. Larger flow velocities would be expected 
on the dusk side since the plasma will accelerate the flow in a westward direction in that sector. 
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Figure 1. The zonal wind (left), meridional wind (middle), and height of the fluid surface (right) at 
1 hr intervals. The initial condition at 0 hr is a 1 m/s eastward wind in the latitude band between 
60° to 70°. The height of the fluid is directly related to the pressure. All subsequent figures will be 
presented in the same format. 
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MIHALJAN'S SOLUTION 



Figure 2. Analytic solution for the adjusted flow based on a linearized set of equations in a 
Cartesian coordinate system. The initial wind was 100 m/s, and the final wind is 80 m/s. The 
initial and final flow velocity will scale linearly together. Thus, an initial wind of 1 m/s will 
produce an adjusted wind of 0.8 m/s. This figure can be compared directly with the solutions 
shown in Figure 1. 
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Figure 3. An enlarged graph showing the adjusted state after 10 hr. The initial wind is eastward at 
1 m/s between 60° and 70 °N. The final zonal velocity agrees closely with the zonal velocity given 
by the linearized analytic solution. The meridional wind is negligible. The height of the free 
surface, which is directly related to the pressure, shows an asymmetry across the channel. The 
spherical geometry causes a focusing at the pole. As a result, the heights are higher on the 
southern side of the belt and lower on the northern side. 
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Figure 4. Adjusted state in the shallow-water model after 10 hr for an initial eastward wind of 100 
m/s. The final zonal wind is approximately 75 m/s which is less than the final value of 80 m/s 
based on the solution of the linearized set of equations. 
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GEOSTROPHIC ADJUSTMENT 
8 HOURS INTEGRATION TIME 



Figure 5. The adjusted state for an initial 300 m/s eastward wind after 8 hr. The final zonal wind 
is close to 210 m/s which is nearly 30 m/s less than the predicted value based on a linear analysis. 
The oscillations which are evident in all three fields are due to gravity waves generated by the 
adjustment process. The waves propagate north and are reflected by the boundary condition at the 
pole. Therefore, the 300 m/s case is near the limit of validity of the numerical solution. 


201 



Figure 6. The adjusted state for an initial westward wind of 100 m/s is shown. The amplitude of 
the final zonal wind is approximately 90 m/s, or 10 m/s greater than expected based on a linear 
treatment. 
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Figure 7. Here the adjusted state for an initial westward wind of 300 m/s is shown. The final state 
also has a zonal wind of 300 m/s. The pressure gradient is very small compared to the cases 
shown in the previous figures, and there is almost no gravity wave energy generated, as can be 
seen in the smoothness and low values of the meridional wind. Thus, a significant asymmetry 
between the dusk and dawn sides of the auroral oval would be expected. We expect that on the 
dusk side, where the winds produced by the Lorentz acceleration are westward, less gravity wave 
energy will be produced by the adjustment process, and the pressure gradients will be smaller. 
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-100 m/s 
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-200 m/s 
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Figure 8. Table comparing the initial wind speed, the wind speed after adjustment in a linearized 
formulation, and the wind speed after adjustment when the nonlinear terms in the momentum 
equation are included. 
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SATELLITE ACCELEROMETER MEASUREMENTS OF NEUTRAL 
DENSITY AND WINDS DURING GEOMAGNETIC STORMS 

Frank A. Marcos 
Air Force Geophysics Laboratory 
Hanscom Air Force Base 
Bedford, MA 01731 
and 

Jeffrey M. Forbes 

Department of Electrical, Computer, and Systems Engineering 
Boston University 
Boston, MA 02215 

A new thermospheric wind measurement technique is reported which is based on a Satellite 
Electrostatic Triaxial Accelerometer (SETA) system capable of accurately measuring accelerations 
in the satellite’s in-track, cross-track and radial directions. Data obtained during two time periods 
are presented. The first data set describes cross-track winds measured between 170 and 210 km dur- 
ing a 5-day period (25-29 March 1979) of mostly high geomagnetic activity. In the second data set, 
cross-track winds and neutral densities from SETA and exospheric temperatures from the 
Millstone Hill incoherent scatter radar are examined during an isolated magnetic substorm occur- 
ring on 21 March 1979. 

A polar thermospheric wind circulation consisting of a two-cell horizontal convection pattern is 
reflected in both sets of cross-track acceleration measurements. Winds are generally trans-polar 
and parallel to the 1600h/0400 MLT meridian with return flows in the morning and evening sec- 
tors near 57° and 70° geomagnetic latitude, respectively. Winds are of order 100-200 ms^ 1 during 
quiet periods, and attain maximum speeds of order 300-600 ms 1 during storm conditions. The 
two-cell pattern is distinctively more well-ordered in geomagnetic rather than geographic coor- 
dinates. The substorm circulation pattern persists almost unattenuated for about 6 hours after the 
magnetic disturbance has returned to quiet levels. 

The density response is highly asymmetric with respect to its day/night behavior. At high 
geomagnetic latitudes (>60°) the daytime density variation, about a 40% increase from quiet 
levels, distinctively reflects the response to an increase in magnetic activity with a time delay of 
less than 3 hours. At lower latitudes the response is smaller (20%) and less well-defined, occurring 
with a time delay of about 6 ± 2 hours near 10-20° latitude. The nighttime density response, on 
the other hand, is not so well defined poleward as it is equatorward of 60 latitude. The time 
delay increases from about 2 ± 2 hours to 6 ± 2 hours from high to low latitudes. An exospheric 
temperature response occurs at Millstone Hill (42°N) with an amplitude of 210 K and time delay 
of between 0 and 2.5 hours. 

Latitude structures of the density response at successive times following the substorm peak sug- 
gest the equatorward propagation of a disturbance with a phase speed of between 300 and 600 
ms -1 . A deep depression in the density at high latitudes (>70°) is evident in conjunction with this 
phenomenon. The more efficient propagation of the disturbance to lower latitudes during the 
night is probably due to the “midnight surge” effect. 
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SATELLITE ACCELEROMETER MEASUREMENTS 



Figure 1. Cartoon showing orientation of accelerometer axes with respect to aerodynamic drag and cross-track wind vectors. 




Figure 2. Raw acceleration data from each accelerometer axis plotted as a function of time. Period 
covered is almost three full orbits. The top frame shows accelerations measured in the along-track 
(Z) direction. Larger accelerations are indicated by larger negative values. Thus satellite minimum 
altitudes (and maximum drags) occur near 80500 and 85900 seconds. 
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First Data Set: period 
of mostly high geomagnetic activity 

KP 



DAY (1979) 


Figure 3. Kp variation during a 5-day period during 1979. Cross-hatched and solid areas represent Quiet (Kp < 3°) and Active 
(Kp > 5-) periods, respectively. Vertical arrows denote periods of the first data set analyzed. 
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2400 


Figure 4. Schematic of 2-cell polar lower thermospheric circulation pattern in geomagnetic latitude 
and local time coordinates. Parallel lines represent bands of satellite tracks during 5-day period 
under study. 
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CROSS -TRACK WINDS 
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SETA-1 CROSS-TRACK WIND DATA 


Figure 5. Cross-axis winds, corresponding to geomagnetically disturbed conditions, plotted vs. GMT, geographic 
latitude (GGLAT), geomagnetic latitude (GMLAT), local time (LT), and altitude (ALT). Unsmoothed (light 
line) and smoothed (heavy line) data are shown. E, W, N refer, respectively to the fact that cross-track 
accelerations are due primarily to Eastward, Westward, and Northward flows. 
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Figure 6. (a) and (b): “cell centers” of double-vortex circulation system for quiet and active 
geomagnetic conditions, respectively, (c) and (d): schematics of 2-cell circulation patterns 
consistent with the quiet and active period data, respectively. 
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2400 


Figure 7. A comparison between areas including “cell centers” for quiet (black) and active 

(white) conditions. 
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AVERAGE WIND VELOCITY VS. KP 


Figure 8. Average wind data plotted vs. the corresponding average Kp values. Solid 
and open circles represent, respectively, measurements taken between 80-90 and 70-80 
GMLAT. Standard deviations of the means are also indicated. 




SUMMARY OF WINDS RESULTS FROM FIRST DATA SET 


1. SATELLITE ACCELEROMETER CROSS-TRACK WIND DATA ARE 
USED TO INFER LOWER THERMOSPHERE CIRCULATION PATTERNS 

2- MAIN FLOW IS ORIENTED PARALLEL TO 1600/0400 GLT MERIDIAN 

3- TWO-CELL SYSTEM COVERS GREATER AREA OF POLAR CAP DURING 
ACTIVE CONDITIONS 

4- WIND INTENSITY TENDS TO INCREASE WITH GEOMAGNETIC ACTIVITY 

Figure 9. Summary of wind results from first data set. 
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Second Data Set: 

isolated substorm 



Figure 10. Variation of the 3-hourly planetary magnetic index a p during the substorm under study. 
Vertical arrows indicate orbits (actually, times of polar crossings) which are explicitly included in 
the analysis of this paper. 
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CROSS-TRACK WINDS 



Figure 11. Cross-track winds from selected orbits (3, 4, 6, 7, 9, 10) as 
a function of geographic (top) and geomagnetic (bottom) latitudes. , 
“E” denotes east-ward oriented winds, “W” denotes westward-oriented 
winds, and TP indicates trans-polar oriented winds. The angle 0 
plotted as a dashed line in the top figure represents the relative 
orientation of the cross-track direction and the geographic meridian. 
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DENSITY ,x ;0 gem 



Figure 12. Temporal response of density at various geographic latitudes for day and night • 
crossings. Exospheric temperatures from Millstone Hill and the a p index are plotted in the bottom 
right hand frame of the figure. For reference, vertical lines are drawn through the “last quiet 
and “first maximum” values of a p . 
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GEOMAGNETIC LATITUDE 


Figure 13. Nighttime density responses for orbits 5, 6, 7, 8 as a function of geomagnetic 
latitude. Vertical arrows indicate approximate latitudes corresponding to density bulge 
associated with propagating disturbance launched during the substorm. 
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CONCLUSIONS 


1- SATELLITE ACCELEROMETER PROVIDES SIMULTANEOUS 
NEUTRAL DENSITY AND CROSS-TRACK WINDS 

2- TRANS-POLAR WINDS IN LOWER THERMOSPHERE ARE 
PARALLEL TO 1600/0400 GLT MERIDIAN 

3. WIND PATTERNS CAN PERSIST FOR SEVERAL HOURS 
AFTER STORM PERIOD 

4. NEUTRAL DENSITY DATA SHOW EVIDENCE OF 
"MIDNIGHT SURGE" EFFECT 

5- GRAVITY WAVE SIGNATURES ARE OBSERVED AT 
HIGH AND MIDDLE LATITUDES 

Figure 14. Conclusions from analysis of both data sets described in this 

paper. 
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TRANSFER FUNCTION ANALYSIS OF THERMOSPHERIC PERTURBATIONS 


H. G. Mayr, I. Harris, F. Varosi*, F. A. Herrero and N. W. Spencer 
Laboratory for Atmospheres, NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 


Applying perturbation theory, a spectral model (Mayr et al., JGR, 1984) in terms of vector 
spherical harmonics (Legendre polynomials, specifically) is used to describe the short term ther- 
mospheric perturbations originating in the auroral regions. The source may be Joule heating, par- 
ticle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is con- 
sidered, allowing for the collisional momentum exchange between species including Ar, 0 2 , N 2 , O, 
He and H. The coupled equations of energy, mass and momentum conservation are solved 
simultaneously for the major species N 2 and O. Applying homogeneous boundary conditions, the 
integration is carried out from the Earth’s surface up to 700 km. In our analysis, the spherical 
harmonics are treated as eigenfunctions, assuming that the Earth’s rotation (and prevailing circula- 
tion) do not significantly affect perturbations with periods which are typically much less than one 
day. Under these simplifying assumptions, and given a particular source distribution in the ver- 
tical, a two-dimensional transfer function is constructed to describe the three dimensional response 
of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this 
transfer function reveals five components: (A) The trapped component which is confined to the 
source region and decays slowly, with the time constant depending on the height of energy deposi- 
tion and the magnitude of eddy diffusion for example. (B) The quasi-horizontally propagating 
gravity wave which is represented by the lower cut-off and the first resonance maximum in the 
transfer function. In the thermosphere, it is also the dominant maximum. The horizontal wave 
length (order of 1000 km) and propagation velocity (700 m/s) are large. (C) The obliquely pro- 
pagating wave generated through partial reflection from the base of the thermosphere. This wave 
appears as a broad secondary maximum in the transfer function; its horizontal wave length and 
propagation velocity (about 350 m/s) is smaller; this wave is important near the source but can- 
not propagate very far horizontally. (D) The ducted wave which is produced in the non-dissipative 
lower atmosphere by total reflection from the Earth’s surface and partial reflection from the 
mesopause temperature minimum. Leaking back into the thermosphere where it originates, this 
wave has a relatively short wavelength but can travel large distances away from the source region 
(pole to equator); the horizontal propagation velocity is about 250 m/s. (E) The waves reflected 
from the surface which appear in the transfer function as broad secondary maxima. In general, 
these waves have relatively short horizontal wavelengths and dissipate rapidly above the 
mesopause. Thus, they cannot propagate far away from the source. To compile the transfer func- 
tion, the numerical computations are very time consuming (about 100 hours on a VAX for one 
particular vertical source distribution). However, given the transfer function, the atmospheric 
response in space and time (using Fourier integral representation) can be constructed with a few 

* Science Applications Research, Riverdale, Maryland 
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seconds of CPU. This model is applied in a case study of wind and temperature measurements on 
the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the 
auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in 
the polar region and then are reflected to propagate toward lower latitudes. 
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h(z) 



WAVE 

RESPONSE 


Figure 1. Block diagram illustrating the analysis. With a fixed height distribution of the excitation source h(z), 
and for frequencies a, and horizontal wave numbers i, the transfer function is constructed from a numerical 
integration of the perturbation equations (model). The model describes the perturbation amplitudes in a static 
background atmosphere with a height dependent temperature distribution, T(z), and corresponding density 
variations. The temporal and horizontal variations of the source are spectrally analyzed (source spectrum) and 
folded into the transfer function to construct (or synthesize) the wave response. 





axi-symmetric source-symmetric 



Figure 2. Assuming that the Earth’s rotation and prevailing circulation do not significantly affect thermospheric 
perturbations with periods much less than one day, one should expect that gravity waves tend to propagate radial 
ly away from a “source pixel” as illustrated in the shading. To first order, the response is “source symmetric”. 
Under this condition, spherical hamonics with zonal wave number m = 0 (Legendre polynomials) which normally 
describe axi-symmetric configurations, can be used to describe an arbitrary three dimensional configuration 
representing a conglomerate of individual source pixels. The two dimensional transfer function, F h (z, co), can 
then be used to describe the three dimensional wave response. 




TRANSFER FUNCTION 


TEMPERATURE 



Figure 3. For a wave period of half an hour, we show the computed transfer function of the 
relative temperature amplitude. The perturbations are excited by a height distributed heat source 
due to Joule heating. Adhering to the order in the abstract, the four propagating wave modes are 
identified: (B) the dominant, quasi-horizontally propagating gravity wave near the lower cutoff at 
the horizontal wave number H = 32, consistant with the theory of Hines (1960). (C) The wave 
reflected from the temperature gradient at the base of the thermosphere (Richmond, 1978). (D) 

The ducted wave returning from the lower atmosphere, which is formed through total reflection at 
the surface and partial reflection from the mesopause temperature minimum at 80 km (Hines, 

1960; Friedman, 1966); this wave appears in the transfer function as a narrow resonance maximum 
near & = 90 and persists throughout the atmosphere. (E) The broad secondary maxima at 
higher wave numbers which are due to reflection from the surface (Francis, 1974). 
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TRANSFER FUNCTION 



Figure 4. Computed transfer function for the temperature amplitude at 300 km, shown versus the 
frequency and the ratio between speed of sound and horizontal propagation velocity (proportional 
toX/co, for w> 2). The frequency range extends from 0 to 140 cycles per day, i.e., down to periods 
of 10 minutes. Note that the individual features of the transfer function line up along constant 
values ofj 6/co, consistent with classical gravity wave theory (Hines, 1960). The labels B through E 
identify the wave modes discussed in Figure 3. At the lower end of the frequency spectrum, label- 
led A, the transfer function shows the property of the thermospheric low pass filter (Volland and 
Mayr, 1970); the global scale and long term variations in the temperature (and density) are 
preferentially excited. This property is determined by the long time constant for dissipative pro- 
cesses such as heat conduction and accounts for the slow decay of the “trapped” density and 
temperature perturbations in the source region. 
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Figure 5. Analogous to Figure 4, we show the computed transfer functions for the vertical and 
horizontal velocities at 300 km. Note that the amplitudes are relatively small at lower frequencies; 
in contrast to the temperature (Figure 4), the low pass filter (labelled A) is not important. With in- 
creasing frequency, the amplitude of the vertical velocity grows and eventually becomes com- 
parable to that of the horizontal velocity. 
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DE-2 WATS 


DAY 81363 UT = 13 hrs. 
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Figure 6. Temperature and wind measurements from the DE-2 WATS experiment (Spencer et al., 
1982) during a magnetic substorm on December 25, 1981. (For quiet conditions, the temperatures 
are about 1100 K). Near 65° latitude, a region of energy deposition is indicated as inferred from 
the magnetic field and particle flux measurements on DE-2 (Farthing, Sugiura et al., 1981; Win- 
ningham et al., 1981). Airglow measurements of the 1500 A emission from DE-1 (Frank et al., 
1981) made four hours earlier showed the auroral oval centered near the magnetic pole. 
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Temperature C4T/T3 



Figure 8. During the transient process of energization, a broad spectrum of frequencies is prominent in the 
source. The frequency spectrum of the resultant thermospheric disturbance is determined by the source spectrum 
and the transfer function. Thus, the frequencies preferentially excited are those which match, through the 
resonances in the transfer function, the characteristic horizontal dimension of the source. This in turn produces 
the local ringing and the wave disturbances which propagate away from the source region. Shown here is a three 
dimensional display of the computed temperature perturbation plotted versus colatitude and time (at t = 0 the 
source is turned on). One perturbation propagates toward the equator, while a second one propagates toward the 
center, i.e. the magnetic pole, where convergence leads to amplification. The latter wave is “reflected” from the 
center and then also propagates equatorwards. Both wave trains are launched when the source is turned on and 
off, thus producing a complicated disturbance. 




COLATITUDE ( 0 R B - L 0 N G = 1 2 0 ) 

Figure 9. Along a satellite pass one would observe only a snapshot of this disturbance. Shown here are latitudinal 
cross sections of the computed vertical velocity and relative temperature perturbation about half an hour (24 
minutes (1) and 36 minutes (2)) after the source is turned on. In this time interval, the numerical results can 
simulate the observations from DE-2 that are shown in Figure 6. The height integrated (localized) energy source 
required to produce this perturbation is 160 erg/cm 2 /sec and is very large when compared with the EUV input. 
During magnetically disturbed conditions, however, such heating rates can be supplied through joule heating (M. 
Sugiura, this workshop). 




EQUATORIAL THERMOSPHERIC MEASUREMENTS OF TEMPERATURES 
AND WINDS AT AREQUIPA, PERU 


John W. Meriwether, Jr. 
Space Physics Research Laboratory 
University of Michigan 
Ann Arbor, MI 48103 

Manfred A. Biondi 
Department of Physics 
University of Pittsburgh 
Pittsburgh, PA 


Arequipa FPI observatory: Facts 

Field-widened 10 cm FPI Automatic operations. Located at Arequipa, Peru (16.5 S, 71.5 W) (at 
the NASA lidar satellite tracking station) 4.5 degrees south of magnetic equator FPI observed 
winds overhead of Jicamarca radar. Observing directions typically Z, N, S, W or N, S, E, W. 

Data selected from observations between new moon, April 1983, and new moon, August, 1983; 55 
nights reduced. Some results from 1984 but not all reduced. Error bars typically 15-20 m/s for 
winds, 50 to 75 degrees for temperature. 

Equatorial kinetic temperature: Results 

Enhancement of FPI temperatures above quiet levels a few hours after the start of magnetic ac- 
tivity. Magnitude of this enhancement about 300 degrees. This is followed by relaxation to pre- 
storm levels. Apparent average offset of RPI temperatures from MSIS by 100-200 degrees for 
relatively quiet times. Definite suggestion of a midnight thermal enhancement for April and August 
data. Magnitude about 100 degrees. Seen in both 1983 and 1984 observations. 

Conclusions 

6300A surface brightnesses: Results 

Definite enhancement of 6300A surface brightness in the south as compared with other direc- 
tions. Probably connected to the tropical airglow arcs. Meridional winds small (<25 m/s) 
throughout night. Indication of northward migration of the observed 6300A enhancement in the 
evening hours as observations approach local winter solstice. This is probably related to the observ- 
ed poleward (to the south) meridional wind (of magnitude 50 m/s) in this period. 

Equatorial thermospheric winds: Results 

Zonal component of winds always eastward, but speed approaches zero sooner near equinox 
than at summer solstice. Typical magnitude at peak is of order 100 to 150 m/s. Suggesstion of 
zonal wind increase after twilight and recovery of 6300A emission for April data. Origin not clear 
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but may be related to midnight thermal enhancement. Meridional wind virtually zero for equinox 
in 1983; shows evening flow towards winter hemisphere in early evening for solstice data. Sugges- 
tion of post-midnight surge in April 1984 data. No major effects associated with mag. storm ac- 
tivity. Suggestion of decrease in zonal component below nominal levels. 
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Hernandez: Mid-Latitude Thermospheric Neutjual Temperatures 



EXPERIMENTAL TEMPERATURE (K) 

Fig. 3. Comparison of the experimental temperatures with the 
MSIS empirical model [Hedin ft ui., 1977J. 


EXPERIMENTAL TEMPERATURE (K) 

Fig. 4. Comparison of the experimental temperature with the 
empirical model of 77W//i>r ft ai (1977a, fe|. 



EXPERIMENTAL TEMPERATURE IK ) 

Fig. 11. Fit of the experimental data to four parameters:. solar 
radio flux, geomagnetic activity, annual variation of the solar 
declination and a semt-annuaJ variation. The correlation coefficient 
is 0.91. 


EXPERIMENTAL TEMPERATURE tX) 

Fig. 5. Comparison of the experimental temperatures with the 
empirical model of Jacchiu ( 19771. 
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POLAR CAP OBSERVATIONS OF THERMOSPHERIC WINDS AND TEMPERATURES 

AT SONDRE STROMFJORD, GREENLAND 

John W. Meriwether, Jr. and P. Shih 
Space Physics Research Laboratory 
University of Michigan 
Ann Arbor, MI 48105 

Conclusions 

Temperatures: ... 

Agreement of averaged temperatures with MSIS looked reasonable for several nights, but tor 
many nights there are differences: 1) midnight period of cooling and 2) temperature increases 
associated with overhead crossings of the auroral belt. The observed rise of the temperature before 
dawn in conjunction with the high 6300A intensities suggests a connection beteween the two ef- 
fects: soft particle precipitation most likely candidate but frictional heating perhaps also a possibili- 
ty. Comparison with TGCM calculations also needed. 

Thermospheric winds: , 

The technique for formulating neutral wind vectors performs well in most cases. The observed 
patterns show evidence for abatement in the midnight sector in the meridional wind component at 
the separatrix between the two cells with a frequency of the order of 20-25 %. Also observed in 
radar observations at Sondre Stromfjord. The observed patterns for magnetically quiet conditions 
show flow characteristic of the auroral belt, westward in evening followed by the midnight surge. 
The observed patterns for active conditions show dominance either by the evening cell or the morn- 
ing cell, but most often the former. We need to correlate our results with IMF measurements. 


253 







Meridional Wind (m/s) 













Intensity (kR) Vertical Winds (m/s) Temperatures (k) 



MAY 3, 1984 


Figure 5 


258 



























ZONAL WIND OBSERVATIONS DURING A GEOMAGNETIC STORM 

N. J. Miller and N. W. Spencer 
Laboratory for Atmospheres 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 

In situ measurements taken by the Wind and Temperature Spectrometer (WATS) onboard the 
Dynamics Explorer 2 spacecraft during a geomagnetic storm display zonal wind velocities that are 
reduced in the corotational direction as the storm intensifies. The data were taken during 
November 24-26, 1982 within the altitudes 275-475 km in the dusk local time sector equatorward 
of the auroral region. Characteristic variations in the value of the Dst index of horizontal 
geomagnetic field strength are used to monitor the storm evolution. The detected global rise in at- 
mospheric gas temperatures indicates the development of thermospheric heating. Concurrent with 
that heating, reductions in corotational wind velocities were measured equatorward of the auroral 
region. Just after the sudden commencement, while thermospheric heating is intense in both 
hemispheres, eastward wind velocities in the northern hemisphere show reductions ranging from 
500 m/s over high latitudes to 30 m/s over the geomagnetic equator. After 10 hours storm time, 
while northern thermospheric heating is diminishing, wind velocity reductions, distinct from those 
initially observed, begin to develop over southern latitudes. In the latter case, velocity reductions 
range from 300 m/s over the highest southern latitudes to 150 m/s over the geomagnetic equator 
and extend into the northern hemisphere. The observations highlight the interhemispheric asym- 
metry in the development of storm effects detected as enhanced gas temperatures and reduced 
eastward wind velocities. Zonal wind reductions over high latitudes can be attributed to the storm 
induced equatorward spread of westward polar cap plasma convection and the resulting plasma- 
neutral collisions. However, those collisions are less significant over low latitudes; so zonal wind 
reductions over low latitudes must be attributed to an equatorward extension of a thermospheric 
circulation pattern disrupted by high latitude collisions between neutrals transported via eastward 
winds and ions convecting westward. 
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1 1 1 1 1 I I I 

STORM TIME -10 0 10 20 30 40 50 60 HOURS 

Figure 1. Magnetic characteristics of the period November 24-26, 1982. The Dst index decreases to 
-21 87 at 1800 UT on the 24th; the three-hourly Kp index reaches 7 near 1400 UT on the 24th and 
remains there until after 0300 UT on the 25th. Storm time is measured from the time of the 
sudden commencement. 
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60° N 



STORMTIME, HOURS 

Figure 2. Gas temperatures measured by the WATS instrument in the dusk local time sector within 
the altitude range 275-475 km. The effects of thermospheric heating appear between 0-20 hours 
storm time over all dusk latitudes equatorward of 60° invariant latitude. The temperature rises 
660° K over 60 °N invariant latitude but only 70° K over ghe geomagnetic equator. 
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STORM TIME, HOURS 


Figure 3. Zonal wind velocities measured over northern latitudes by the WATS instrument at the 
same time as the measurements represented in Figure 2. Velocities are relative to the corotation 
velocity, hence zero means corotation. Over all northern latitudes equatorward of 60° invariant 
latitude, velocities display a decrease in value immediately after the sudden commencement. 
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Figure 4. Zonal wind velocities measured as in Figure 3 but over southern latitudes. After 20 hours 
storm time the data display a velocity decrease that affects the zonal wind over latitudes from 60° 
south invariant latitude to 20° north invariant latitude. 
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MEASUREMENTS OF THERMOSPHERIC RESPONSE TO AURORAL ACTIVITIES 


S. Okano* and J. S. Kim 

Department of Atmospheric Science 
State University of New York at Albany 
1400 Washington Avenue 
Albany, New York 12222 

The Joule heating produced by auroral electrojets and its thermospheric response can be studied 
by monitoring the thermospheric temperatures by optical methods; simultaneously, the concurrent 
auroral electrojet activities can be investigated by using geomagnetic records obtained from sta- 
tions along a meridian close to the observation site of optical measurements. 

We report, in this paper, the measurements of thermospheric response to auroral activities 
which were made at Albany (42.68 °N, 73.82 °W), New York on September 2, 1978 (UT) when an 
isolated substorm occured. The thermospheric temperatures were measured by using a high- 
resolution Fabry-Perot interferometer that determines the line profiles of the [OI] 6300A line 
emission. The intensities and latitudinal positions of auroral electrojets were obtained by the 
analysis of magnetograms from the IMS Fort Churchill meridian chain stations. 

* Present Address: Upper Atmosphere and Space Research Laboratory, Tohoku University, 
Sendai, Japan 
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Figure 2a. The plots of the one-minute averages of X m (magnetic north) components 
observed at stations in the Fort Churchill chain on September 2, 1978 (UT). 
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Figure 2b. The plots of the one-minute averages of Z (downward) components observed at 
stations in the Fort Churchill chain on September 2, 1978 (UT). 
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AXm,AZ (V) 

Figure 3. AX m and AZ values at 0720 UT on September 2, 1978 from each station were 
plotted and the best-fitted curves were drawn through AX m and AZ data, respectively, in order 
to determine the position and intensity of the westward electrojet. 
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Figure 4. The position of electrojet and (AX m max ) 2 value that were obtained by the analysis of magnetograms 
from the IMS Fort Churchill meridian chain stations and thermospheric temperature measured from Albany, New 
York are plotted as a function of universal time from top to bottom, respectively. 


INITIAL RESULTS OF THE GLOBAL THERMOSPHERIC MAPPING STUDY (GTMS) 

W. L. Oliver, J. E. Salah, R. G. Musgrove, and J. M. Holt 
MIT Haystack Observatory, Westford, Massachusetts 01886 U.S.A. 

V. B. Wickwar 

SRI International, Menlo Park, California 94205 U.S.A. 

G. J. Hernandez 

NOAA, Boulder, Colorado 80303 U.S.A. 

R. G. Roble 

NCAR, Boulder, Colorado 80307 U.S.A. 

The Global Thermospheric Mapping Study (GTMS) is a multi-technique experimental study of 
the thermosphere designed to map simultaneously its spatial and temporal morphology with a 
thoroughness and diversity of measurement techniques heretofore unachieved. Three-day cam- 
paigns at the summer and winter solstices are planned to study the seasonal variations in ther- 
mospheric structure. The GTMS is designed around the Incoherent Scatter Radar Chain in the 
western hemisphere. The European incoherent scatter radars and the worldwide communities of 
Fabry Perot interferometers, meteor wind radars, partial reflection drifts radars, MST radars, and 
satellite probes are included to extend the spatial coverage and types of measurements available. 
Theoretical and modeling support in the areas of thermospheric and ionospheric structure, tides, 
and electric fields are included to aid in program planning and data interpretation. 

The initial GTMS campaign was conducted on 26-28 June 1984, during a period of sharply 
declining solar activity. Solar activity was low on the three observation days (F10.7 = 97, 98, 96) 
and magnetic conditions were unsettled to active (A= 10, 12, 20). All six incoherent scatter radar 
facilities collected data. Each collected F region data day and night while Saint Santin and 
Millstone Hill additionally collected E region data during daylight hours. Initial results from Son- 
drestrom and Millstone Hill are presented. Good quality Fabry Perot data were collected at Fritz 
Peak and San Jose dos Campos. Weather conditions produced poor results at Arequipa and 
Arecibo. Initial results from Fritz Peak are presented. Mesosphere/lower-thermosphere observa- 
tions were conducted under the ATMAP organization. The magnetometer chains also were opera- 
tional during this campaign. Initial TGCM predictions have been made for assumed solar- 
geophysical conditions, and selected results are presented. 

Sondrestrom ionospheric results show a decrease in electron density from June 26 to 27 to 28, 
isolated periods of ion frictional heating during the nights throughout the field of view, and strong 
ionospheric convection velocity reversals well beyond midnight. Neutral wind results show a major 
period of enhanced northward velocity in response to strong heating to the south of the radar on 
June 28, the most disturbed of the three days. 

Millstone Hill ionospheric results show lower densities on June 27 than on June 26 but higher 
densities on June 28. Frictional heating effects and velocity reversals are restricted to the most 
northerly observations. Millstone Hill neutral atmosphere results show a higher temperature and 
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atomic oxygen density over a wide latitude range on June 28 than on the previous two days, with 
particularly strong effects at lower latitudes. These are believed to be the first incoherent scatter 
radar determinations of the latitudinal variation of the neutral atomic oxygen density. 

Fritz Peak Fabry Perot data show a strong latitudinal gradient in neutral temperature, with 
temperatures 5 degrees in latitude to the north of the station exceeding those to the south by 
200 K. Early on the morning of June 28 the Fabry Perot winds showed a sudden and strong shift 
toward the equator with the winds to the north exceeding 300 m/s but the winds to the south at- 
taining a speed of only 130 m/s. 

TGCM simulations were made to predict the global steady-state behavior of thermospheric 
temperature, winds, density, and composition for pre-assumed values for solar flux (F10.7 = 120- 
somewhat high) and magnetic activity (Ap = 4 - somewhat low). A cross-tail magnetospheric 
potential and an auroral particle input of 2 ergs/cm **2/sec were included but tides were not. 
These simulations predict a strong summer-solstice neutral temperature gradient over North 
America with a temperature maximum over Greenland. Meridional wind predictions over Son- 
drestrom show a strong diurnal flow with poleward velocities of 170 m/s during the day and 

equatorward velocities of 300 m/s at night. A large quantity of such results are available from the 
simulations. 

As more data are received and analyzed, composite maps of global upper atmospheric and 
ionospheric structure will be compiled and compared with models and simulations which will be 
updated as pertinent information on solar/geophysical inputs becomes available. Magnetometer 
data will be used to estimate global Joule heating rates to act as inputs to the models. The models 
will rely to various degrees on inputs from the measurements and on feedback from other models 
as thermospheric, ionospheric, and electrodynamic properties of the upper atmosphere all are in- 
terdependent. 

A second GTMS campaign will be conducted on 15-17 January 1985. Optical coverage, especial- 
ly at high northern latitudes near the Radar Chain, should be much improved during this period. 

The GTMS data-taking campaigns will be followed by a Workshop at MIT on 17-19 July 1985. 
In addition to the GTMS solstice campaigns, the Workshop will include the Equinox Transition 
campaign, a radar/optical campaign coordinated by AFGL in September 1984, designed to 
observe the summer-winter thermosphere circulation transition and the competition between high- 
latitude and sub-solar forcing in controlling circulation structure. Consideration will be given to 
presenting these results at an AGU session in December 1985 or at an NCAR/URSI Radar 
Workshop on coordinated radar campaigns. 
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GTMS - EXPERIMENTAL ORGANIZATION 


OVERALL COORDINATION - MILLSTONE HILL (W. OLIVER, J. SALAH) 


TECHNIQUE COORDINATORS - 

INCOHERENT SCATTER RADARS 
SONDRE STROMFJORD 
MILLSTONE HILL 
ARECIBO, JICAMARCA 
EISCAT , SAINT SANTIN 
FABRY PEROT INTERFEROMETERS 
METEOR RADARS 
PARTIAL REFLECTION RADARS 
MST RADARS 
SATELLITE PROBES 


W. OLIVER 
V. WICKWAR 

V. OLIVER 

W. SWARTZ 
D. ALCAYDE 

J. MERIWETHER 
R. ROPER (GLOBMET) 
J. FORBES (ATMAP) 
J. FORBES (ATMAP) 
A. HEDIN 


GTMS - THEORY/MODEL/SIMULATION ORGANIZATION 

1 - THERMOSPHERE SIMULATION (TGCM) ROBLE ET AL. 

2 - THERMOSPHERE MODEL COMPARISON (E.G., MSIS) HEDIN ET AL. 

3 - TIDES MODELING/COMPARISON FORBES ET AL. 

4 - IONOSPHERE SIMULATION SCHUNK ET AL. 

5 - ELECTRIC FIELD MODELING/COMPARISON RICHMOND ET AL. 

6 - NE/TE MODEL COMPARISON BRACE ET AL. 

7 - NE MODEL COMPARISON CHIU ET AL. 
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GTMS 


26-28 JUNE 1984 


SOLAR/GEOPHYSICAL INDICES 


MONTHLY 

AVERAGES 




DAILY VALUES 


MONTH 

PI 0 . 7 

DAY 

FI 0. 7 

APR 

K-FREDERICKSBURG 

K-ANCHORAGE 

MARCH 

1 22.0 

JUNE 24 

100.3 

16 

5 -4 -2-2-2 -2-2-3 

5-6-3-3-2-3-2-3 

APRIL 

1 28.7 

25 

101 .3 

13 

4-3-2-2-2-2-3-3 

4-3-3-2-3-3-3-2 

MAY 

1 28.3 

26 

96.8 

10 

3-3-3-2-1 -1-2-3 

3-3-2-3-2-3-2-3 

JUNE 

100.3 

27 

98.2 

12 

4-4-3-2-1 -0-2-2 

5-4-4-4-3-2-3-2 

JULY 

89-3 

28 

96.2 

20 

3-3-4-3-3-3-4-4 

4-4-5-4-4-6-4-3 

AUGUST 

83-7 

29 

97-0 

12 

3-2-2-3-1 -3-3-3 

4-3-3-3-1 -3-3-3 

SEPT. 

(78.1 ) 

30 

97-8 

1 1 

3-3-2-2-2-2-3-3 

3-4 -3 -3 -4 -4-2-3 


SITES INVOLVED IN THE GLOBAL THERMOSPHERIC MAPPING STUDY 


SITE 


STATE/COUNTRY LAT LON EXPERIMENT TYPE 


Saint Santin 

SOUCY 

Bologna 

EISCAT 

Kiruna 

Kyoto 

Adelaide 

Christchurch 

Pai rbanks 

Poker Plat 

Calgary 

Saskatoon 

Fritz Peak 

Sunset 

Urbana 

Atlanta 

Ann Arbor 

Laurel Ridge 

Jicamarca 

Millstone Hill 

Arequipa 

Durham 

Arec ibo 

Sondr estrom 

San Jose dos Campos 


Prance 44 
West Germany 52 
Italy 44 
Scandanavia 69 
Sweden 68 
Japan 35 
South Australia -35 
Hew Zealand -43 
Alaska 64 
Alaska 64 
Alberta 51 
Saskachewan 52 
Colorado 40 
Colorado 40 
Illinois 40 
Georgia 33 
Michigan 42 
Pennsylvania 40 
Peru -12 
Massachusetts 43 
Peru -1 6 
New Hampshire 43 
Puerto Rico 18 
Greenland 67 
Brazil -23 


2 

IS 





1 0 





MST 

1 1 



MR 



19 

IS 





20 


PP 




135 



MR 



1 38 


PP 


PR 


1 72 




PR 


21 2 


PP 




212 





MST 

246 


pp 




254 




PR 


255 


PP 




255 





MST 

272 



MR 


MST 

276 



MR 



277 


PP 




281 


PP 




283 

IS 




MST 

289 

IS 

PP 




289 


PP 




289 



MR 



294 

IS 

PP 



MST 

309 

IS 

PP 




314 


PP 





LAT: LATitude, + in northern hemisphere, - in southern 

LON: east LONgitude 

IS: Incoherent Scatter 

PP: Fabry Perot interferometer 

MR: Meteor wind Radar 

PR: Partial Reflection drifts radar 

MST: Mesosphere-Stratosphere-Troposphere radar 
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GTMS 


26-28 JUNE 1984 


IS 


EISCAT 


DAYS OF OBSERVATION: 26-28 JUNE 
OBSERVING SCHEME: CP-3 - WIDE LATITUDE SCAN 
CYCLE TIME: 30 MINUTES 

MEASURED QUANTITIES - IONOSPHERIC (F REGION) 

NE, TE, TI, E FIELD - VS LATITUDE 

MEASURED QUANTITIES - NEUTRAL (UPPER THERMOSPHERE) 
TEMPERATURE - VS LATITUDE 
MERIDIONAL WIND - VS LATITUDE 
OXYGEN DENSITY - VS LATITUDE 

DATA QUALITY: PRESUMED GOOD 


GTMS - 26-28 JUNE 1984 - IS - SAINT SANTIN 


DAYS OF OBSERVATION: 26-28 JUNE- 

OBSERVING SCHEME: VERTICAL TRANSMITTED BEAM (SAINT SANTIN) 

3 RECEIVERS SCANNING TRANSMITTED BEAM 
(NANCAY, MENDE, MONPAZIER ) 

(MONPAZIER RESTRICTED TO F REGION ON 28-TH) 

CYCLE TIME: DAY - 60 MINUTES 
NIGHT - 30 MINUTES 

MEASURED QUANTITIES - IONOSPHERIC (LOCAL) 

F REGION - NE, TE, TI , E FIELD 
E REGION - NE, TE, TI , COLL. FREQ., ION DRIFT 

MEASURED QUANTITIES - NEUTRAL (LOCAL) 

F REGION - TEMPERATURE 

MERIDIONAL WIND 
OXYGEN DENSITY (DAY) 

E REGION - TEMPERATURE (DAY) 

WIND VECTOR (DAY) 

N2 DENSITY (DAY) 

DATA QUALITY: GOOD 
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GTMS 


26-28 JUNE 1984 


IS 


ARECIBO 


DAYS OP OBSERVATION - 26-28 JUNE 

OBSERVING SCHEME - CONICAL BEAM SWING 

INTERLACED POWER AND SPECTRAL MEASUREMENTS 

CYCLE TIME: 16 MINUTES FOR ANTENNA SWING 

1-2 MINUTES FOR LINE-OF-SIGHT MEASUREMENTS 

MEASURED QUANTITIES - IONOSPHERIC (LOCAL P REGION) 

NE, TE, TI, E FIELD 

MEASURED QUANTITIES - NEUTRAL (LOCAL UPPER THERMOSPHERE) 
TEMPERATURE 
MERIDIONAL WIND 
OXYGEN DENSITY (DAYTIME) 

DATA QUALITY - GOOD 


GTMS - 26-28 JUNE 1984 - 

DAYS OP OBSERVATION: FARADAY ROTATION 

IS SPECTRAL DATA 


IS 

- JICAMARCA 

(NE) 

- 26-28 

JUNE 


- 27-28 

JUNE 


OBSERVING SCHEME - 

SPLIT BEAM, EAST AND WEST 

CYCLE TIME: 5 MINUTES 

MEASURED QUANTITIES - IONOSPHERIC (LOCAL F REGION) 

NE, TE, TI, E FIELD 

MEASURED QUANTITIES - NEUTRAL (LOCAL UPPER THERMOSPHERE) 
TEMPERATURE 

?? OXYGEN DENSITY (DAYTIME) 

DATA QUALITY: FARADAY ROTATION - GOOD 

IS SPECTRAL - UNCERTAIN 
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GTMS 


26-28 JUNE 1984 


SONDRESTROF 


DAYS OF OBSERVATION - 26-28 JUNE 
OBSERVING SCHEME - 

11 -POSITION MEASUREMENT: 5 PAIRS STRADDLING FIELD LINE 

(FOR E FIELD VS LATITUDE) 
1 POSITION UP B FIELD 
(FOR LOCAL WINDS'! 

N-S ELEVATION SCAN (SCALAR QUANTITIES VS LATITUDE') 

UP B FIELD AFTER 11 -POSITION AND ELEVATION SCAN 

CYCLE TIME: 22-5 MINUTES (6 MINUTES UP B) 

MEASURED QUANTITIES - IONOSPHERIC (F REGION) 

NE, TE, T±, E FIELD - VS LATITUDE 

MEASURED QUANTITIES - NEUTRAL (UPPER THERMSOPI1ERE ) 
TEMPERATURE - VS LATITUDE 

MERIDIONAL WIND - LOCAL 
?? OXYGEN DENSITY 

DATA QUALITY - GOOD 


GTMS - 26-28 JUNE 1984 - IS - MILLSTONE HILL 


DAYS OF OBSERVATION - 26-28 JUNE 


OBSERVING SCHEME - 

DOUBLE ELEVATION SCAN: ONE ALONG MILLSTONE MERIDIAN 

ONE CANTED 30 DEGREES TO WEST 

ZENITH MEASUREMENTS 

DAYTIME: ADDED 3-POSITION E REGION MEASUREMENT 

CYCLE TIME: NIGHT - 30 MINUTES 

DAY - 60 MINUTES 


MEASURED QUANTITIES - IONOSPHERIC 

F REGION - NE, TE, TI, E FIELD - VS LATITUDE 
E REGION - NE, TE, TI, COLL. FREQ., ION DRIFT - LOCAL 


MEASURED QUANTITIES - NEUTRAL 
F REGION - TEMPERATURE 

MERIDIONAL WIND 
OXYGEN DENSITY 
E REGION - TEMPERATURE 
WIND VECTOR 
N2 DENSITY 


VS LATITUDE 
VS LATITUDE 
VS LATITUDE (DAY) 
LOCAL (DAY) 

LOCAL (DAY) 

LOCAL (DAY) 


DATA QUALITY - GOOD 
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GTMS 


26-28 JUNE 1984 


FABRY PEROT OBSERVATIONS 


6300-A NIGHTTIME OBSERVATIONS 


1 


2 


SAO JOSE DOS 
JUNE 25 
JUNE 26/27 
JUNE 27 
JUNE 28 
JULY 01/02 
JULY 02/03 
JULY 03/04 


CAMPOS: 

- 1850-2330 LT 

- 1900-0430 LT (VERY LOW INTENSITY AFTER 00 LT ) 

- NO OBSERAVTIONS (BAD WEATHER) 

- NO OBSERVATIONS (BAD WEATHER) 

- 1950-0200 LT 

- 1830-0400 LT (LOW INTENSITY) 

- 1830-0300 LT 


AREQUIPA - DATA TAKEN, WEATHER WAS CLOUDY 


3 - ARECIBO - DATA TAKEN, SKY WAS HAZY, POOR QUALITY 

4 - FRITZ PEAK: 

JUNE 27 - 0300-1100 UT (TEMPERATURE ONLY) 

JUNE 28 - 0300-1100 UT (TEMPERATURE + HORIZONTAL WIND) 

5 - LAUREL RIDGE - TAPE RECORDER MALFUNCTION 


6 - ANN ARBOR - HO DATA 


7 - CALGARY - NO DATA 

8 - FAIRBANKS - NO DATA, INOPERATIVE IN SUMMER 

9 - SONDRESTROM - NO DATA, INOPERATIVE IN SUMMER 


OTHERS ?? 


GTMS - 26-28 JUNE 1984 - D&E-REGION /MESOSPHERE OBSERVATIONS 


1 - ENTIRE ATMAP NETWORK - REPORTS INCOMPLETE 

METEOR RADARS 

PARTIAL REFLECTION DRIFTS RADARS 
MST RADARS 

2 - INCOHERENT SCATTER RADARS - DAYTIME ONLY 

MILLSTONE HILL 
SAINT SANT IN 

3 - FABRY PEROT INTERFEROMETERS (GREEN LINE) - NO OBSERVATIONS 

(ONLY F REGION RED LINE DATA WERE TAKEN) 
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GTMS 


26-28 JUNE 1984 


SATELLITE DATA 


??? 


GTMS - 26-28 JUNE 1984 - MAGNETOMETER DATA 


MAGNETOMETER CHAINS WERE OPERATING 

KAMIDE AND AKASOPU PLAN TO CALCULATE GLOBAL MAPS OP 
JOULE HEAT PRODUCTION AND IONOSPHERIC CURRENTS 
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GTMS/EQUINOX-C IRCULATION WORKSHOP 


WORKSHOP TO STUDY GLOBAL FEATURES OF 
THERMOSPHERE STRUCTURE AS OBSERVED DURING: 

( 1 ) GTMS SUMMER (JUNE 1984 ) AND WINTER (JANUARY 1985 ) 
SOLSTICE CAMPAIGNS 

( 2 ) EQUINOX-CIRCULATION (SEPTEMBER 1984 ) CAMPAIGN 

WORKSHOP DATES/PLACE: 17-19 JULY 1985 AT MIT 


GTMS/EQUINOX-C IRCULATION AGU SESSION 

DECEMBER 1985 MEETING IN SAN FRANCISCO 
SESSION TO BE PLANNED AT JULY 1985 WORKSHOP 
ABSTRACTS DUE IN TO AGU IN SEPTEMBER 1985 
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Illustration of the latitude-altitude coverage attained by a 
north-south elevation scan down to about 5 degrees at Millstone 
Hill. Plotted is the logarithm (base 10) of the electron density 

(cm**-3). 

MILLSTONE HILL RADAR - 27,JUN,1984 
START UT = 0:24:44 AZ = 0.1 EL = 5.3 

END UT = 0:41: 3 AZ = 180 .1 EL = 5.7 
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DEGREES 


Thermospheric temperatures measured in four directions from Fritz 
Peak during the night of 26/27 June 1984 - 


F.P.O. 6300 A TEMPERATURES 840627 (DAY 179) 



300 . 500 . 700 . 900 . 1100 . 
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MEflIOIONAu 'WIND < M/SCC I 


MERIDIONAL WINO (rt/SEC) f 62.5. -50.0) 



TGCM simulation of the steady-state meridional wind flow at F 
region altitudes above Sondre Stromfjord during 26-28 June 1984. 
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OBSERVATIONS OF VERTICAL WINDS AND 
THE ORIGIN OF THERMOSPHERIC GRAVITY WAVES 

LAUNCHED BY AURORAL SUBSTORMS AND WESTWARD TRAVELLING SURGES 

David Rees 

Department of Physics and Astronomy 
University College London 
Gower Street, London WC1E 6BT, UK 

Several sequences of observations of strong vertical winds in the upper thermosphere will be 
discussed in this paper, in conjunction with models of the generation of such winds. In the auroral 
oval, the strongest upward winds are observed in or close to regions of intense auroral precipita- 
tion and strong ionospheric currents. The strongest winds, of the order of 100 to 200 m/sec are 
usually upward, and are both localized and of relatively short duration (10 to 20 min). In regions 
adjacent to those displaying strong upward winds, and following periods of upward winds, 
downward winds of rather lower magnitude (40 to about 80 m/sec) may be observed. Strong and 
rapid changes of horizontal winds are correlated with these rapid vertical wind variations. Con- 
sidered from a large scale viewpoint, this class of strongly time-dependent winds propagate global- 
ly, and may be considered to be gravity waves launched from an auroral source. During periods of 
very disturbed geomagnetic activity, there may be regions within and close to the auroral oval 
where systematic vertical winds (upward or downward) of the order of 50 m/sec will occur for 
periods of several hours. Such persistent winds are part of a very strong large scale horizontal 
wind circulation set up in the polar regions during a major geomagnetic disturbance. This second 
class of strong horizontal and vertical winds corresponds more to a standing wave than to a gravi- 
ty wave, and it is not as effective as the first class in generating large scale propagating gravity 
waves and correlated horizontal and vertical oscillations. A third class of significant (10 to 30 
m/sec) vertical winds can be associated with systematic features of the average geomagnetic energy 
and momentum input to the polar thermosphere, and appear in statistical studies of the average 
vertical wind as a function of Universal Time at a given location. 

INTRODUCTION 

Recent papers by Hernandez 1 , Meriwether et al 2 , and Rees et al 3 ’ 4 have described the observa- 
tion of significant (20 to 150 m/sec) vertical winds from mid-latitude 1 and from high (auroral) 
latitude stations 2 ’ 3 > 4 , using ground-based Fabry-perot interferometers to observe the ther- 
mospheric OI 630 nm emission line. Spencer et al 5> 6 have observed similar vertical wind 
magnitudes in-situ from their satellite-borne WATS (Wind and Temperature Sensor) instruments 
on Atmospheric Explorer and on Dynamics Explorer. 

Such wind magnitudes are at least one order of magnitude larger than the systematic vertical 
winds (3 -5 m/sec) which can be associated with the diurnal ‘breathing’ of the thermosphere 7 ’ 8 at 
upper thermospheric levels (300 to 400 km) in response to solar UV and EUV heating at middle 
and low latitudes. Generally, the observations of strong vertical winds indicate that there is a rapid 
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time-variation of such winds, and that the horizontal scale of the regions of strongest upward and 
downward flow is also relatively limited, of the order of 100 to a few hundred km. 

The observations indicate that the strong vertical winds are the result of intense excitation of the 
middle and upper thermosphere. The source is strong, localized and may vary rapidly with time. 
Correlated horizontal and vertical wind changes 1 - 3 suggest that, after the initial excitation, the 
propagation of the wind disturbance follows the natural thermospheric response to strong excita- 
tion, a form of gravity wave propagation, where the disturbances may propagate throughout the 
entire global thermosphere. 

RESUME OF OBSERVATIONS 

The most frequently observed class of strong (> 50 m/sec) vertical wind disturbance 3 is 
associated with auroral substorms and westward-travelling-surges, and similar energetic auroral 
phenomena. On examining the data from several hundred nights when observations of ther- 
mospheric winds are available from Kiruna Geophysical Institute (67 N, 22 E geographic), a loca- 
tion within the auroral oval, it would appear that such disturbances are observed, on average, 
about one in two or three nights of clear sky observations, for the geomagnetic activity conditions 
of the period late 1980 to late 1983. 

The observations indicate that when any substorm or westward travelling surge occurs equator- 
ward of the observing station, or propagates from equatorward to overhead, strong upward ver- 
tical winds will be observed immediately prior to the arrival of the surge or substorm, as seen in 
the auroral emissions (particularly in OI 630 nm). Associated with the strong upward winds is an 
outward explosion of the horizontal winds. 

There is a high degree of correlation between the vertical and horizontal wind changes, which 
will be described in reference to three particular examples. Events which occur poleward of the 
observing station create similar changes of the horizontal wind, and although the vertical wind 
changes are not usually observed as clearly for such events, this has to be purely a result of less 
favourable observing geometry. 

I. Nov 23/24 1982 

Figure 1 shows the combined horizontal and vertical wind observations from KGI during the 
night of Nov 23/24 1982. There was a major substorm and westward travelling surge about 19 
UT, which propagated toward and overhead Kiruna from the southeast. This propagation can be 
inferred from the OI 630 nm intensity observations. Following a general increase in the magnitude 
of westward and equatorward winds between 18.00 UT and about 18:40 UT, there was an short- 
lived period of extremely strong upward winds (up to 165 m/sec) starting about 18.50 UT, which 
coincided with a sudden reversal in the meridional wind as observed to the north of Kiruna. 

Note: all of the observations of the ‘horizontal’ winds were made at a zenith distance of 60 
degrees. For a mean emission altitude of 240 km, this observing region is some 400 km away from 
Kiruna along the specific viewing direction. 

After 19.20, the strong upward vertical winds subsided, and the meridional wind, to the north 
of Kiruna, reverted to a strong equatorward value (400 m/sec). The degree of correlation between 
changes in the various wind components, and the nature of the ‘outward explosion’ of the winds 
during the disturbance is illustrated further in Figure 2. There were at least two surges of the 
equatorward wind observed to the south of Kiruna, during the period 19.00 to 20.30 UT, which 
could be associated with this period of intense auroral activity. 
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Later in the night, up to midnight UT, there were several significant fluctuations of the vertical 
wind, of the order of 50 m/sec. These fluctuations are difficult to relate to local auroral activity, 
and may reflect the propagation of waves from more distant auroral disturbances. 

II. Feb. 4/5 1983 

Wind observations during a second disturbed night (Feb 4/5 1983) are shown in Figure 3. There 
were several fluctuations of the vertical wind up to about 20 UT, and then there was a short 
period of violent (150 m/sec) upward wind centered about 21 UT. Immediately before and after 
the period of strongest upward winds, there was a short period of 50 m/sec downward winds. As 
with the wind observations of the Nov 23/24 period, there was again a strong reduction of the 
equatorward wind observed to the north of Kiruna during the period around that when the 
strongest upward winds occurred. 

A significant feature of the auroral activity of both these periods was that the auroral activity 
was mainly equatorward of Kiruna. 

III. Dec. 7/8 1982 

To contrast with observations which are typical of strong events where the main auroral and 
electrojet activity occurs equatorward of Kiruna, the data obtained on Dec. 7/8 1982 will be 
shown. Between 19.00 and 20.00 UT, the wind component observed in the northwest direction 
changed by about 600 m/sec, with rather smaller changes in all other directions. The vertical wind 
during this event was, however, strongly downward (50 m/sec). The main auroral activity and the 
centre of the auroral electrojet were north of Kiruna. 

This difference in location of the main auroral activity relative to Kiruna explains the different 
vertical wind behaviour which was observed caused by the two rather similar geomagnetic events. 
The comparison of the two events also defines the spatial extent of the region of strong upward 
winds rather well. 

When the major region of heating is overhead or to the north (poleward) of Kiruna, by perhaps 
200 km, the equatorward and poleward wind explosion is observed. This is associated with a 
region of downward winds overhead Kiruna, as the atmosphere relaxes on leaving the region of in- 
tense heating. When the strong heating is about 200 km equatorward of Kiruna, rather than 
directly overhead or poleward, the effects of the heating in driving upward and poleward winds 
are directly observable. 

Part of the explanation is due to the combination of the observing geometry and the geometry 
of the earth’s magnetic field. OI 630 nm is generated from regions (200 to 300 km) which are 
significantly higher than those where the particle and frictional heating maximises during auroral 
substorms (120 to 160 km). The region of maximum heating would thus be expected to be 
somewhat poleward (~50 km) of the region of maximum OI 630 nm emission, assuming that they 
were related. In practice, it might appear from some of the time-variations observed in OI 630 nm 
and the wind variations that the regions of highest heating were rather further poleward (100-200 
km) of the maximum OI 630 nm emissions. 

From the auroral intensity data obtained during these rather representative events, the latitudinal 
extent of the regions generating the strongest upward winds appeared to be about 200 km. The 
longitudinal extent is probably of the order of 1000 km, but will be related to the specific mor- 
phology of a given auroral event. 
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PERIODS OF SYSTEMATIC LARGE VERTICAL WINDS DURING MAJOR 
GEOMAGNETIC DISTURBANCES 

In Figure 5, some of the wind observations made from Kiruna during the afternoon of Dec. 17 
1982 are shown. Very fast westward winds occurred during the early part of the afternoon, peak- 
ing in excess of 900 m/sec, mainly westward and northwestward in direction. Throughout the 
period from 13:30 UT to 17:00 UT (3.5 hours), there was a mean downward vertical wind of the 
order of 50 m/ sec. Kiruna was located on the poleward side of the main part of the greatly ex- 
panded auroral oval during this extremely disturbed period. The combination of extremely rapid 
westward winds, driven 9 by ion drag in the auroral oval, and downward winds appear to be 
part of a form of standing wave, with the atmosphere relaxing as it leaves a region of intense ac- 
celeration and heating (the dusk auroral oval). 

Systematic mean vertical winds of the order of 10 to 20 m/sec have also been found 3 in 
statistical analyses of many days of thermospheric wind observations from both Skibotn (Norway) 
and Svalbard (Spitzbergen, 16 E 79 N). These mean vertical winds appear to be consistent and 
reproducible signatures of the geomagnetic forcing of the auroral and polar thermosphere, and are 
distinct from the strong time-dependent winds discussed above. 

THEORETICAL EXPLANATION AND SIMULATION OF INTENSE VERTICAL WINDS 

Figure 6 displays the simulated horizontal (vector) and vertical (scalar) wind component distribu- 
tions over the northern, winter, polar region during a period of relatively quiet geomagnetic condi- 
tions. This is a quasi-steady state simulation, the only time-dependence being due to the rotation 
of the earth, carrying the offset geomagnetic polar regions around the geographic poles, and thus 
varying the location, and solar photo-ionisation and insolation, of the geomagnetic polar regions 
as a function of UT. 

There is weak upwelling at low and middle latitudes on the dayside, and weak downwelling on 
the nightside at middle latitudes, except for a small, post-midnight, region of weak upwelling 
associated with the ‘post-midnight bulge’, and its temperature/thermospheric wind anomaly. The 
geomagnetic polar region contains regions of significant upwelling (3 to 10 m/sec) around the 
boundaries of the auroral oval, and over the central polar cap. There are also regions of weak 
downwelling in the poleward boundaries of the dusk and dawn parts of the auroral oval. 

In Figure 7, the horizontal and vertical wind distribution is shown for a second quasi-steady 
state simulation which uses a geomagnetic input within the polar regions 4 corresponding to 
moderately disturbed conditions (Kp~4). Upward wind velocities are considerably enhanced, with 
peak values of the order of 40 m/sec in the dayside cusp. Regions of moderate downward winds 
(10 m/sec) border the regions of upwelling. The peak upward winds of 40 m/sec in the dayside 
cusp are generated by a combination of convergence from the dusk and dawn auroral oval, fric- 
tional heating within the dayside auroral oval, and a soft ‘cusp’ electron source of the order of 4 
ergs/cm 2 . 

Figures 8, 9 and 10 show, with the same presentation as Figures 6 and 7, the circulation changes 
which develop during a simulated geomagnetic (auroral) substorm 4 in the region near magnetic 
midnight. ....... 

The simulated disturbance is started at 20.4 UT (Figure 8), with a large input of energy into a 
limited region of the auroral oval close to magnetic midnight. There is an immediate response in 
the vertical wind, with upward velocities of the order of 200 m/sec throughout the region where 
the peak energy input (100 ergs/cm 2 ) has been dumped. Horizontal wind changes occur with a 
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slight delay, but within 30 min (at 20.9 UT, Figure 9), there is a major expansion of the at- 
mosphere horizontally, as well as vertically, away from the region of intense energy input. The 
equatorward and poleward wind surges generated by the large input can be easily seen propagating 
away from the vicinity of the simulated substorm. Large equatorward winds have reached a 
latitude of 50 degrees in the midnight region, from their source latitude of 70 degrees, withinn 30 
min of the substorm onset. 

The major ‘substorm’ source of energy is cut off at 20.9 UT, allowing a recovery from the 
disturbance. By 21.6 UT (Figure 10), the wind distribution within the polar cap and auroral oval 
has returned to a state close to the initial state. However, at lower geomagnetic latitudes, large 
wind disturbances continue to propagate equatorward for several more hours. By 21.6 UT, the 
disturbances have reached 30 degrees latitude, and within a further hour (2.2 hours after onset), 
the propagating waves from the northern polar region meet with those generated in the southern 
hemisphere near the equator. 

Figures 9 and 10 display the extent to which these propagating waves have correlated vertical 
and horizontal wind oscillations. At a particular mid-latitude ground-based station, the vertical 
and horizontal winds will both oscillate quasi-periodically as the series of gravity waves propagate 
past the station. 

The wind disturbances set up by a single geomagnetic pulse, with a simple geometrical and time- 
dependent variation of the geomagnetic energy inputs, are quite complex, and there is also a com- 
plicated local time/longitudinal structure. The simulated disturbances correspond quite well to the 
time-dependent observations reported 1 ’ 3 , and have many of the features of the theoretical gravi- 
ty waves which have been generated in previous linear two-dimensional models. 

During extremely disturbed periods, the auroral oval expands equatorward, there is a con- 
siderable increase in the cross-polar cap electric potential, and the total deposition of geomagnetic 
energy in both polar regions may exceed 10 12 watts. The effects of such conditions on the ther- 
mosphere can be estimated from Figure 1 1 , which is taken from a simulation generated to corres- 
pond to an event observed on Dec 12 1981. This event was qualitatively similar to the one observ- 
ed on Dec 17 1982, described previously. 

Extremely large sunward horizontal winds occur in the dusk and dawn parts of the auroral oval, 
particularly in the dusk oval (near 1 km/sec), and there are also strong anti-sunward winds of 
close to 800 m/sec over much of the polar cap. The vertical winds generated under such condi- 
tions are qualitatively similar to those shown in Figure 7, but with a factor of 2 to 4 increase in 
magnitude. For example, the downward winds of 10 m/sec which occur on the poleward boun- 
dary of the dusk auroral oval are increased to about 40 to 60 m/sec. It would appear that this is 
the simplest mechanism which can account for the persistent 50 m/sec downward winds observed 
from Kiruna during the disturbance of Dec 17 1982. The persistent downward winds exist as a 
‘standing wave’, where the atmosphere is relaxing after intense acceleration and heating in the 
dusk auroral oval. 

Under less disturbed conditions, the persistent diurnal variation of vertical wind behavior 
observed at polar and auroral locations can also be explained in a similar way. There is con- 
siderable day to day variation in this behaviour, which can be related to the complex time- 
dependent variations of individual geomagnetic disturbances. 

SUMMARY 

Vertical winds of the order of 5 to 10 m/sec or larger are a persistent feature of the polar upper 
thermosphere. The signature of geomagnetic heating of the upper thermosphere under relatively 
quite geomagnetic conditions can explain average upward and downward winds of such 


301 



magnitudes which are observed in systematic studies. Such magnitudes exceed the average diurnal 
variations associated with the thermospheric ‘breathing’ in response to solar UV and EUV heating 
at low and middle latitudes. 

Violent excursions of the vertical wind, of short duration, typically 10 to 20 min, where upward 
velocities of 200 m/sec may be observed at times, and with large correlated changes of the 
horizontal wind, can be associated with intense auroral substorms and similar energetic events 
(Westward Travelling Surges, etc). Such events generate a localized explosion of the upper 
themosphere. Away from the source region, the propagating waves from such events are essential- 
ly the thermospheric gravity waves launched from an auroral source. These waves propagate 
globally, and leave the mid-latitude thermosphere ‘ringing’ for many hours, as the gravity waves 
propagating form both hemispheres interact. During periods of moderate to intense geomagnetic 
activity, the wave patterns may become so complex that it is impossible to identify individual 
disturbances by their correlated horizontal and vertical wind changes (10 to 30 m/sec at mid- 
latitudes, rather than the 100 m/sec observed at auroral latitudes). 

During extremely disturbed periods, standing waves may be generated, with embedded regions 
of large (upward and downward) vertical winds (50 m/sec), associated with, or adjacent to, 
regions of horizontal winds of the order of 500 m/ sec. Although these may be the signature of the 
deposition of intense geomagnetic heating of the thermosphere, their long duration, sometimes ex- 
ceeding 4 to 6 hours, implies that they are less efficient generators of gravity waves than the source 
associated with auroral substorms in the magnetic midnight region. 

It is not possible to identify from the FPI observations alone the combination of heating sources 
which contribute to generating the strong vertical winds. The 3-D T-D model allows various mixes 
of heat sources to be tested. If the source were kilo-volt auroral electrons, a deposition rate in ex- 
cess of 100 erg/cm 2 /sec would be required to explain 100 m/sec upward winds (with the 10 min 
observed duration). Satellite observations of precipitation include such values, but probably too 
rarely to explain the frequency of strong upward vertical winds. 

A significant or large contribution from supra-thermal or low energy electrons would decrease 
the demand on total energy input considerably, since these electrons heat the thermosphere at 
levels where the generation of large vertical winds is more efficient (per unit energy deposition 
rate). 

There are objections to this source as a sole generator of strong vertical winds. Firstly, even 
modest (10 - 20 ergs/cm 2 /sec) of 5 to 10 ev electrons require very large number fluxes of elec- 
trons, which may not be consistent with potential source regions. Secondly, there is little support- 
ing data from satellite or ground-based auroral morphology studies. In particular, the FPI itself 
tends to observe the largest vertical winds in regions which are poleward, but immediately adja- 
cent, to regions where the substorm- or surge-associated enhancement of OI 630 nm emission 

occurs. . . , , T , , 

A third contributor to the energy requirements could be intense frictional and Joule heating in 

the region of, and immediately adjacent to, the substorm- and surge activity. Intense field-aligned 
currents are known, from recent satellite observations, to complete the magnetospheric circuit 
associated with substorm and surge current systems in adjacent regions which are not strongly ex- 
cited by auroral precipitation. The strong horizontal currents required to complete the circuit of 
upward and downward field-aligned currents have to flow through regions where the ionospheric 
conductivity is not particularly enhanced. This implies that the Joule and frictional heating rates 
might be quite intense, and a valuable contributor to the total heating demanded by the observed 
strong vertical winds. Measurements of the state of the ionosphere, and its current and electric 
field systems, in regions adjacent to substorms and surges have, as yet, provided rather inadequate 
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results to fully test and complete the menu of geomagnetic heat sources contributing to these 
violent and rather fascinating wind disturbances. 
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Figure la. Thermospheric wind observations in the N, NE, E, S, W, NW and Vertical Directions 
from the FPI at Kiruna Geophysical Institute on the night of 23/24 Nov 1982. The major focus of 
interest is the correspondence between changes observed in each of the viewing directions in the 
period 18.00 to 20.00 UT before, during and after the strong geomagnetic disturbances around 
19.00 UT. Observed winds are indicated by the crosses ( + ), and the OI 630 nm emission intensity 
by the continuous line. 
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Figure 2a. Thermospheric wind observations in the N, NE, E, S, W, NW and Vertical Directions from 
the FPI at Kiruna Geophysical Institute on the night of 23/24 Nov 1982. The period 18.00 to 20.00 UT 
before, during and after the strong geomagnetic disturbances around 19.00 UT is examined in detail, 
with the instantaneous wind distribution being shown at times between 18.00 and 20.00 UT, to place 
the vertical winds in the context of the horizontal circulation. 





Figure 3a. Thermospheric wind observations in the N, NE, E, S, W, NW and Vertical Directions 
from the FPI at Kiruna Geophysical Institute on the night of 4/5 Feb 1983. The major focus of 
interest is the correspondence between changes observed in each of the viewing directions near the 
strong geomagnetic disturbance around 21.00 UT. There was a ‘great red’ aurora between 17.00 
and 18.00 UT. Observed winds are indicated by the crosses ( + ), and the OI 630 nm emission 
intensity by the continuous line. 
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Figure 4a. Thermospheric wind observations in the N, NE, E, S, W, NW and Vertical Directions 
from the FPI at Kiruna Geophysical Institute on the night of 7/8 Dec 1983. The major focus of 
interest is the correspondence between changes observed in each of the viewing directions in the 
period 18.00 to 20.00 UT before, during and after the strong geomagnetic disturbances around 
19.00 UT. Observed winds are indicated by the crosses ( + ), and the OI 630 nm emission intensity 
by the continuous line. 
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Figure 5a. Thermospheric wind observations in the West, North, 
Northwest and Vertical Directions from the FPI at Kiruna 
Geophysical Institute on the night of 17 Dec. 1982. The major focus 
of interest is the correspondence between the strong westward winds 
in the period 13.30 to 16.00 UT during the major geomagnetic 
disturbances which started at 08.00 UT. Observed winds are indicated 
by the crosses ( + ). 
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Figure 5d. Figure 5 continued. 
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Figure 6. Vertical and Horizontal wind distributions over the Northern Polar region for relatively quiet 
geomagnetic conditions. Even at a low level of geomagnetic energy and momentum input, there are upward winds 
over much of the winter auroral oval and polar cap regions. At lower latitudes, the upward winds of the dayside, 
and the downward winds of the nightside, can be observed. The magnitudes of these mid-latitude winds are less 
than 5 m/sec at 320 km. 
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Figure 7. Vertical and Horizontal wind distributions over the Northern Polar region for rather disturbed 
geomagnetic conditions. Upward winds of 40 m/sec occur in the dayside polar cusp region, where there is strong 
convergence and ‘soft’ particile heating. There are upward winds of 5 to 20 m/sec around much of the auroral 
oval. 
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Figure 8. Vertical and Horizontal wind distributions over the Northern Polar region during a geomagnetic 
substorm. The intense heating from particle and Joule heating effects in the midnight region of the auroral oval 
generate strong localized upward winds, and force an outward explosion of the thermosphere away from the 
region ot the substorm and westward travelling surge. This Figure shows the onset phase, when the rapid upward 
vertical winds have been generated, but before the outward expansion has had time to develop. 
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Figure 9. Vertical and Horizontal wind distributions over the Northern Polar region for geomagnetic substorm 
conditions. The situation, simulated in Figure 8, has developed for a further 40 min, and now the outward 
explosion from the region of the substorm is clearly defined. At this stage, the gravity waves expanding from the 
auroral source region are well defined. 
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Figure 10. Vertical and Horizontal wind distributions over the Northern Polar region for geomagnetic substorm 
conditions. This Figure shows the initial recovery phase after the disturbance. The outward explosion is no longer 
evident, and the peak upward winds have subsided. At lower latitudes, in the nightside part of the thermosphere, 
the rapid equatorward propagation of the gravity waves launched by the intense auroral disturbance can be seen. 
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Figure 1 1 . Vertical and horizontal wind distributions over the northern polar region for very disturbed geomagnetic 
conditions. This simulates the wind disturbances during a period when the IMF is extremely strong (20 nT), and 
with a large southward component. At such times, the auroral oval expands, and there may be a cross-cap poten- 
tial of 120 - 150 kV. Km/sec winds are driven in the auroral oval and over the polar cap, and regions of 
systematic upward and downward wind flows exist within, and around the boundaries of the regions of intense 
geomagnetic energy and momentum inputs. 


THE TIME-DEPENDENT DYNAMICAL RESPONSE OF THE THERMOSPHERE 
TO MAJOR GEOMAGNETIC DISTURBANCES 

David Rees 

Department of Physics and Astronomy 
University College London 
Gower Street, London WC1E 6BT, UK 

Several recent observations of thermospheric dynamics, made in the polar regions during ex- 
tremely disturbed geomagnetic periods are reviewed. In general, the magnitude and the variability 
of winds in the thermospheric polar regions increases with magnetic activity, as measured by any 
of the conventional indices. However, none of the conventional indices is a particularly good aid 
to predicting wind magnitudes. In very general terms, two major factors may be considered in 
describing the wind system. The magnitude of the IMF and, in particular, its southward compo- 
nent, determine the size of the auroral oval, and the strength of the cross-polar cap potential. This 
determines the size of the auroral oval, the magnitude of the sunward winds in the auroral oval 
and of those blowing anti-sunward over the polar cap, and is probably the major factor in deter- 
mining the rate of geomagnetic energy deposition in the thermosphere. Superimposed on this 
enhanced polar circulation system will be the effects of discrete auroral substorms. From a global 
view-point, the effect of substorms is to generate a series of strong disturbances which propagate 
from their source region, usually near magnetic midnight in the auroral oval. The energy associated! 
with discrete substorms is, however, usually a rather small proportion of the total global 
geomagnetic input during disturbed periods. Some of these substorm-related effects are discussed 
in more detail in an accompanying paper. A very large asymmetry in the development of polar 
wind disturbances can be related to the direction and magnitude of the ‘Y’ component of the 
IMF. The location of very fast anti-sunward wind jets in the polar cap follows the ion drift asym- 
metries generated by positive (dusk side maximum) or negative (dawn side maximum) values of the 
IMF Y component. In the dawn auroral oval and during moderately to very disturbed conditions, 
strong sunward winds are only observed when the IMF Y component is zero or negative. When 
the IMF Y component is positive, the winds of the dawn auroral oval are relatively small, or pure- 
ly equatorward, even during extremely disturbed conditions. These observations of thermospheric 
wind disturbances will be evaluated by comparison with global simulations of the thermospheric 
response to theoretical and semi-empirical models of the polar electric field, and of the effects of 
magnetospheric particle precipitation. 
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Figure 1. Neutral wind response observed at Kiruna on Dec 12, 1981, by a ground- 
based fabry-Perot interferometer. Zonal Wind velocity component. Winds are 
measured in locations 400 km to the east and to the west of Kiruna. Altitude of OI 
630 nm emission is about 240 km. 
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Figure 2. As Figure 1, but for the meridional wind component. Meridional components are 
measured about 400 km to the north and to the south of Kiruna. 
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Figure 3. As Figure 1, but for Dec 17, 1982, zonal wind. 
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Figure 5. As Figure 1, but for Feb 12/13, 1982, zonal wind. 
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Figure 6. As Figure 5, but for the meridional wind component. 
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Figure 9. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for relatively quiet geomagnetic 

conditions. (J2 320, 180, NUB, polar) 
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Figure 10. Thermospheric winds (320 km altitude) from the UCL 3-D I-D model lor relatively disturbed 

geomagnetic conditions. (JQ 320 180, NUB, polar) 
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Figure 11. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for geomagnetic storm 

conditions. (N4, 320, 156, NUB, polar) 
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Figure 12. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for relatively quiet geomagnetic 

conditions. (J2 320, 018, NGB, CART, 50-80) 
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Figure 13. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for relatively disturbed 
geomagnetic conditions. (JQ 320, 018, NGB, CART 50-80) 
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Figure 14. Thermospheric winds (320 km altitude) from the UCF 3-D T-D model for relatively quiet geomagnetic 

conditions. (M2 320, 018, NGB, CART, 50-80) 


18E 

1289.6 
1 254 . 4 


1219.2 
1 184.0 
1 148.8 
1 1 13.6 
1078.5 


1043.3 
1 008 . 1 
972.9 
937.7 

DEG K 


UCL 3DTD MODEL TEMPERATURE PLOT M'23200 1 8 . NGB LONG . = 


80. 

78. 


500 

M/S 


58. 

54. 

50. 


74. 


70. 


66 . 


REES, page 337 


UNIVERSAL TIME 


ELEC. FLD: LEF-1 FLs 25 to 48 of run 2 

ELEC. DENS: CHIU NO PARTICLES OR ENHANCEMENTS 

DATE: OCT 16 


18.0 


24.0 




UCL 3DTD MODEL 


Uo 

U> 

00 



LATITUDES : 
50 : 90 


REES, page 338 


500 

M/3 


.0 


JQS320000 . NUB 
TEMPERATURE PLOT 

UT= 0.0 


ELEC 


ELEC 


DATE 


SOFT 


NORThERN 

HEMISPHERE 


FLD: B2-A2 
DENS: CHTU-PIONS 
DEC 21 

PARTICLE INPUT 
IN CUSP AND POLAR CAP 
HOURS 9.6-33.6 


1290.2 

1232.0 

1204.6 

1163.6 

1 143.0 
1 125.9 

1 105.4 

1081 .4 

1054.0 

1019.8 

982.2 

DEG K 


Figure 15. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for relatively disturbed 
geomagnetic conditions. (JQ 320, 000 NUB, polar). Simulation for a positive Y component of the IMF. 
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Figure 16. Thermospheric winds (320 km altitude) from the UCL 3-D T-D model for relatively disturbed 
geomagnetic conditions. (KB2 320, 000, NUB, polar), simulation for a negative Y component of the IMF 
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It has recently been demonstrated that the dramatic effects of plasma precipitation and convec- 
tion on the composition and dynamics of the polar thermosphere and ionoshpere include a 
number of strong interactive, or feed-back, processes. To aid the evaluation of these feed-back 
processes, a joint three dimensional time-dependent global model of the earth’s thermosphere and 
ionosphere has been developed in a collaboration between University College London and Shef- 
field University. This model includes self-consistent coupling between the thermosphere and the 
ionosphere in the polar regions. Some of the major features in the polar ionosphere, which the ini- 
tial simulations indicate are due to the strong coupling of ions and neutrals in the presence of 
strong electric fields and energetic electron precipitation will be reviewed. The model is also able to 
simulate seasonal and UT variations in the polar thermosphere and ionoshperic regions which are 
due to the variations of solar photoionization in specific geomagnetic regions such as the cusp and 
polar cap. 

INTRODUCTION 

Two three-dimensional time-dependent models of the global thermosphere have been developed 
as diagnostic tools to aid the evaluation of empirical thermospheric data 1 ’ 2 . These models have in- 
cluded descriptions of the major energy and momentum inputs to the thermosphere, such as solar 
UV and EUV heating, and the geomagnetic polar inputs caused by the combined effects of 
energetic particle precipitation and the magnetospheric convective electric field imprinted on the 
polar regions 3 - 4 > 5 . 

The descriptions of the geomagnetic inputs have been somewhat improved by iterative 
ionospheric/thermospheric modelling 6 ’ 7 . It has been possible to model the polar ionosphere as in- 
fluenced by magnetospheric precipitation and solar illumination, and as modified by ionospheric 
convection and thermospheric winds which have themselves been driven by ion drag resulting from 
the combination of enhanced plasma densities and strong convection. This ionosphere, plus the 
particle precipitation energy and convection electric field has then been used to generate the 
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resulting patterns of thermospheric winds, and the consequent changes of thermospheric density, 
composition and temperature. 

The precipitation and convection models used in these initial attempts were, however, very 
idealized, and it was only possible to iterate between the ionoshperic and thermospheric models by 
a laborious data base exchange once per 24 hours UT. Fully time-dependent, and UT-dependent 
computations could not be carried out by this method, although the results represented, for the 
polar ionosphere and thermosphere, a major advance over the use of the global ‘Chiu’ 
ionosphere 8 , which has no contribution which can be recognized as due to the influence of the 
magnetospheric processes. The compositional effects of the geomagnetic input were not passed 
back to the ionospheric model, rather the MSIS empirical model was used to predice ther- 
mospheric composition under representative solar, seasonal and geomagnetic conditions. 

More complex precipitation patterns 9 and convection patterns 10 , using a one-dimensional solu- 
tion to the ionospheric effects of the precipitation, based on Roble and Rees (1977), have been 
used in an interim attempt to study the response of the polar thermosphere to severe geomagnetic 
disturbances. These attempts have been a useful means of representing the enhancement of 
ionospheric conductivity, ion drag and frictional heating both in the simulation of average 
geomagnetic input conditions 12 - 13 and for storm-time simulations 14 - 15 . Plasma transport effects 
due to convection and winds are, however, crucial at F-region altitudes. Additionally, the feedback 
of induced changes of thermospheric composition, which changes rapidly and locally during in- 
tense geomagnetic distrubances, is an important omission in previous ionospheric modelling 
studies. 

To overcome these limitations, a strategy has been developed which permits the combination of 
3-D T-D global thermospheric model and a polar/auroral ionosphere model to be run with the fre- 
quent interchange of the inter-dependent parameters. Firstly, however, the description of the 
geophysical processes which occur within the ionosphere have had to be updated, compared with 
the descriptions used in previous work 6 - 7 . 

The computational and data base strategy, permitting the frequent interchange of essential 
parameters, and the detailed physical content of the new ionospheric code will be described in 
detail in a series of papers now in preparation, however, it is appropriate to cover the major 
elements which have been changed since the prior work (i.e. 6 - 7 > 12 ). 

DESCRIPTION OF THE JOINT MODEL 

I. THE IONOSPHERIC CODE 

The polar ionospheric code developed at Sheffield University for simulating structures of the 
polar ionosphere, the plasmasphere, and the causes of complex phenomena such as the mid- 
latitude trough, has been described in a series of papers 6 - 7 > 16 . The earlier model dealt primarily 
with the problems of the upper F region, and its extension either into the plasmasphere, or into 
the high-altitude polar ionosphere. The light atomic ions, essentially 0 + , He+, and H + , were thus 
the most important components of the ionospheric code. 

In dealing with the lower parts of the ionosphere, below about 200 km, many aspects of the 
original code have had to be modified. Molecular ions, which have much more rapid recombina- 
tion rates than atomic ions, tend to dominate over the atomic ions in the increased presence of 
molecular nitrogen and molecular oxygen. As a result, residence times are thus much shorter, so 
that transport effects, particularly due to horizontal transport, become less important. However, 
strong vertical movements of ionization, resulting from fast horizontal convection, and also from 
rapid horizontal and vertical wind motions, must be considered. 
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Next, the strategy of data interchange between two vastly different styles of models has to be 
optimised: 

The thermosphere model exists on a regular global grid, using spherical coordinates in the 
horizontal dimension, and with pressure coordinates in the vertical dimension. The global grid 
rotates with the earth. 

The ionosphere model exists as a series of one-dimensional, time-dependent solutions to the con- 
ditions within a hypothetical flux tube. The contents of these flux tubes are influenced by 
precipitation, solar photoionization, and recombination and charge exchange processes (etc), as the 
flux tubes are convected through regions where the geophysical conditions vary greatly. 

A typical flux tube in the polar regions is convected throughout the dayside and nightside parts 
of the polar cap and auroral oval, experiencing highly variable solar photo-ionization, particle 
precipitation and convection velocity. The flux tube will also be convected through regions where 
the temperature and composition of the thermosphere, in addition to the wind velocity, will be 
highly variable. For convenience, the convection paths followed by individual flux tubes are used 
as the fundamental grid of the ionospheric model. 

Even when the only time-dependence in the model is that induced by the UT-dependent location 
of the geomagnetic polar regions with respect to solar illumination, each flux tube has to be 
followed for many hours. As the flux tube convects throughout the high latitude regions, the solu- 
tion will eventually become convergent and stable to the input conditions, which can vary over ex- 
treme ranges as flux tubes convect into and out of sunlight, or regions of strong convection 
and/or precipitation. 

A sequence of points thus have to be followed along each convection path, so that at any 
specific Universal Time, there is a reasonably uniform distribution of populated flux tubes along 
each of the convection paths which are used in the computations. 

II. THE THERMOSPHERIC CODE 

The UCL three dimensional time-dependent global thermospheric model has been described in 
detail in a number of recent papers 1215 . The only significant change created in the interactive and 
self-consistent modelling approach is that, in the high latitude regions, the earlier computation of 
the ionosphere from either the CHIU global model (applicable to quiet geomagnetic conditions 
only), or by the PIONS ionospheric model 1415 , is now replaced by the ionospheric computations 
described conceptually above. For the gain of realism, the additional computation involved in this 
procedure uses only marginally greater computer resources than the previous ionospheric computa- 
tion from the CHIU model, or from implementing the PIONS routine. 

PRESENT STATUS 

The present status of this work is that early diagnostic tests of the joint code have been run, to 
test the numerical method and code, and to examine the sensitivity of the resulting ionosphere to 
the range of geomagnetic and solar/seasonal conditions. 

Several of the early interesting results from these tests will be discussed, with their significance 
for the future adaption of the joint code to the routine computation of self-consistent global 
models of the thermosphere and ionosphere. 
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SIMULATIONS WITH STEADY-STATE GEOMAGNETIC INPUTS 


Some major results of the work of Quegan et al 6 , and Fuller-Rowell et al 7 will be summarized. 
The ‘Sheffield’ model, used a geometry of co-incident geographic and geomagnetic poles (UT- 
independent), and were carried out for the northern polar region at winter solstice. Separately, the 
polar ionosphere and the global thermosphere models were run to stability, and the relevant data 
sets were interchanged (ionospheric density distribution transferred to the thermosphere model, 
and thermospheric winds etc transferred to the ionospheric model). This data set interchange was 
iterated once (after 24 hours). The resulting thermospheric and ionospheric data were convergent 
after this single iteration. 

The polar ionosphere and thermosphere wind system resulting from these early computations 6 
are shown at 320 km altitude in Figure 1. For comparison, Figure 2 shows the polar ionosphere 
and the thermospheric wind system produced by the CHIU ionosphere model. 

In the polar ionosphere, the auroral enhancement of plasma density is very pronounced, and a 
‘hole’ has appeared in the polar cap as a result of plasma stagnation with this particular polar 
electric field. Around the boundaries of the auroral oval, features similar to the trough associated 
with the plasmapause have appeared 6 ’ l6 . 

In comparison with empirical data, two major factors must be considered. Firstly, although the 
values used to describe the magnetospheric precipitation are by no means extreme, individual flux 
tubes are exposed to this precipitation for many hours within both the dusk and the dawn parts of 
the auroral oval. Secondly, there is no feed-back into the thermospheric composition in this early 
simulation, so that one important consequence of intense enhancement of ion-neutral momentum 
coupling, namely the strong enhancement of molecular nitrogen density in thermospheric regions 
which are strongly heated by frictional or particle heating, does not occur due to use of the MSIS 
global model. 

The very long exposure to persistent (if moderate) particle fluxes, and the lack of thermospheric 
compositional response, generates ionospheric plasma densities within the auroral oval which are 
probably higher than would really occur. Naturally time-dependent ionospheric convection, where 
individual flux tubes are carried into and out of regions of auroral precipitation, and the enhance- 
ment of nitrogen density in regions of strong thermospheric heating, causing a major increase in 
the effective ionospheric recombination rate, are likely to suppress the high modelled auroral elec- 
tron density values. 

Features such as the ionospheric ‘polar hole’, which have been observed in satellite data, are 
due, in the model simulations, to stagnation regions, where plasma remains caught up in polar cap 
regions on the night-side where there is no precipitation and where there is no solar photo- 
ionization. 

In reality, such conditions may rarely occur, and Killeen et al 17 have recently shown that regions 
of major polar cap plasma depletions may be the signature of intense ion-neutral coupling in 
regions of very fast anti-sunward ion flow. 

TIME-DEPENDENT JOINT SIMULATIONS 

The structure of both thermospheric and ionospheric parts of the model has been created so 
that there is an easy interchange of data sets between the two models. At present, one or other 
model is being run independently, in virtually the same mode used in the original study 6 . The ex- 
ception is that thermospheric composition can be passed between the models, and a sequence of 
simulations are in progress to examine the effects of closing feedback loops such as the 
geomagnetic effects on thermospheric composition. 
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The structure of the models is organized in such a way that, with a minimum time constant of 
15 minutes, the geomagnetic inputs to the polar ionosphere and thermosphere can be varied. 
Ultimately, this will allow the realistic simulation of strongly time-dependent events. At this stage, 
the major limitation to exploitation of the joint ionospheric/thermospheric model will then be the 
limited availability of truly global data sets describing the geomagnetic input to the polar 
ionosphere and thermosphere. 

SIMULATIONS OF THE POLAR IONOSPHERE 

The first simulations using the new ionospheric code have used, as part of the input, the ther- 
mospheric wind system derived from 3D TD model simulations for approximately equivalent 
seasonal, solar activity and geomagnetic activity. Simulations have been carried out for the sum- 
mer and winter solstices and for the equinox. As a result of these simulations, it is possible to see 
the effects of variable solar photo-ionization, and of seasonally-variable thermospheric composi- 
tion in the polar regions. 

Although the simulations have been carried out for a full range of Universal Times, the results 
will only be displayed at a single UT for each of the four simulations to be described. Two of the 
simulations use the ‘A2’ polar convection field (Figure 3 10 ) to describe the northern hemisphere 
polar convection when the IMF Y component is negative. The last simulation, for northern winter 
solstice, uses the ‘B2’ polar convection field (Figure 4), corresponding to a condition when the 
IMF Y component is positive. The convections paths, followed by individual flux tubes, respon- 
ding to these convection patterns, as seen from the sun, are shown in Figures 5 and 6 respectively. 

For all of these simulations, a background thermospheric wind model, which is approximately 
that corresponding to the situation for which we are now generating the ionospheric model, has 
been used. The neutral gas composition and temperature for these initial simulations have been 
taken from MSIS 17 , preparatory to running the full interactive model. 

A. THE POLAR IONOSPHERE AT WINTER SOLSTICE: 

Northern Hemisphere, A2 convection, IMF BY negative 

Figures 7 and 8 show the polar plasma densities at 18 UT at altitudes of 160 and 320 km respec- 
tively. At each altitude, the auroral enhancement of plasma density is pronounced, as is the mid- 
latitude peak of plasma density in the early afternoon hours on co-rotating flux tubes below the 
latitude of the auroral oval. 

B. THE POLAR IONOSPHERE AT SUMMER SOLSTICE: 

Northern Hemisphere, A2 convection field (IMF BY negative) 

In Figures 9 and 10, the ionospheric structure for the same UT is shown for the summer solstice 
(N. Hemi.). In this simulation, the peak electron densities in the auroral oval are considerably 
lower than in the winter simulation. This is due to the greater amount of molecular nitrogen, at 
the equivalent altitude and pressure level, in the upper thermosphere at the summer solstice, caus- 
ing an increase in the effective ionospheric recombination rate, which more than offsets the in- 
creased solar photo-ionization in this summer period. As a result, the electron densities fall to 
about 50% of the peak values obtained in the winter solstice simulation. 
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In Figure 1 1 , the corresponding plasma density distribution from the Chiu model is shown for 
the northern summer polar region. The peak number densities are not too different, however, the 
entire distribution has changed, so that the effects of plasma precipitation are most marked, while 
the mid-latitude densities of the theoretical model are significantly below those of the Chiu model. 

DISCUSSION 

IONOSPHERIC STRUCTURES IN THE POLAR REGIONS 

The dominating structure seen in the polar plasma densities is the enhancement of electron den- 
sity in the vicinity of the auroral oval, at all altitudes, and during summer and winter solstice con- 
ditions. The plasma density enhancements are due to the prolonged exposure of individual flux 
tubes to magnetospheric precipitation as these flux tubes convect sunward within both the dusk 
and the dawn parts of the auroral oval. Although the precipitating fluxes are modest, the duration 
of exposure is prolonged (several hours). The peak plasma densities appear to be rather higher 
than occur in the real world, and we anticipate that this is due to the prolonged exposure of in- 
dividual flux tubes to magnetospheric precipitation, with inadequate feed-back due to composi- 
tional change. 

In the real world, under quiet and moderately disturbed geomagnetic conditions, the flux tubes 
may exprience much more variable convection and precipitation, which reduce the average ion 
generation. Also, thermospheric feed-back mechanisms, not yet allowed to operate fully (composi- 
tional changes), may act to deplete plasma densities on flux tubes which spend long periods within 
the auroral oval. 

SEASONAL VARIATIONS 

The seasonal variations of plasma density are in the sense which has been reported by Foster 19 
from observations with the Millstone Hill Incoherent Scatter Radar. That is, the dayside electron 
densities are higher in winter than in summer, and the peak auroral oval plasma densities are 
likewise higher in winter than in summer. In the present simulations, this is due to the higher con- 
centration of molecular species in the MSIS summer polar thermosphere. Since this is also a 
feature of the UCL 3-D T-D model 13 - 15 - 20 , this seasonal anomaly in polar plasma density will, un- 
doubtedly, also be a feature of the fully coupled ionosphere and thermosphere model. 

VARIATIONS WITH THE Y COMPONENT OF THE IMF 

The major variations which are induced by changing the polar convection field are only seen at 
the higher altitudes, where recombination is slower, and where the transport of plasma causes a 
strong modification of the plasma distribution. Figures 8 and 10 (BY negative) illustrate the pat- 
tern created when the highest anti-sunward ion flow regions is the dusk side of the polar cap. 

A tongue of enhanced plasma density reaches back across the polar cap on the dusk side when 
BY is negative. This is very reminiscent of the structures described by Weber et al 21 . They have 
reported anti-sunward convecting patches of enhanced ionospheric plasma density and related 
features of polar cap airglow features. 

When the region of highest anti-sunward flow moves to the dawn side of the polar cap (BY 
positive), the plume of enhanced plasma density follows this change in the polar convection, and 
creates a quite asymmetrical plasma distribution compared with Figures 8 and 10. 
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In the joint modelling studies which have been carried out so far, the wind structures (3-D T-D 
model) and the thermospheric composition and temperature structures (MSIS for comparable con- 
ditions) which have been imposed as input conditions to the ionospheric model are moderately 
close to the conditions we would expect to develop. This does not, however, allow the full feed- 
back processes to develop, particularly in thermospheric temperature and composition. The major 
effects which would be anticipated as a result of closing the feed-back loop are associated with the 
enhancement of recombination coefficients in regions where there is strong ion-neutral coupling 
and frictional heating, or where there is strong particle heating. These fully-coupled simulations 
will be carried out in the near future, and it will then be possible to examine the conclusions men- 
tioned earlier relating to the regions of strong depletion of polar cap plasma densities. 

A significant problem which has to be solved in the near future is how to introduce a better 
method of preventing the plasma densities within the auroral oval from reaching unrealistically 
large values. Possibly, the explanation lies in the feed-back coupling of enhancement of molecular 
neutral and thus ionic species in those regions where there is intense heating for long periods, as 
plasma is forced to remain on convection paths which remain for long periods in regions such as 
the dusk auroral oval. Another possibility is that the temporal variations of precipitation and con- 
vection ensure that, in the real world and under slightly to moderately disturbed conditions, 
specific flux tubes are continuously moving into and out of regions of significant precipitation. In 
this way, the enhancement of plasma density would be rather more modest. 

Accounting for this latter mechanism in a theoretical model is not a trivial process, but 
represents one of the outstanding goals of simulating in a realistic way the structure of the polar 
ionosphere under a range of realistic seasonal, solar and geomagnetic conditions, where the major 
feed-back mechanisms between the ionosphere and polar thermosphere are taken into account. 
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Figure 1. Thermospheric wind system at 320 km produced by the ‘Sheffield’ polar electron density model, in 
combination with the Heppner model B2 polar electric field. The polar presentation is in geographic coordinates, 
and the contours of the electron density distribution are also shown. The enhancement of electron density around 
the dusk and dawn parts of the auroral oval is distinct. 
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Figure 2. Thermospheric wind system at 320 km produced by the ‘Chiu’ global ionospheric model, with the same 
polar electric field as used to produce Figure 1 . The much lower electron densities can be seen from the scalar 
background values. The winds in the auroral oval are significantly lower, approximately half the value of those 
produced by the ‘Sheffield’ ionosphere. 


12 h 



(P 

MODEL A-2 


Figure 3. The A2 model of the polar electric field, from Heppner and Maynard 1983. This field is 
applied to the northern hemisphere when the IMF BY component is negative, and to the southern 
polar region when the IMF BY component is positive. It applied to moderately disturbed 
geomagnetic conditions 3 < kp < 4 (approx). 
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Figure 5. Convention paths traced by ionospheric flux tubes following the A2 polar electric field, as seen in 

solar and geomagnetic coordinates. 
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Figure 7. Polar plasma densities at 160 km altitude produced by the winds and electric fields (convection) due to 
the A2 electric field and the same precipitation used in the original ‘Sheffield’ study. This study has been carried 
out for conditions of (northern) winter solstice, so that there is little solar photo-ionization within the auroral oval 
and polar cap regions. It can be seen that the plasma density enhancement is strictly limited to the vicinity of the 
auroral oval and its associated precipitation. 
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Figure 8. Same as Figure 7, but at 320 km. Transport effects are now very important, and there is an intense 
plume of ionisation carried anti-sunward over the polar cap, away from the polar cusp region, by the strong anti- 
sunward ion convection velocities in the dusk region of the polar cap. 
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Figure 9. Same as Figure 7, at 160 km altitude, but for (northern) summer solstice conditions. Despite the much 
higher solar photo-ionisation source, the actual electron densities, within the auroral oval, and on the dayside, at 
sub-auroral latitudes, are lower than those generated in the winter solstice model. 
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Figure 10. Same as Figure 9, but at 320 km altitude. 
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Figure 11. Thermospheric wind system at 320 km produced by the ‘Chiu’ global ionospheric model, with the 
same polar electric field as used to produce Figure 2, but for the June solstice. The electron densities are rather 
lower, at mid-latitudes than for the Dec. solstice. 





PARTICLE PRECIPITATION INTO THE THERMOSPHERE 

(Invited Review) 

Patricia H. Reiff 

Department of Space Physics and Astronomy 
Rice University, Houston, TX 77251-1892 

Particle precipitation plays an important role in thermospheric dynamics, often being both the 
most important ionization source and the most important heat source, comparable to Joule 
heating rates in the auroral zones [Harel et al., 1981; Rees, et ah, 1983; see Reiff, 1983a] and 
typically exceeding solar ultraviolet as an ionization mechanism in the nightside auroral zones and 
winter polar caps (Figure 1, from Wallis and Budzinski, [1981]). Rees [1963] has shown that, 
roughly speaking, one electron-ion pair is produced by each 35 eV of incident electron energy 
flux; thus, over half of the incident electron energy flux goes into heating rather than into ioniza- 
tion. Precipitating ions also can produce ionization, also requiring roughly 35 eV per pair [Rees, 
1982]; however, since ion energy fluxes are typically much weaker than electron fluxes, they have 
often been neglected. 

The altitude at which the electrons and ions deposit their energy is strongly dependent on parti- 
cle energy, with ionization occurring at progressively lower altitudes for increasingly more energetic 
electrons. For 300 keV electrons, the altitude of maximum ionization occurs at 72 km; 95 km for 
40 keV; ~ 107 km for 10 keV, ~ 170 km for 1 keY, and > 216 km for 0.4 keV [Figure 2, from 
Rees, 1963]. Ions deposit their energy at altitudes higher than electrons of comparable energy 
[Rees, 1982]. The altitude of ionization becomes important in determining resultant conductivity 
profiles, with lower altitudes favoring a higher ratio of Hall to Pedersen conductance; this in turn 
affects the height distribution of Joule heating rates. 

The bulk of the particle precipitation occurs on field lines which map to the plasma sheet and 
ring current regions. The auroral zone precipitation is greatly enhanced over its typical loss cone 
fluxes by the existence of electric fields parallel to the magnetic field line [Knight, 1973]. The 
auroral electric fields are associated with upward field-aligned (Birkeland) currents. The lower- 
latitude ring current precipitation is caused principally by pitch-angle scattering in the equatorial 
plane. (See review by Reiff, [1983b]). 

It is not currently possible (nor is it likely in the near future) to have a complete, detailed, 
global distribution of particle precipitation on a nearly-continuous basis, although progress is being 
made using several complementary techniques. In-situ spacecraft and rocket particle observations 
are extremely useful and yield results on very good spatial scales along the flight track (km to m, 
for rocket flights); global models, however, must be constructed statistically, and as such lose 
much small-scale detail (Figure 3, from Spiro et al. [1982]). Using a number of spacecraft (Tiros 
and DMSP), separated in local time, one can construct a global energy flux estimate with roughly 
15-minute time resolution [David Evans, personal communication, 1984]; however, obtaining a 
two-dimensional pattern from combining different in-situ data sets is a cumbersome technique and 
cannot at present be done accurately with less than a two-hour time resolution [Simons et al., 
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1985]. Some of the statistical models have become very detailed, however, with complete spectra 
(not just energy flux and average energy) available at each local time-latitude bin (Figure 4 from 
Hardy et al., [1985]). Incoherent radar studies can also be used to infer precipitating electron 
fluxes [Vondrak and Baron, 1977], but they also are limited in field of view. 

Optical techniques show perhaps the most promise for developing an accurate two-dimensional 
pattern on reasonable time scales. Images from high altitude (such as Dynamics Explorer and the 
proposed International Solar Terrestrial Research Program) have the advantage of global coverage 
at ten-minute time resolution, at the expense of spatial resolution (-50 km) [Frank et al., 1982]. 
Lower-orbiting spacecraft, such as Hilat and DMSP (which will add ultraviolet imaging to its op- 
tical data starting in 1989), have much better spatial resolution (~km); however, a single pass 
generally cannot image the entire auroral zone, and is repeated only with 45-minute time resolu- 
tion (e.g., Meng and Huffman, 1984). Ultraviolet and X-ray imaging are particularly helpful, since 
they are responsive to the more energetic particles, and, combined with optical data, can yield a 
reasonable estimate of the energy distribution of the precipitating particles, [e.g., Luhmann, 1977; 
Imhof et al., 1982], Calibration of the optical images by “ground truth” (images which include 
radar facilities) or “space truth” (in situ particle detectors) is vital and is currently underway. 

In summary, the particle precipitation into the thermosphere is both an important ionization 
source and an important heat source; since the globally-integrated value can vary over more than a 
factor of ten, and the instantaneous local rate can vary over nearly three orders of magnitude 
(Figure 5, from Spiro et al., [1982]) global maps of precipiation rates are extremely important for 
predicting thermospheric “weather.” 
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Figure 1 . Contour plots of E p and E H (mhos) for 
the two K p cases produced by particle, 
background, and solar sources for 1700 UT on 
April 1 . (a) E p for 0 < K p < 3 0 , (b) E H for 0 < 
K p < 3 0 , (c) E p for 3 0 < K p , (d) E H for 3 0 < K p 
[from Wallis and Budzinski, 1981]. 



Figure 2. Ionization profiles for incident isotropic monoenergetic 
electrons of energies between 0.40 and 300 keV, respectively 
[from Rees, 1963]. 
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PRECIPITATING ENERGY FLUX 



Figure 3. Contours of constant precipitating 
electron energy flux sorted according to auroral 
electrojet (AE) index. The outermost contours 
correspond to a bin-average energy flux of 2.5 x 
KF 1 erg cm- 2 s -1 , with adjacent contours 
representing factor of 2 increases in energy flux 
[from Spiro et al., 1982]. 
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Figure 5. Total precipitating electron energy input (both hemispheres) per hour magnetic local 
time.(a) Data sorted according to AE index, (b) data sorted according to Kp index, and (c) 
results calculated from data given in Tables A1 and A2 of Wallis and Budzinski [1981]. 
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A METHOD FOR EXTRACTING MERIDIONAL WINDS FROM IONOSONDE 
MEASUREMENTS BY USING IONOSPHERIC MODELS 


P. G. Richards and D. G. Torr 
Center for Atmospheric and Space Sciences 
Utah State University, UMC 34 
Logan, UT 84322 


There has been great progress in modelling and measuring the dynamics of the neutral upper at- 
mosphere in recent years [Hernandez and Roble, 1984a, b]. However, future progress will depend 
on the availability of global measurements of neutral winds. 

At present, neutral winds can be deduced from the ion velocities measured by incoherent scatter 
radar and also by Fabry-Perot monitoring of the doppler shifted red line of atomic oxygen. 
Neither of the above methods can be relied on to provide adequate global coverage because radars 
are expensive to build and operate and Fabry-Perot measurements are restricted to clear skies and 
nighttime. 

The purpose of this paper is to draw attention to a relatively cheap means of supplementing the 
data base of neutral winds provided by the radar and optical measurements with data obtained by 
ionosondes. 

Rishbeth [1972] in his review of F-region dynamics, derived the relationship between the height 
of the F2 layer and the component of the neutral wind parallel to the magnetic field of the Earth. 
In this paper, we examine the sensitivity of the height and density of the F2 layer over Boulder, 
Colorado on 30 July 1982 to changes in meridional wind speed using a comprehensive in- 
terhemispheric numerical model [Richards and Torr, 1985] that solves the continuity and momen- 
tum equations for H + and 0 + , the energy equations for Te and Ti, and the 2-stream photoelec- 
tron equations to obtain electron heating rates. For the neutral atmosphere temperature and den- 
sities, we have used MSIS-83 model of Hedin [1983]. 
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Figure I. Variation in the height of the F2 layer over Boulder Colorado on 30 July 1982 (broken 
line) compared with the model calculation (full line). 
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Figure 2. Observed (broken line) and calculated (full line) electron density variations 
corresponding to the conditions of Figure 1 . 
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Figure 3. The full line represents the wind variation used to obtain the model results of Figures 1 
and 2. The broken curve is the wind variation used to obtain the model results in Figures 4 and 5. 
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Figure 4. Variation in the height of the F2 layer over Boulder Colorado on 30 July 1982 (broken 
line) compared with the model calculation (full line) using the wind variation depicted by the 
broken line in Figure 3. 
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Figure 5. Observed (broken line) and calculated (full line) electron density variations using the 
wind variation depicted by the broken line in Figure 3. 
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THE QUADRUPOLE IONOSPHERE 
H. Rishbeth 

Department of Physics, University, 

Southampton S09 5NH England 

The paper discusses the principal features that might exist in the terrestrial ionosphere if the 
geomagnetic field were to assume a quadrupole form during a magnetic polarity reversal. Two 
possible configurations are considered, the “axial quadrupole’’ (Fig. 1) and the “lateral 
quadrupole’’ (Fig. 2). Interesting phenomena are anticipated in “magnetic equatorial” regions 
where the field is horizontal, and fast magnetospherically-driven plasma convection might occur at 
latitudes where the field is steeply inclined. 

The general effect of changes of field strength on conductivity is considered; a weaker field 
raises the ionospheric conducting layer and enhances the conductivity; a stronger field lowers the 
height of the conducting layer and decreases the conductivity. 
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EISCAT VELOCITY PATTERNS 
FOR THEORETICAL PLASMA CONVECTION MODELS 

H. Rishbeth 

Rutherford Appleton Laboratory 
Chilton, Didcot, Oxon 0X11 OQX 

J. J. Sojka 

Centre for Atmospheric & Space Sciences Utah State University 
Logan, Utah 84322 

Theoretical line-of-sight velocities, as would be observed by the EISCAT radar, are computed 
for idealized models of plasma convection in the polar ionosphere. The calculations give the 
velocity as a function of range and Universal Time. For several variants of the Volland and Heelis 
convection models (Fig. 1, Fig. 2), the paper examines how the maxima, minima and reversals of 
velocity depend on beam azimuth. The analysis is designed to be applied to data from the UK- 
POLAR experiment, an example of which is shown in Fig. 3. 
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RESULTS OF VARIOUS STUDIES MADE WITH THE NCAR 
THERMOSPHERIC GENERAL CIRCULATION MODEL (TGCM) 

(Invited Review) 

R. G. Roble 

High Altitude Observatory 
National Center for Atmospheric Research 
Boulder, CO 80307 

The NCAR thermospheric general circulation model has been used for a variety of ther- 
mospheric dynamic studies. It has also been used to compare model predictions with 
measurements made from various ground-based Fabry-Perot interferometer stations, incoherent 
scatter radar stations and the Dynamics Explorer satellites. In this report the various input and 
output features of the model are described. These include the specification of solar EUV fluxes, 
and descriptions of empirical models to specify auroral particle precipitation, ion drag, and 
magnetospheric convection. Results are presented for solstice conditions giving the model pertur- 
bation temperature and circulation response to solar heating forcing alone and also with the inclu- 
sion of magnetospheric convections for two different dawn-dusk potential drops, 20 and 60 kV 
respectively. Results at two constant pressure levels Z = + 1 at 300 km and Z= -4 at 120 km are 
presented for both the winter and summer polar cap regions. The circulation over the northern 
hemisphere polar cap in both the upper and lower thermosphere are presented along with a figure 
showing that the circulation is mainly a non-divergent irrotational flow responding to ion drag. 

The results of a study made on the southern hemisphere polar cap during October 1981 where 
Dynamics Explorer satellite measurements of winds, temperature and composition are compared to 
TGCM predictions are also presented. A diagnostic package that has been developed to analyze 
the balance of forces operating in the TGCM is presented next illustrating that in the F-region ion 
drag and pressure provide the main force balance and in the E-region ion drag, pressure and the 
coriolis forces provide the main balance. The TGCM prediction for the June 10, 1983 total solar 
eclipse are next presented showing a thermospheric disturbance following the path of totality. 
Finally, results are presented giving the global circulation, temperature and composition structure 
of the thermosphere for solar minimum conditions at equinox with 60 kV magnetospheric convec- 
tion forcing at high latitudes. 
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22 JANUARY 1979 SOLAR UV FUUX 



SOLAR UV FLUX-RATIO JANUARY 1979 TO JULY 1976 



Figure 1. Solar EUV spectrum measured by H. E. Hinteregger from the Atmospheric Explorer satellite on 22 
January 1979 and ratio of the Solar EUV flux measured in January 1979 to that on July 1976. 
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Figure 2. Calculated global mean exospheric temperature profiles for various solar F10.7 values. The dashed 

lines are predictions from the MSIS model. 
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Figure 3. Mean thermospheric circulation for equinox and solstice conditions for various levels of 

geomagnetic activity. 
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Figure 5. Upper figures, the magnetospheric convection potential pattern in kV for quiet 
geomagnetic activity over the northern hemisphere polar cap for two UT’s. The lower figures 
represent the direction of the ion drift (E x B) associated with the potential pattern. Max. arrow 
represents a speed of 375 ms -1 . 
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Figure 6. Contours of constant precipitating electron energy flux measured by the Atmospheric Explorer satellite 
and sorted according to auroral electrojet (AE) index. The outermost contours correspond to a line-average 
energy flux of 2.5 x lO 1 erg cnr 2 s _1 , with adjacent contours representing factor of two increases in energy flux 
(Spiro et al., 1982) 
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Figure 7. Polar plots giving the direction and magnitude of the calculated southern hemisphere 
(summer) high-latitude circulation and contours of perturbation temperature (K) along the Z = 

+ 1 (300 km) constant-pressure surface at December solstice for (a) solar heating only; (b) solar 
heating plus magnetospheric convection with a crosstail potential of 20 kV; and (c) solar heating 
plus magnetospheric convection with a cross-tail potential of 60 kV. The wind speed associated 
with the maximum arrow is 10 m/s in (a), 200 m/s in (b), and 380 m/s in (c) (Roble et al., 1983) 
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Figure 8. Same as Figure 7, except along the Z = -4 (130 km) constant-pressure surface. The 
wind speed associated with the maximum arrow is 75 m/s in (a), 77 m/s in (b), and 89 m/s in (c) 
(Roble et al., 1983) 
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Figure 9. Same caption as in Figure 7 except for northern hemisphere (winter). The wind speed 
associated with the maximum arrow in 181 m s _1 in (a), 200 m s _1 in (b), and 339 m s _1 in (c) 
respectively. 
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(c) LOCAL TIME 

Figure 10. Same caption as in Figure 7 except along the Z = -4 (130 km) constant pressure 
surface in the Northern Hemisphere (winter). The wind speed associated with the maximum arrow 
is 108 m s-' in (a), 117 m s -1 in (b) and 120 m s _1 in (c), respectively. 
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Figure 1 1 . Meridonal cross sections giving contours of the calculated perturbation temperature (K) 
{ (a) and (b) j and meridional wind speed (ms -1 , positive northward) { (c) and (d) } , 
respectively. The local time of the meridional slice is 12.00 L.T. (a) and (c) are the perturbation 
temperature and meridional wind for the case of solar heating only, and (b) and (d) are the 
perturbation temperature and meridional wind for the case of solar heating plus magnetospheric 
convection with a cross-tail potential of 60 kV. 
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Figure 12. Same caption as Figure 11 except along the O.OOLT meridian. 
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Figure 13. Contours of perturbation temperature (K) along the Z = + 1 (300 km) surface for two 
universal times, Figure 6a, 0600 UT, and Figure 6b, 1800 UT, and along the Z = -4 (130 km) sur- 
face at 0600 UT (Figure 6c) and 1800 UT (Figure 6d) for the case where the geographic and 
geomagnetic poles are displaced, and the crosstail potential for the magnetospheric convection 
model is 60 kV. The circulation is indicated by the arrows giving direction, and the length is the 
magnitude of the winds. The length of the maximum arrow represents 395 m s _1 (Figure 6a), 

373 m s _1 (Figure 6b) 100 m s _1 (Figure 6c), and 109 m s _1 (Figure 6d), respectively. Local noon is 
indicated by an open circle on the boundary. 
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Figure 14. Contours of the calculated divergence and curl of the wind difference field shown in 
Figures 9b and 9d respectively. Positive contours represent divergence and counterclockwise 
rotation for the divergent and curl fields, respectively. Units are s -1 when these fields are 
multiplied by 10~ 5 for Figures 10a and 10b and 10~ 6 for Figures 10c and lOd. Local noon is 
indicated by an open circle on the boundary. 
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Figure 15. Vectors indicating the magnitude and direction of the ion drift, in the absence of neutral 
collisions, associated with magnetospheric convection over the southern hemisphere polar cap. The 
vectors are applied at the upper boundary of the model, and the TGCM calculates the ion drift at 
different altitudes, considering collisons with neutrals. Displaced geographic and geomagnetic poles are 
considered. The solid line through the geographic pole shows the geographic position of the satellite 
track at various universal times. 
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Figure 16. Contours of the TGCM-calculated perturbation temperature (degrees Kelvin) and wind 
vectors over the southern hemisphere polar cap at four different universal times: (a) 0300 UT, (b) 
0900 UT, (c) 1500 UT, and (d) 2100 UT. The winds and perturbation temperature are values along 
the Z = + 1 constant-pressure surface (5.4 x 10 5 /xbar or 50 Pa) in the TGCM, that is at an altitude 
of approximately 300 km. The DE 2 satellite track is shown as the solid line in the figure on a 
0900-2100 LT pass over the geographic pole. The maximum arrows represent wind velocities of (a) 
375 m s _1 , (b) 368 m s -1 , (c) 404 m s 1 , and (d) 409 m s _1 , respectively. 
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Figure 17. Contours of the TGCM-calculated atomic oxygen mass mixing ratio ^(0) along the Z 
= + 1 (300 km) constant-pressure surface over the southern hemisphere polar cap for four 
different universal times: (a) 0300 UT, (b) 0900 UT, (c) 1500 UT, and (d) 2100 UT. 
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Figure 19. Contours of temperature variation (in percent) from an average value as a function of 
universal time over the southern hemisphere polar cap. (a) Average of DE 2 FPI measurements 
during October-November 1981. (b) TGCM model predictions, (c) MSIS-83 model predictions 
along the satellite track. The average temperature over the latitude-UT grid is indicated in the 
upper right corner. The solid curve is the geographic position of the south geomagnetic pole as a 
function of UT. The TGCM model predictions include heating by a magnetospheric cusp and 
magnetospheric convection with a cross-tail potential of 60 keV. In this figure and all succeeding 
figures with a similar presentation, data in the top portion of the figure have been obtained at a 
local solar time of 0900 hours, while data in the bottom part of the figure have been obtained at a 
local solar time of 2100 hours. This corresponds to the 12-hour local time change as the DE 2 
spacecraft crosses the geographic pole in each continuous perigee pass. 
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Figure 20. Contours of the percent by which composition departs from an average value as a function of 
universal time over the southern hemisphere polar cap. (a) Average variation of atomic oxygen and (d) average 
variation of molecular nitrogen observed by the NACS instrument on board the DE 2 satellite during October- 
November 1981 for geomagnetic quiet conditions (Ap<20); (b) TGCM-calculated variation of atomic oxygen; (e) 
molecular nitrogen; (c) MSIS-83 predicted variation of atomic oxygen; and (f) molecular nitrogen for similar 
geophysical conditions. The solid curve indicates the projection of the diurnal locus of the south geomagnetic 
pole on the local time plane of the orbit. The TGCM model includes heating by a magnetospheric cusp and 
magnetospheric convection with a cross-tail potential of 60 eV. 
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Figure 21. Neutral wind vectors as a function of universal time over the southern hemisphere polar cap. 
(a) Average of DE2 FPI/WATS measurements during October-November 1981. (b) TGCM model 
predictions along the satellite track. The satellite track is from 2100 to 0900 LT (bottom to top of the 
figure), and the bulk of the arrows indicate a wind directed to the right of the satellite track. The 
reference arrow for each figure is 500 ms -1 . The solid curve indicates the projection of the diurnal 
locus of the south geomagnetic pole on the local time plane of the orbit. The TGCM model includes 
heating by a magnetospheric cusp and magnetospheric convection with a cross-tail potential of 60 keV. 
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Figure 22. The calculated (a) wind vectors and perturbation temperature 
contours along the Z = + 1 constant-pressure surface over the southern 
hemisphere polar cap at 0300 UT on October 21, 1981. The maximum 
arrow in (a) represents a wind velocity of 375 m s 1 . The Figures 3b 
through 3f give vectors of various forces acting on the neutral wind: (b) 
the ion-drag, (c) Coriolis, (d) advection, (e) pressure, and (f) viscous 
forces. Note that the length of the maximum arrow represents a 
different magnitude in each of the figures: (b) 0.62, (c) 0.05, (d) 0.12, 

(e) 0.14, and (f) 0.05 m s' 2 . 
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Figure 26. (a) Contours of the calculated eclipse function at 0300, 0400, 0500, 0600, and 0700 UT at 300 km 
indicating the temporal progression of the eclipse shadow on June 11, 1983. (b) Calculated temperature (degrees 
Kelvin) and wind (meters per second) difference fields caused by the eclipse shadow at 0300, 0500, and 0700 UT 
and along the model Z = 0 constant pressure surface at approximately 300 km. 
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Figure 27. Contours of (a) calculated temperature (degrees Kelvin) and wind 
(meters per second) difference fields, (b) vertical velocity (meters per second), 
positive upward, (c) percent deviation of molecular nitrogen mixing ratio, and (d) 
percent deviation of atomic oxygen mixing ratio at 0600 UT and along the model Z 
= 0 constant pressure surface near 300 km during June 11, 1983 solar eclipse. 

The heavy dashed line along the 140 °E longitude indicated the cross section of the 
contours presented in Figure 28. 
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Figure 28. Contours of calculated (a) temperature difference (degrees Kelvin), (b) 
vertical wind (meters per second), positive upward, (c) zonal wind (meters per 
second), positive eastward, (d) meridional wind (meters per second), positive 
northward, (e) percent change of atomic oxygen mixing ratio, and (f) percent 
change of molecular nitrogen mixing ratio along the 140 °E longitudinal slice at 
0600 UT during the June 11, 1983 solar eclipse. 
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Figure 29. The calculated global circulation, temperature and compositional structure along the Z = +2 (300 
km) constant-pressure surface for the case of solar heating plus magnetospheric convection with a cross-tail 
potential of 60 kV: (a) contours of perturbation temperature (K) and arrows giving wind direction and the length 
of the wind speed with the maximum arrow representing 336 m s~'; (b) contours of ^ 0 ; (c) contours of ^ 02 . The 
same fields along the Z = -4 (120 km) constant-pressure surface are represented in (d), (e) and (f), respectively. 
The maximum arrow in (d) represents a wind speed of 79 m/sec. 








MILLSTONE HILL MEASUREMENTS AND TGCM SIMULATION 
FOR THE 30 MAY 1984 ANNULAR SOLAR ECLIPSE 


J. E. Salah, W. L. Oliver - MIT Haystack Observatory, Westford, MA 01886 
B. A. Emery, R. G. Roble, E. C. Ridley - NCAR, Boulder, CO 80307 
H. C. Carlson - AFGL, Hanscomb Air Force Base, Bedford, MA 01731 

On 30 May 1984, the Millstone Hill incoherent scatter radar was operated to gather data on the 
effects of the annular solar eclipse on the structure and dynamics of the ionosphere and ther- 
mosphere. The eclipse path was about 3 deg. south of Millstone which experienced a maximum 
obscuration of 86% at about 1705 UT. Both the zenith and steerable antennas at Millstone were 
used in the experiment to collect data on the temporal evolution of the eclipse effects. This experi- 
ment represented the first opportunity at Millstone to collect data during an eclipse in the absence 
of a major magnetospheric disturbance which had previously made the unravelling of eclipse ef- 
fects difficult. In addition, the configuration of the experiment and analysis of the data included a 
detailed examination of the effects on the neutral atmosphere. A major catalyst for this study was 
the opportunity to compare the results with the predictions made from the Thermospheric General 
Circulation Model (TGCM) at NCAR, as a calibration point for the model. The plans for the ex- 
periment were formulated at the February 1984 Incoherent Scatter Data Users Workshop; the 
Arecibo and Sondrestrom radars also participated as part of a radar chain experiment. The 
analysis of the data from these stations is in progress; this paper presents the initial results from 
Millstone Hill. 

The observations, confirming the TGCM predictions, indicate a drop in electron density by 60% 
at the peak obscuration time in the altitude range 130 - 250 km. At higher altitudes, a smaller 
decrease (14%) is observed, in contrast to the model predictions and the difference is due to 
assumed flux boundary conditions. Although a difference in mean electron temperature exists be- 
tween model and temperature, a good agreement is found in the relative drop ( — 25%) during the 
eclipse. Within somehwat large fluctuations it is possible to discern a 50° K drop in the exospheric 
temperature measurements as predicted from the TGCM. 

From three dimensional ion drift measurements across the eclipse path a neutral wind velocity, 
governed primarily by the meridional component, of about 50 m/sec at eclipse time is observed. 
Within the uncertainties of the measurements, little change in the winds is observed prior to and 
following the eclipse, and the results are similar to previous summer day observations where a 
reversal occurs near 1400 EST. Electric fields derived from the data indicate magnitudes of the 
order of 1 mV/m or less, primarily north-westward except between 1530 and 1730 U.T. where the 
field is south-eastward. While the TGCM predicted similar overall wind magnitudes and directions 
as observed, it has been difficult to verify the convergence of wind perturbations onto the eclipse 
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path as would be expected and as indicated by the model. This is because of geometric considera- 
tions of the ion drift vectors relative to the eclipse path, and because of the sensitivity of the radar 
measurements in the F-region to both the meridional and zonal components. It is nonetheless 
possible to state that the overall dynamic behavior predicted by the TGCM is consistent with the 
observations. 
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ANNULAR SOLAR ECLIPSE OF 1984 MAY 30 



Figure 1. 
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Figure 8a. NCAR TGCM update. 
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Figure 12. 
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Figure 14. NCAR TGCM predictions of Neutral Wind Perturbations for Eclipse - 17 U.T. 
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Figure 16. 
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A DYNAMO THEORY PREDICTION FOR SOLAR CYCLE 22: SUNSPOT NUMBER, 
RADIO FLUX, EXOSPHERIC TEMPERATURE, AND TOTAL DENSITY AT 400 KM 

K. H. Schatten and A. E. Hedin 
Laboratory for Atmospheres 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 

Using the “dynamo theory” method to predict solar activity, a value for the smoothed sunspot 
number of 109 ± 20 is obtained for solar cycle 22. The predicted cycle is expected to peak near 
December, 1990 ± 1 year. Concommitantly, F 107 radio flux is expected to reach a smoothed value 
of 158 ± 18 flux units. Global mean exospheric temperature is expected to reach 1060 ± 50K 
and global total average total thermospheric density at 400 km is expected to reach 4.3 x lO" 15 gm 
cm -3 ±25%. 
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SMOOTHED R 



Figure 1. Predicted smoothed sunspot number and radio flux from 1984 to 1995 (solid curve). Estimated early 
average uncertainties (dashed curves) are smaller than in many other predictions. The previous prediction for solar 
cycle 21 (shown) compares well with the observed sunspot numbers (dots). 
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Figure 2. Smoothed maximum sunspot numbers versus polar field 
bending angle (a) and field strength and faculae count (b) at or near 
the previous sunspot minimum. The dots are labeled by the year that 
the polar field value was obtained. The arrows show the 1983 (top) 
value for the polar field bending and the 1984 (bottom) polar field 
strength. 
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Figure 3. Migration of spot zone for different R m (heights of sunspot 
maximum). Abcissae: time measured in solar rotation from sunspot 
maximum; ordinates: latitude (Waldmeier, 1939). 
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Figure 4(a). Left-The solar activity minimum to minimum duration (m-m), is graphed vs. solar maximum sunspot 
number, (b) Middle-The rising portion of the cycles’ duration (m-M) is graphed, (c) Right-The trailing end of the 
sunspot cycles’ duration (M-m). 





Figure 5. Shown are the predicted sunspot numbers from 1984 to 1995 with various authors’ methods. 
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A THREE DIMENSIONAL MODEL OF THE VENUSIAN 
THERMOSPHERE WITH SUPERROTATION 

D. R. Stevens-Rayburn 
Applied Research Corporation 
Landover, Md. 20785 

H. G. Mayr and I. Harris 
Planetary Atmospheres Branch 
Goddard Space Flight Center 
Greenbelt, Md. 20771 

An improved three-dimensional spectral model of the thermosphere of Venus is described. The 
model solves the Navier-Stokes equations and includes nonlinear effects for an arbitrary number 
of atmospheric species. A two-dimensional axisymmetric model of the superrotation of the ther- 
mosphere is also presented. This model addresses the Pioneer-Venus mission finding, which sug- 
gested the thermospheric rotation rate to be much higher than that of the planet as seen from the 
asymmetric distribution of hydrogen and helium. Both models include the effects of an anisotropic 
eddy diffusion that is consistent with atmospheric mixing length theory. 
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Figure 3. Thermospheric dynamics model results: 

(a) Equatorial cross-section of temperature versus local time at 170 km. 

(b) Equatorial cross-section of densities versus local time at 170 km. 


448 






LATITUDE LATITUDE 


(a) 


Altitude 170 km. 

HORIZONTAL WIND VELOCITY VECTOR MAP 
SPECIES OXYGEN 


VELOCITY SCALE ( 100 M/S ) 



V K l\ I V 1 /t~/i /' 

v \ \ \ ! t //. 
\ \ i t / /. 

t / / s' . 




* y 


s' / 


N. VH 


s' S H 


-GO 


-90 


/ / l \ ^ 

^ / / I \ 

'/ / 1 l \ 


y y '/ 1 1 \ 1 V \ 

^ * * « /-» “ 1 c* 



LOCAL TIME 


(b) 


VERTICAL VELOCITY ( CM / S ) 



Figure 4. Thermospheric dynamics model results: 

(a) Horizontal wind velocity vector map at 170 km. 

(b) Contour map of vertical velocity for oxygen at 170 km. 
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Figure 5. Superrotation model results: 

(a) Zonal velocity components versus altitude. 

(b) Zonal velocity versus latitude. 
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JOULE HEATING AND FIELD- ALIGNED CURRENTS: 

PRELIMINARY RESULTS FROM DE-2 

M. Sugiura 

Laboratory for Extraterrestrial Physics 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 

There are three main processes by which energy is transferred from the magnetosphere to the 
thermosphere. These are, not necessarily in the order of importance: (a) charge exchange of the 
ring current particles; (b) precipitation of charged particles; and (c) joule dissipation by the 
magnetosphere-ionosphere current systems. The importance of this last process has been recog- 
nized and the rate of joule heating has been estimated by many workers (references given in Table 
1} . 

Observations of the electric (E) and magnetic (B) fields from DE-2 are providing a new set of 
data on field-aligned currents. One of the remarkable features found in these observations is the 
high correlation between an orthogonal pair of the E and B field components (Siguira et al., 1982, 
1984; Sugiura, 1984; general references on field-aligned currents are given in these papers). Figure 
1 shows an example, where AB is the perturbation magnetic field obtained by substracting a model 
internal field from the observed field. The north-south component of the electric field is well cor- 
related with the east-west component of the magnetic field. The correlation is 0.94, using Vi sec- 
ond averages. Figure 2, which gives high-pass filtered data of the same set as Figure 1, shows that 
the correlation between the orthogonal E and B components is equally good for the fluctuating 
components; the correlation coefficient is 0.93. 

Figure 3 demonstrates this reasoning. Figure 3 shows the spacecraft X, Y, and Z components of 
AB and the spacecraft X component of E during a northern polar pass. The interplanetary 
magnetic field was northward (B z >0), and hence the size of the polar cap was small. Figure 4 
shows how the angle, a, between the orbit plane and the dipole meridian plane changes along the 
pass; the definition of the angle a is illustrated in Figure 5. The interval between 1417 and 1423 
UT was further divided into 1 -minute intervals and the correlation coefficient between AB Z and E x 
was calculated for each of these 1-minute intervals. The results are shown in Table 2. We conclude 
that the proportionality between the orthogonal B and E components is invariant with respect to 
rotation of the coordinate system. Hence the missing E z must be proportional to AB X to the same 
extent that E x is correlated to AB Z . 

Comparing AB Z with E x for the evening auroral oval pass in Figure 3, AB Z changes are very 
small, but are correlated with E x as seen in the correlation coefficient, r, for intervals D and E. 
The AB Z amplitudes are small in these intervals because the ionosphere was not illuminated and 
hence the conductivity was low. The precipitating electrons did not have energies large enough to 
produce appreciable ionization in the E region of the ionosphere (according to the LAPI plasma 
observations); hence field-aligned currents could now flow, and the electric field was large, being 
not short-circuited. The solar elevation angle along the pass is indicated in Figure 3. 


451 



Figure 6 shows a pass on the following day when IMF B z became southward. Now the size of 
the polar cap has expanded, and there are appreciable field-aligned currents in the auroral oval in 
darkness (near midnight) because of enough ionization from auroral electrons. Note the good cor- 
relation in interval C where a changed rapidly. 

From the DE-2 observations of E and AB, we conclude that the electric field and the perturba- 
tion magnetic field are orthogonal and that their magnitudes are proportional. Their directional 
relation is such that E x AB is downward. This relation between E and AB is frequently observed, 
especially on the day side. 

It can be shown that when the above relation between the electric and magnetic fields holds, the 
following equations constitute a particular solution of the current continuity equation. 

£ p -k//x o = 0 (1) 

V. J H = 0 (2) 

where £ p is the ionospheric Pedersen conductance, k the constant of proportionality between the 
orthogonal perturbation magnetic field and electric field components, fi 0 the permeability of free 
space, and J H the height-integrated Hall current. If (1) holds, then with the proportionality be- 
tween the electric and perturbation magnetic fields, we obtain a relation that the vertical compo- 
nent of Poynting’s vector equals joule dissipation, namely that the vertically downward energy 
flow equals joule dissipation; that is, the divergence of the horizontal energy flux is zero. The 
energy flux from below the ionosphere is expected to be negligible because the perturbation 
magnetic vector is parallel to the electric field; this is because the magnetic field below the 
ionisphere is mainly from the Hall current. It seems reasonable that the divergence of the Hall cur- 
rent is zero to the same extent as the plasma convection flow is horizontally divergence free. If this 
is true, then the field-aligned current sheets are parallel to the convective flow pattern or the 
equipotential surfaces. 

If we accept the above solution, then the Polynting vector Ex AB//x 0 provides the energy flux 
from the magnetosphere that is expended in the ionosphere by joule dissipation. During the 
magnetic storm of September 6, 1982 the peak energy flux reached 200 mW/m 2 in the cusp 
region. In an auroral oval pass the peak energy flux is typically a few tens of mW/m 2 under 
moderately dislturbed conditions. 

In recent years, observational data have accrued concerning the relationship between the in- 
terplanetary magnetic field and the size of the polar cap and also about the evolution of a 
substorm or a magnetic storm. It is suggested that these findings be incorporated in future model 
calculations. 
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Table 1 



DAY 81251 

INTERVAL 

UT 

1417- 1418 

1418- 1419 

1419- 1420 

1420- 1421 

1421- 1422 

1422- 1423 


8 SEPTEMBER 1981 


ANGLE a 

MLT 

CORR. 

2P 

DEGREES 

HRS 

COEFF. 

MHOS 

151-137 

1 3.4-1 4.4 

0.99 

11.3 

137-109 

14.4-16.5 

0.98 

7.4 

109-68 

16.5-19.3 

0.99 

6.0 

68-40 

19.3-21.0 

0.90 

2.8 

40-27 

21.0-21.8 

0.97 

4.1 

27-20 

21.8-22.2 

0.97 

2.7 


Table 2 



DE - 2 


OCT 1981 



ALT (KM) 358 347 

M LT (HR) 8.5 8.8 

HAT (DEG) 81.2 78.4 


338 

9.0 

75.1 


331 

9.1 

71.5 


324 

9.2 

67.6 


319 

9.2 

63.7 


Figure 1 . The north-south component of E, the north-south and east-west components (in dipole 

coordinates) of the peturbation B field. 
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DE- 2 FILTERED 



UT 02:51 02:52 02=53 

Figure 2. Data in Figure 1 were passed through a hi 
components. The orthogonal components of 1 
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DE-2 8 SEPTEMBER 1981 
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ILAT(DEG) 48.5 64.2 79.2 81.7 67.2 51.2 

ML AT (DEG) 51.1 67.0 81.9 78.3 62.8 46.9 

Figure 3. A northern polar pass of DE-2 on 8 September 1981 during IMF B z > 0. The 
correlation between the orthogonal pair of B and E components (i.e., AB Z vs. E x ) is high 
even when the angle a (Fig. 4) changes rapidly as in intervals A and B. 
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Figure 4. Angle a along the pass in Fig. 3. 
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Figure 5. Angle, a, between the orbit plane and the dipole meridian plane; a is 
measured from the dipole south direction to the direction of the spacecraft velocity. 
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Figure 6. A northern pass when IMF B z < on September 1981. The angle a changed 
rapidly from intervals B to C. The field-aligned current was appreciable on the nightside 
because of auroral electrons. 
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OBSERVATIONS OF COMPOSITION FROM PIONEER VENUS 

H. A. Taylor, Jr. 

NASA/Goddard Space Flight Center 
Laboratory for Atmospheres 
Greenbelt, Maryland 20771 

Low latitude distributions of atmospheric neutral hydrogen have been derived at Venus for the 
period 1979-1980. In-situ measurements of H + , 0 + , O, and C0 2 obtained from the ion and 
neutral mass spectrometers on the Pioneer Venus orbiter are combined with the appropriate 
chemical equilibrium relationship to determine the abundance of neutral hydrogen which is very 
difficult to measure directly. The measurements are all obtained below 165 km on the nightside 
and below 200 km on the dayside, based on evidence for chemical equilibrium prevailing up to 
those altitudes. During the period examined nearly three complete diurnal cycles were available 
and a comparison of the year-to-year variation in hydrogen content is made across the dawn 
region where the distributions of light gases are most pronounced. The dawn bulge in H (and also 
in He) which was reported from the first diurnal cycle by Brinton et al. (1979) is found to persist. 
The hydrogen concentrations peak near 2-5 x 10 7 /cm 3 , which is about 400 times greater than 
dayside concentrations. Superimposed upon the diurnal variation are strong day-to-day variations 
in which n(H) changes by as much as a factor of five. Such variations are linked to pronounced 
changes in the ion and neutral composition which sometimes occur in association with solar wind 
disturbances passing the planet. The interaction of the solar wind and the planetary environment 
somehow results in large changes in the relative abundances and scale heights of the ion and 
neutral species, thus modifying the derived values of n(H). These variations in the ion distributions 
are not surprising owing to the strong dependence of the nightside ionization upon convection 
from the dayside and associated sensitivity of this convection to changes in solar wind pressure 
and interplanetary magnetic field variations. The variation exhibited by the neutrals, however, ap- 
pears to require some other explanation owing to the limited momemtum transfer between the ions 
and neutrals. Allowing for these short term perturbations, there appears to be no clear evidence 
for interannual variation in n(H) during the period examined, apparently consistent with the very 
small change in solar EUV flux over the same interval. 
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Figure 1 . Local time distributions of n(H) derived from in-situ ion and neutral composition 
measurements at Venus during 1979-1980. Clusters of points are for individual orbits at specific 
local times. Gaps occur due to data interruptions and mode incompatibilities. The dawn bulge in 
hydrogen is evident featuring a sharp night-to-day gradient near 0500-0600. Sharp day-to-day 
variations in n(H) are attributed to solar wind perturbations and associated interaction with the 
ion and neutral constituents used in the derivation of hydrogen. 
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Figure 2. A comparison of the individual neutral and ion species variations observed in association 
with the passage of a strong solar wind disturbance on March 1 relative to a quiet day on March 
4. Note for the neutrals in the upper panel that each of the species is more abundant on the 
disturbed day while in the lower panel the ion abundances are lower on the disturbed day. The 
reduction in nightside ion abundances with increased solar wind pressure occurring on March 1 is 
understandable on the basis of the day-to-night convection necessary for maintenance of the 
nightside ionosphere. The response of the neutrals is unexplained but apparently reflects 
compound results of energetic and dynamic perturbations from undetermined sources included in 
the overall solar wind disturbance event. 
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WAMDII OBSERVATION OF AN AURORAL ATMOSPHERIC WAVE EVENT 

R. H. Wiens and G. G. Shepherd 
C.R.E.S.S. York University, Toronto, Canada M3J 1P3 

K. V. Paulson 

Dept, of Physics, University of Saskatchewan, Saskatoon, Canada S7N OWO 

Field tests of WAMDII 1 (Wide Angle Michelson Doppler Imaging Interferometer) on February 
23/24, 1984 in Saskatoon produced wind images of the aurora that show well-defined wave struc- 
tures in X5577. Wavelengths and velocities of these structures suggest an acoustic-gravity wave in- 
terpretation, but their short duration warrants their being termed an event. The dopplergrams are 
presented with emphasis on the interpretation of the waves and a discussion of possible sources. 

'G. G. Shepherd and WAMDII Coinvestigators, submitted to Geophys . Res . Let., 1984 
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